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RESUMEN 

Se presentan resultados, interpretaciones y recomendaciones correspondientes al 
ensaye de diez trabes acarteladas de concreto reforzado de tamano casi real, disenadas para 
fallar a cortante ante carga estatica. Los resultados experimentales muestran que estos 
elementos presentan un mecanismo resistente a cortante que difiere con el observado en vigas 
prismaticas, debido a un evidente trabajo de arco propiciado por la geometria de las trabes. 
Asimismo, se presenta una comparacion entre dos metodos analiticos simples que permiten 
estimar las deflexiones esperadas en estos elementos para las condiciones del primer 
agrietamiento diagonal y falla por cortante. 



SUMMARY 

This paper presents the final research results and interpretations about the testing of ten 
prototype reinforced concrete haunched beams. Beams were designed to present a shear 
failure under static loading. Results show that these elements have a different behavior with 
respect to prismatic beams, because haunched beams favor an arching action in the haunched 
length as a main resisting failure mechanism. In addition, two simple analytical methods to 
compute deflections in haunched beams at the first diagonal cracking and shear failure are 
compared in this study. 
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1 INTRODUCCION 

El uso de trabes acarteladas de concreto reforzado, TACR, es atractivo como 
alternativa estructural para el diseno de edificios de dimensiones importantes, ya que 
proporcionan algunas ventajas sobre las trabes de concreto de seccion constante, entre las 
cuales se pueden mencionar las siguientes (Tena-Colunga, 1994): 

• Aumentar la rigidez lateral. 

• Reducir el peso. 

• Reducir la altura de entrepiso. 

• Disminuir la cantidad de concreto. 

• Facilitar la ubicacion de las diferentes instalaciones. 

La principal desventaja de estos elementos es que su habilitado y cimbrado es mas 
laborioso respecto a las trabes de seccion constante. 




lc « L/3 1 = (L-(21c)) lc « L/3 



L 

Figura 1 Geometrfas propuestas en Mexico para TACR: a) hoy dia, b) inicios del siglo XX. 

El uso de las TACR en edificios ubicados en las zonas de alto peligro sismico de la 
ciudad de Mexico, no es nuevo. De hecho, fueron muy utilizadas desde la primera mitad del 
siglo pasado. En fechas recientes se han vuelto a proponer estos elementos, los cuales, en 
cuanto a geometria, poseen algunas variaciones respecto a las utilizadas en la primera mitad 
del siglo pasado. Basicamente estas diferencias son dos: 

1. En los edificios de reciente construccion, se han propuesto trabes en las cuales las 
cartelas tiene aproximadamente un tercio de la longitud efectiva del claro, situacion 
claramente diferente a las trabes de inicios del siglo pasado, en las cuales la 
longitud de las cartelas era muy inferior a un tercio del claro. 

2. Los angulos de acartelamiento observados en edificios de reciente construccion, 
son del orden de cinco a ocho grados. En los edificios de principios del siglo 
pasado, los angulos de acartelamiento eran muy superiores a los observados 
recientemente. 
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La figura 1 muestra las diferencias geometricas entre ambos tipos de TACR; es 
evidente que los acartelamientos propuestos en los inicios del siglo pasado, semejan mas a 
unas mensulas en los extremos de las trabes. 





Figura 2 Edificios de reciente construccion ubicados en la ciudad de Mexico que utilizan TACR. 

Existen libros especializados que proporcionan ecuaciones teoricas para estimar la 
contribucion del concreto, V cR , en la resistencia a cortante de las TACR. Estas ecuaciones 
estan basadas en los planteamientos clasicos de trabes de seccion constante (prismaticas), 
pero modificados para considerar la variacion del peralte efectivo en la longitud de las cartelas 
(Park y Paulay, 1997; MacGregor, 1997; Nilson, 1999). Estas ecuaciones deben validarse con 
evidencia experimental, ya que como senala Nilson en su libro “existe muy poca investigation 
referente a la resistencia a cortante de vigas con altura variable” (Nilson, 1999). 

La poca evidencia experimental del comportamiento a cortante estatico de las TACR 
con proporciones similares a la de la figura 1-a, se puede resumir en dos trabajos 
fundamentales: 1) trabes con refuerzo por cortante (Debaiky y El-Niema, 1982) y, 2) trabes sin 
refuerzo por cortante (MacLeod y Houmsi, 1994). Ambos trabajos concluyen que la geometria 
de las TACR modifica el comportamiento a cortante respecto del ampliamente reportado en las 
trabes de seccion constante. 
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Como consecuencia de la escasa informacion disponible sobre el comportamiento de 
las TACR, los reglamentos vigentes de diseno de estructuras de concreto mas utilizados en 
Mexico, como son las recientes Normas Tecnicas Complementarias para el Diseno de 
Estructuras de Concreto version 2004 (NTCC-04), y el reglamento del Instituto Americano del 
Concreto version 2002 (ACI-318-02), no proporcionan recomendaciones especfficas para 
diseno de trabes de seccion variable. De esta manera, se podrfa afirmar que el diseno de estos 
elementos ha dependido fundamentalmente de la experiencia y juicio del calculista estructural. 

Ante la evidencia de que las TACR son utilizadas en la practica mexicana como una 
solucion estructural, y aunado a la escasa informacion existente para disenarlas correctamente, 
se tomo la decision de desarrollar un programa experimental que permita comprender su 
comportamiento a flexion y cortante ante cargas estaticas y reversibles, para determinar cuales 
son las variables mas significativas a considerar en su diseno. Como primera etapa de este 
proyecto experimental, se decidio iniciar el estudio de su comportamiento a cortante estatico, 
tomando como base los experimentos realizados previamente (Debaiky y El-Niema, 1982; 
MacLeod y Houmsi, 1994), pero utilizando la geometrfa y armados observados en los edificios 
de reciente construccion en la ciudad de Mexico (figura 2). El estudio incluye la comparacion de 
las deflexiones medidas en los experimentos, contra las que se obtienen al utilizar dos metodos 
elasticos que consideran de manera sencilla la condicion agrietada de los elementos. 



2 DESCRIPCION DE LOS ESPECIMENES DE PRUEBA 
2.1 Geometri'a y cargas 

La geometrfa de las TACR se determino segun lo observado en edificios de reciente 
construccion en la ciudad de Mexico, las cuales concuerdan con las recomendaciones de 
Mendez (1991) para proporcionar TACR con cartelas lineales. El claro efectivo de las 
trabes fue L= 2.80 m, las cartelas se colocaron en los tercios extremos del claro efectivo 
(/c= 93.3 cm). Los angulos de acartelamiento, a , fueron 0° (prismatica), 3.07°, 6.12°, 9.1 3° y 
12.10°. El acartelamiento se obtuvo al variar el peralte total del tercio central de las trabes 
(h min = 45, 40, 35, 30 y 25 cm), y manteniendo constante el peralte total de los extremos 
{h m ax = 45 cm). La geometrfa de todas las trabes es tal que se garantiza la condicion de viga 
esbelta que consideran las NTCC-04 en la seccion 2.5.1. 1, ya que L/h max > 5.0. Ademas, 
con la finalidad de no magnificar el mecanismo de arco caracterfstico en las TACR 
(Debaiky y El-Niema, 1982; MacLeod y Houmsi, 1994), la relacion entre claro de cortante, 
a, al peralte efectivo maximo fue mayor a 2.50 (Wang y Salmon, 1979). El ancho es 
constante en todos los especimenes y fue de 22 cm. Las trabes se ensayaron simplemente 
apoyadas y las cargas monofonicas, V, se aplicaron a 1 0cm de los vertice que se forma en 
la interseccion de los tramos acartelados con el tramo prismatico, hacia el centro del claro 
(en lo sucesivo vertices). La figura 3 muestra la geometrfa general de los especimenes de 
ensaye y la ubicacion de las cargas. 
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Figura 3 Geometrfa y cargas de los especimenes de prueba. 

2.2 Diseno de los especimenes de prueba 

Para cada uno de los cinco acartelamientos se construyeron dos elementos: a) uno sin 
refuerzo transversal, solo con cuatro estribos para sujetar el armado longitudinal, y b) otro con 
un armado minimo en el claro de cortante, igual al establecido en la seccion 2. 5. 2. 2 de las 
NTCC-04 para trabes prismaticas. Las propiedades indice de los materiales fueron, para el 
concreto, un esfuerzo a compresion r c = 250 kg/cm 2 , y para el acero de refuerzo un esfuerzo de 
fluencia fy= 4200 kg/cm 2 . La contribucion a cortante del concreto, V cR , se calculo segun las 
recomendaciones de la seccion 2. 5. 1.1 de las NTCC-04, para trabes de seccion constante sin 
preesfuerzo provistas con una cuantia de acero a tension mayor o igual a 0.015 (p >0.01 5): 



= 0.5 bdfT' 



(D 



donde b es el ancho de las trabes, y d el peralte efectivo de la seccion. Para cumplir el armado 
minimo se colocaron en el claro de cortante estribos de dos ramas del numero 2.5 (0=0.79 cm), 
separados cada 18.5 cm. La contribucion del acero de refuerzo transversal, V sR , tambien se 
calculo segun las recomendaciones del las NTCC-04 (seccion 2. 5. 2. 3), para estribos 
perpendiculares al eje del elemento: 



V. 



sR 




( 2 ) 



5 

donde A v es el area de la ramas de los estribos, d el peralte efectivo de la seccion, y s la 
separacion entre estribos. De esta manera, la resistencia nominal a cortante de las trabes es la 
suma de las contribuciones del concreto y del acero: 



V„=V cR +V sR 



( 3 ) 



La resistencia nominal a cortante se estimo suponiendo dos peraltes efectivos crfticos 
en la longitud de las cartelas: 1) el que se ubica a una distancia de un peralte total del apoyo, y 
2) el peralte minimo. La figura 4 muestra esquematicamente los peraltes seleccionados. 
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A un peralte Peralte 

del apoyo mmimo 




Figura 4 Peraltes seleccionados en el calculo del cortante nominal. 

Una vez calculada la resistencia nominal a cortante de los especimenes, se calcularon 
las areas de acero longitudinal para garantizar la falla por cortante. La resistencia a flexion se 
calculo utilizando el nuevo bloque equivalente de compresion del concreto de las NTCC-04, el 
cual es similar al del reglamento ACI-318-02. En el calculo de la resistencia a flexion se utilizo el 
peralte efectivo del tramo central prismatico {d min ), ya que allf se presentan las maximas 
demandas de momento. Con la finalidad de no introducir variables extras, los diez especimenes 
cuentan con el mismo armado longitudinal. En todos los calculos se considero un recubrimiento 
de 4cm y factores de reduccion unitarios, por lo que, como se explico anteriormente, las 
resistencias obtenidas fueron nominales y no de diseno. El refuerzo longitudinal propuesto para 
los especimenes fue de cuatro varillas del numero 8 (^=2. 54cm), en el lecho inferior 
(A s =20.26cm 2 ), y tres varillas del mismo diametro en el lecho superior (A s '=1 5.20cm 2 ). A pesar 
de las altas cuantfas de acero provistas en las trabes, teoricamente cualquier seccion de estas 
cumple el criterio de seccion subreforzada para trabes doblemente armadas (Park y Paulay, 
1997). El criptograma de identificacion de las trabes es: TASCa,-R y , donde el contador / indica el 
angulo de acartelamiento a: 0=0°, 1=3.07°, 2=6.12°, 3=9.13°y 4=12.10°. El contador y identifica 
el tipo de refuerzo transversal por cortante: y=0 indica ausencia de refuerzo transversal, y y=1 es 
para los elementos con el refuerzo mmimo. El armado de cada uno de los especimenes de 
ensaye se condensa en la tabla 1 ; dicha tabla se complementa con los armados mostrados en 
las figuras 5, 6 y 7. 



Tabla 1 Tipificacion de los armados de los especimenes de prueba. 


Trabe 


a 


Refuerzo flexion 


Refuerzo cortante 


Superior 


Inferior 


Cartelas 


Centro 


Vertice 


TASCaO-RO 


0° 


3#8 


4#8 


1 E#2.5 (armado) 


1 E#2.5 (armado) 


- 


TASCal-RO 


3.07° 




it 


1 E#2.5 (armado) 


1 E#2.5 (armado) 


- 


TASCa2-R0 


6.12° 




it 


1 E#2.5 (armado) 


1 E#2.5 (armado) 


- 


TASCa3-R0 


9.13° 




it 


1 E#2.5 (armado) 


1 E#2.5 (armado) 


- 


TASCa4-R0 


12.10° 




it 


1 E#2.5 (armado) 


1 E#2.5 (armado) 


- 


TASCaO-RI 


0° 




it 


7E#2.5 @ 18.5 cm 


1 E#2.5 @ 18.5 cm 


3E#2.5 @ 18.5 cm 


TASCal-RI 


3.07° 




it 


7E#2.5 @ 18.5 cm 


1 E#2.5 @ 18.5 cm 


3E#2.5 @ 18.5 cm 


TASCa2-R1 


6.12° 




it 


7E#2.5 @ 18.5 cm 


1 E#2.5 @ 18.5 cm 


3E#2.5 @ 14 cm 


TASCa3-R1 


9.13° 




it 


7E#2.5 @ 18.5 cm 


1 E#2.5 @ 18.5 cm 


3E#2.5 @ 7.5 cm 


TASCa4-R1 


12.10° 




it 


7E#2.5 @ 18.5 cm 


1 E#2.5 (armado) 


3E#2.5 @ 4.5 cm 







Figura 6 Armado del elemento TASCa4-R1. a= 12.10°. 



A los especimenes con refuerzo transversal se les coloco un refuerzo adicional en la 
zona del vertice, con la finalidad de evitar una falla local debida a la fuerza vertical originada por 
el cambio brusco de direccion del acero longitudinal inferior (Park y Paulay, 1997). Este 
detallado se puede proponer de muchas maneras, en este caso se supuso que la fuerza vertical 
esta asociada a la fluencia del refuerzo inclinado, y que esta se concentra en el vertice. De esta 
manera, se coloco un estribo a cada lado del estribo de vertice, a una separacion s cd , calculada 
segun la ecuacion 2 (figura 8). 




SECCION 1 SECCION 2 



Figura 7 Secciones transversales tipicas. 
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Figura 8 Refuerzo en la vecindad del vertice en los elementos reforzados por code. 

Los elementos se instrumentaron internamente con deformimetros electricos (strain 
gages), de tal manera que se pudiera cuantificar la tension y compresion del acero longitudinal 
en la longitud acartelada, asi como la contribucion del acero de refuerzo transversal de los 
elementos correspondientes. La instrumentacion externa solamente se diseno para medir la 
deflexion vertical en ambos vertices y en el centro del claro. 



3 RESULTADOS EXPERIMENT ALES 

Se identificaron durante los ensayes tres cortantes caracterfsticos: 1) el que produjo el 
primer agrietamiento diagonal, V agr , 2) el cortante de falla, V u , y 3) el cortante asociado al 
colapso, V c i ps . Los cortantes asociados al primer agrietamiento y falla se verificaron 
analfticamente mediante el cambio de rigidez secante (AK) de la curva cortante aplicado contra 
desplazamiento al centro del claro, entre cada incremento de cortante /, segun la ecuacion 4. 
En la figura 9 se muestra la grafica AK, contra el cortante aplicado V, para el elemento TASCa2- 
RO. La razon por la cual el cortante de colapso no se identified como el cortante de falla, es que 
para ese estado de carga, el patron de dano no corresponde al de un estado limite de falla 
(figura 10), si no que mas bien, lo que se pretendio determinar fue la capacidad remanente de 
carga y/o deformacion de las TACR, por lo que una vez alcanzada la falla, se aplico de nuevo la 
carga a hasta el colapso de los especimenes. 

V. -V. , 

AK, = ' i > 1 (4) 




V {Ton) 



Vanacion de rigidez secante en la trabe TASCa2 — R0 

Figura 9 Comprobacion analftica de los cortantes del primer agrietamiento y falla. 
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a) b) 

Figura 10 Tipificacion del dano en las TACR: a) cortante ultimo, y b) colapso. Elemento TASCa3-R1. 

0=9.13°. 



En la figura 1 1 se muestran las graficas V-S para los elementos TASCa2-R0 y TASCa3- 
R1 . En estas graficas se observa que se presento un comportamiento lineal en las TACR hasta 
el instante del primer agrietamiento diagonal. De la figura 11 tambien es evidente que las TACR 
poseen una notoria capacidad de deformacion sosteniendo la carga de falla, situacion 
claramente distinta a las trabes prismaticas, donde una vez ocurrida la grieta diagonal de falla, 
la trabe pierde toda capacidad de sostener la carga. 

3.1 Agrietamiento 

El agrietamiento presentado en las TACR confirmo el mecanismo de arco reportado por 
Debaiky y El-Niema (1982) y MacLeod y Houmsi (1994). Este mecanismo se manifesto por la 
presencia, en cada cartela, de puntales inclinados de compresion (figura 12). Los puntales 
tienden a formarse entre el punto de aplicacion de la carga y la longitud media de las cartelas. 
Obviamente, la demanda fue mayor en los puntales de las trabes con refuerzo por cortante, tal 
y como se observa en la figura 12-b, donde la reaccion del puntal deformo considerablemente 
la geometrta inicial de la trabe. Ademas, sin importar la presencia de los estribos, se presento 
siempre mas de una grieta diagonal de importancia y el agrietamiento se distribuyo en casi todo 
el claro de cortante. 



g 





a) b) 

Figura 11 Curvas tfpicas V-S : a) sin estribos (a=6.12°), y b) con estribos (a=9.13°). 
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a) b) 

Figura 12 Dano presentado en las TACR: a) elemento TASCa2-R0 (a=6.12°), y b) Elemento TASCa2- 

R1 (a=6.1 2°). 

Con base en el patron de agrietamiento observado, se concluyo que el detallado en el 
vertice de las trabes con refuerzo por code fue adecuado, ya que el agrietamiento no se 
propago hacia el centro de la trabe, como fue el caso de las trabes sin refuerzo transversal. 
Asimismo, en ningun elemento de la serie R1 se presento una falla local en el vertice. 

3.2 Resistencia a cortante 

La tabla 2 compara la resistencia a cortante experimental, V u , contra la resistencia 
nominal, V n , calculada a partir de las resistencia medida del concreto en cada trabe y del acero 
de refuerzo transversal (4592 kg/cm 2 ). La comparacion se realiza para los dos peraltes 
seleccionados en la etapa de diseno. Tambien se indica el cortante del primer agrietamiento 
diagonal, V agr . La resistencia del concreto es la que corresponde a cilindros de control 
ensayados en un margen de mas o menos dos dfas del ensaye de la trabe correspondiente. 

De la tabla 2 se observa que en los elementos sin refuerzo transversal (R0), se 
obtuvieron mejores predicciones con el peralte minimo que con el otro criterio. En el caso de los 
elementos con refuerzo transversal (R1), se observa una tendencia dividida, ya que en los 
elementos TASCal-RI (3.07°) y TASCa2-R1 (6.12°) el peralte minimo proporciono la mejor 
estimacion, mientras que para los elementos TASCa3-R1 (9.13°) y TASCa4-R1 (12.10°), el 
peralte ubicado a un peralte total del apoyo proporciono las mejores predicciones. Un aspecto 
de suma importancia es que al disenar por cortante las TACR bajo la optica de un diseno por 
secciones, el utilizar el peralte minimo no necesariamente proporciona un criterio conservador, 
ya que, como se detallara mas adelante, se requiere considerar la fuerza vertical que introduce 
el refuerzo longitudinal inclinado. 
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De la misma tabla 2, es evidente que la resistencia ultima a cortante es funcion del 
angulo de acartelamiento, ya que al aumentar el angulo de las cartelas, el volumen de concreto 
disminuye, lo que implica que la resistencia a cortante tambien disminuya. En promedio las 
TACR exhibieron un 60% de la resistencia a cortante de las trabes prismaticas de referenda. 



Tabla 2. Comparacion entre los cortantes experimentalmente y nominales. 


Elemento 


a 


f'c 

(kg/cm 2 ) 


EXPERIMENTALES (ton) 


NOMINALES (ton) 


Vagr 


Vu 


A un peralte total 
del apoyo 


Peralte 

minimo 


TASCaO-RO 


0° 


334 


4.50 


7.50 


8.24 


8.24 


TASCal-RO 


3.07° 


321 


5.75 


6.75 


7.60 


7.10 


TASCa2-R0 


6.12° 


295 


5.00 


6.00 


6.83 


5.86 


TASCa3-R0 


9.13° 


236 


2.75 


3.75 


5.71 


4.40 


TASCa4-R0 


12.10° 


281 


2.50 


3.00 


5.78 


3.87 


TASCaO-RI 


0° 


315 


7.00 


25 


17.98 


17.98 


TASCal-RI 


3.07° 


269 


11.00 


20 


16.34 


15.25 


TASCa2-R1 


6.12° 


292 


8.75 


17 


15.60 


13.37 


TASCa3-R1 


9.13° 


288 


4.00 


12 


14.51 


11.18 


TASCa4-R1 


12.10° 


212 


4.00 


8 


12.64 


8.47 



3.3 Capacidad de deformacion 

Aunque las TACR presentaron resistencias ultimas a cortante inferiores a las trabes 
prismaticas de referenda, las trabes de seccion variable presentaron una mayor capacidad de 
deformacion (en promedio un 35% mayor), en la carga ultima, lo cual ya habia sido reportado 
de manera cualitativa por MacLeod y Houmsi, (1994). Este aumento se relaciona 
principalmente con la capacidad de redistribuir el agrietamiento en la longitud acartelada 
(figuras 10 y 12), mas que en la reduccion de la inercia en la zona acartelada. La mejor 
distribucion del agrietamiento permitio que el cortante ultimo no ocurriera ante la aparicion de la 
primera grieta diagonal importante, como sucedio con las trabes de seccion constante. Sin duda 
este comportamiento es benefico, ya que segun se constato en los ensayes, las TACR 
presentan fallas por cortante marcadamente menos subitas, en comparacion con las trabes 
prismaticas. 

La figura 13 muestra el comportamiento de las TACR para un mismo angulo de 
acartelamiento con y sin estribos. Como era de esperarse, se observa una respuesta similar 
hasta el primer agrietamiento diagonal de las trabes sin refuerzo transversal, y un notorio 
aumento en la capacidad de deformacion total en las trabes con refuerzo transversal. 
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Deflexiones medidas en las trabes TASCoil— Rj Deflexiones medidas en las trabes TASCoi4— Rj 

a) b) 

Figura 13 Comparacion de la respuesta en TACR: a) sin estribos («=3.07°), y b) con estribos 

(a=12.10°). 



3.4 Refuerzo longitudinal 

Las lecturas de los deformimetros electricos colocados en el acero longitudinal 
inclinado (tension), mostraron que las deformaciones asociadas a la falla en los especimenes 
sin refuerzo transversal, fueron mucho menores a la deformacion de fluencia, ya que se 
observo en promedio un 10% de esta, segun los ensayes a tension de las probetas de control. 





Deformaciones en acero de refuerzo. Trabe TA5Co!2-R1 

a) b) 

Figura 14 Comportamiento tfpico del acero longitudinal en las TACR: a) sin estribos, y b) con estribos. 

Por otra parte, las deformaciones en las trabes con refuerzo transversal mostraron dos 
tendencias marcadas: 1) deformaciones del orden del 70% de la fluencia en la trabe prismatica 
y en la de menor acartelamiento (a=3.07°), y 2) deformaciones iguales o mayores a la fluencia 
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en los especimenes restantes. Por otra parte, en ningun caso se observaron deformaciones 
significativas del acero en compresion. En la figura 14 se muestra el comportamiento tipico 
observado en el acero longitudinal de las TACR. 

3.5 Refuerzo transversal 

El trabajo desarrollado por el refuerzo transversal de los elementos R1 se cuantifico con 
las deformaciones medidas al cortante de falla en los estribos de una sola de las cartelas. En el 
elemento TASCa4-R1 se instrumentaron todos los estribos de la longitud acartelada, y en los 
elementos restantes, solamente se monitorearon los estribos ubicados en el vertice y cerca de 
la longitud media de la cartela (figura 15). 

La figura 16 muestra las deformaciones medidas en los estribos de los elementos 
mostrados en la figura 15. A partir de las mediciones, y del hecho que en algunos elementos la 
falla no ocurrio en el lado donde se coloco la instrumentacion, se generalizan algunas 
tendencias observadas en el refuerzo transversal de las TACR: 

1. Los estribos ubicados en el vertice presentaron demandas de deformacion 
superiores a la fluencia. Sin lugar a dudas este comportamiento se asocia a la 
fuerza vertical que introduce el quiebre del acero longitudinal. 

2. Los estribos mas demandados son aquellos que se encuentran ubicados entre la 
longitud media de las cartelas y el vertice. Esto es congruente con la concentracion 
del dano en esta zona a causa de los puntales de compresion. 

3. En promedio se observo un 95% de la fluencia en todos los estribos que fueron 
cruzados por grietas diagonales importantes o en una zona de concentracion de 
agrietamiento. 



V=8 ton 




V=17 ton 




Figura 15 Estribos cruzados por el agrietamiento en el cortante ultimo, a) TASCa4-R1 , y b) TASCa2-R1 
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Deform a cion es en estnbos. Trabe TASCa4-R1 



Deformaciones en estnbos. Trabe TASCa2-R1 



a) b) 

Figura 16 Comportamiento del acero transversal en las TACR: a) a=12.10°, y b) a=6.12 c 



4 CALCULO DE DEFLEXIONES 



Se calcularon las deflexiones al centro del claro asociadas a los cortantes del primer 
agrietamiento diagonal y falla de los elementos, las cuales se compararon contra las obtenidas 
con dos metodos analiticos. Los metodos analiticos fueron: 1) el propuesto por Tena-Colunga 
(1996) que utiliza elementos viga-columna de seccion variable segun la teoria elastica de 
Bernoulli-Euler, y 2) utilizando elementos finitos elasticos lineales tipo placa, los cuales, 
mediante una correcta discretizacion del mallado incluyen directamente la variacion del eje 
neutro y por consiguiente el acoplamiento entre los elementos mecanicos (El-Mezaini et al., 
1991; Balkaya, 2001). Asimismo, se modelo de manera indirecta la condicion agrietada segun 
las limitaciones de cada metodo. La comparacion se limita al calculo de deflexiones, las 
diferencias que pueden obtenerse en el calculo de elementos mecanicos en vigas aisladas y en 
marcos de concreto reforzado, se discuten en Tena (2003). 

4.1 Modelado con elementos viga-columna de seccion variable 



Se utilizo el modelo mostrado en la figura 17. La condicion del primer agrietamiento 
diagonal se modelo utilizando el modulo de elasticidad del concreto medido en los cilindros de 
control (E c ) y la inercia gruesa; asimismo, se utilizo un modulo de Poisson v<=0.20. Ademas, se 
realizaron analisis considerando y sin considerar la inclinacion del eje neutro en los tramos 
acartelados, el cual se define con el angulo <^de la figura 17. 

Para considerar la condicion agrietada, se realizo una reduccion de la rigidez a flexion 
(E c / g ). Esta reduccion se hizo multiplicando la rigidez a flexion por un factor de reduccion de 
agrietamiento, R apr (ecuacion 5). 



V V 




10.00 



10.00 



Figura 17 Modelo utilizado para el analisis de las TACR con elementos viga-columna 
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R 



agr 



agr 



I T 



( 5 ) 



El factor de reduccion por agrietamiento se calculo dividiendo la inercia de la seccion 
transformada-agrietada, l agr , entre la inercia de la seccion transformada de la seccion total, l T , 
en lugar de dividirla entre la inercia de la seccion gruesa, l g , como lo han recomendado algunos 
autores (Paulay y Priestley, 1992) y los reglamentos de mas uso en Mexico (NTCC-04; ACI- 
318-02). Se utilize esta forma de reduccion, ya que, en este caso, se emplearon secciones 
fuertemente armadas por flexion, lo cual obviamente proporciona altas inercias transformadas, 
condicion claramente distinta que al emplear cuanttas ttpicas por flexion, como en las que se 
basan las recomendaciones de Paulay y Priestley y los reglamentos de referenda. Cabe 
mencionar que para secciones ttpicas por flexion, la ecuacion 5 proporciona reducciones 
similares a las recomendadas por los autores y reglamentos mencionados, teniendo como 
ventaja el significado ftsico de dividir secciones con inercias transformadas, y no una seccion 
transformada entre una gruesa. La figura 18 muestra las secciones que representan las 
secciones comentadas en este parrafo. 





SECCION GRUESA SECCION TRANSFORMADA 
(Ig) (It) 



SECCION TRANSFORMADA 
AGRIETADA 



Figura 18 Diferentes idealizaciones para una seccion de concreto reforzado 



El factor R agr se calculo para el mayor peralte total de las cartelas, ya que al ser el 
refuerzo constante en toda la longitud, es en esa seccion donde se obtienen las menores 
cuanttas de acero y por consecuencia las mayores reducciones. Los analisis se realizaron 
reduciendo por el mismo factor toda la longitud de los modelos. 

Los resultados obtenidos demuestran que, al menos para el calculo de deflexiones, no 
tiene una importancia significativa tomar en cuenta la inclinacion del eje neutro de las cartelas, 
por lo que es posible ahorrase el calculo de rotacion de matrices para este fin. 



4.2 Modelado con elementos finitos 



Se construyeron mallas de elementos finitos (figura 19), las cuales fueron 
representativas de cada geometrta. La malla corresponde a la distribucion de la cuadrtcula 
dibujada en los especimenes y que sirvieron para ubicar el agrietamiento progresivo durante los 
experimentos. Para modelar la condicion no agrietada se utilizaron las mismas propiedades del 
concreto empleadas en el modelado con elementos viga-columna de seccion variable. 
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1 @ 3.30 = 3.30 cm 

23 @ 5 = 1 15 cm 1 @ 6.65 = 6.65 cm 




La condicion agrietada de los elementos se realizo disminuyendo significativamente el 
modulo de elasticidad de los elementos de la malla que fueron cruzados por el agrietamiento, 
de esta manera, se obtuvieron modelos con “agrietamiento localizado”. Esta tecnica, tambien 
llamada “agrietamiento de modulos”, ya ha sido empleada para determinar el efecto del 
agrietamiento en elementos estructurales construidos con materiales petreos, tal es el caso de 
la mamposteria (Tena-Colunga y Abrams, 1990; Meli, 1998). 



V = 6 ton. 




Figura 20 Agrietamiento localizado en la trabe TASCa2-R0 (a=6.12°) 



Los porcentajes asignados a los elementos agrietados fueron del 1%, 3% y 5% de E c . 
Vale la pena mencionar que si se introduce un modulo de rigidez nulo en los elementos 
agrietados, se puede crear una condicion de singularidad en la matriz de rigidez, sobre todo 
cuando los patrones de agrietamiento son muy extendidos. La figura 20 muestra el 
agrietamiento al cortante de falla en el elemento TASCa2-R0. Los resultados obtenidos 
mostraron que es suficientemente representative asignar un 5% de rigidez remanente a todos 
los elementos agrietados. 

4.3 Comparacion de metodos 

Las deflexiones calculadas con los dos metodos explicados anteriormente, se 
compararon contra las curvas secantes V-S para los estados del primer agrietamiento diagonal 
y falla obtenidos experimentalmente. Los resultados muestran que ambos metodos tienen la 
misma capacidad de prediccion, la cual tiende a mejorar conforme aumenta el angulo de 
acartelamiento, siendo mas notoria esta tendencia en las trabes con refuerzo transversal 
modeladas con elementos finitos (figura 21). Es evidente que ambos modelados subestimaron 
las deflexiones experimentales, pronosticando en promedio solo el 60% de las deflexiones 
medidas. 
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Curvas secantes en V, y V, en la trabe TASCa2— R0 Curvas secantes en V a9 , y V u en la trabe TASCa3 — R1 

a) b) 

Figura 21 Comportamiento tfpico en las predicciones analfticas de la deflexion al centra del claro en las 
TACR: a) elementos sin refuerzo transversal, y b) con refuerzo transversal 



Aunque el modelado no lineal con ambos metodos presenta varias limitaciones, tales 
como ignorar la presencia del acero longitudinal y, sobre todo, el deslizamiento de las barras y 
apertura de grietas, la mayor discrepancia entre los valores analiticos contra los obtenidos 
experimentalmente se asocia en parte a las incertidumbres en la estimacion de las propiedades 
elasticas iniciales (E c yv), y a una incompatibilidad entre el modo de falla a cortante y la teoria 
elastica para deformaciones pequenas, en la cual se basan los metodos utilizados. El modo de 
falla observado pudiera tener mejor sustento bajo la optica de la Teoria de Puntales y Tensores; 
sin embargo, no debe esperarse necesariamente una mejor prediccion de las deflexiones en las 
TACR con un modelo de armadura que tome en cuenta el agrietamiento de los puntales de 
compresion, segun se recomienda en el Apendice A del reglamento del ACI-318-02, ya que 
tambien tendria limitantes similares a la observada con los elementos viga-columna y 
elementos finitos, sobre todo si se utiliza una teoria basada en deformaciones pequenas. Una 
mejor aproximacion en cuanto a la estimacion de las deflexiones con fallas tan abruptas, como 
las de code, tendria mas sustento con formulaciones analfticas mas complejas, utilizando 
elementos finitos no lineales, donde se consideren modelos con agrietamiento distribuido 
(smeared-crack models), o se discreticen los agrietamientos como discontinuidades con base 
en las teorfas de la mecanica de fracturas. 



5 RESISTENCIA A CORTANTE DE LAS TACR 

La resistencia nominal a cortante (V n ) se determino segun el enfoque tradicional, en el 
cual, la resistencia a cortante de una seccion es la suma de las contribuciones del concreto, V cR 
y del acero de refuerzo transversal, V sR (ecuacion 6). Este enfoque, aunque criticado, aun 
prevalece en la mayorfa de los reglamentos del mundo debido a que representa una 
metodologfa sencilla que proporciona disenos seguros, aunque no explica de manera 
satisfactoria las observaciones experimentales. 

V n =V cR +V sR ( 6 ) 

El calculo de V n se realizo al cortante de falla mediante el diagrama de cuerpo libre 
(DCL), mostrado en la figura 22. Del DCL es obvio que el acero longitudinal inclinado introduce 
una fuerza vertical que no debe ignorarse en el mecanismo resistente a cortante. 
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Figura 22 Diagrama de cuerpo libre a la falla de las TACR 

5.1 Contribucion del concreto 

A partir del DCL de la figura 22 es posible determinar la contribucion del concreto en la 
resistencia a cortante: 

V, R =V„+T X na-'£V l ( 7 ) 

donde V u es el cortante aplicado en la condicion de falla, T es la fuerza de tension en el acero 
longitudinal inclinado, 2t/, es la fuerza de tension de los estribos cruzados por el agrietamiento, 
y a es el angulo de acartelamiento. 

Debido a que la contribucion del concreto se determino con los elementos sin refuerzo 
transversal, el termino 2\/,=0 de la ecuacion 7 es nulo. Por otra parte, el termino V cR , 
propiamente deberfa representar solamente la contribucion a cortante del bloque no agrietado 
en compresion, pero es practica comun suponer que en este termino se incluyan las 
aportaciones de los mecanismos de dovela y trabazon del agregado (Nilson, 1 999). 

La figura 23 muestra el comportamiento de V cR en la historia de carga de las trabes sin 
refuerzo transversal. Se observa claramente que la contribucion del concreto en cada trabe 
tiende a ser mas importante al aumentar el angulo de acartelamiento, en especial una vez 
ocurrido el primer agrietamiento diagonal. Bajo la optica del diseno por secciones, este 
comportamiento se asocia a la fuerza vertical que introduce el acero longitudinal inclinado, la 
cual, para una misma fuerza de tension, es mayor al aumentar el acatartelamiento, lo que obliga 
al concreto a mantener el equilibrio hasta la falla del elemento. 

Se realizo una regresion lineal con los valores calculados de V cR al cortante de falla 

normalizandolos por el factor f cbdmin , contra la tangente del angulo de acartelamiento. Se 

selecciono el peralte efectivo mfnimo, debido a que como se mostro en la tabla 2, con este 
peralte se obtuvieron las mejores estimaciones de la resistencia a cortante en las trabes sin 
refuerzo transversal. La figura 24 muestra esta regresion. 

La contribucion del concreto en las TACR se puede escribir de la siguiente manera: 

VcRTA = 0 ( 8 ) 
donde d crit es el peralte crftico a cortante, el cual se determina con la ecuacion 9: 
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donde h max , h min , a y d min ya se han definido anteriormente, l c es la longitud de la cartela y r el 
recubrimiento del lecho de tension. 






Com porta mien to de V oH en 
la trobe TASCotO-RQ 

a) 



Com portam lento de V cH en 
la trobe TASCa2-RQ 

b) 



Comportamiento de V cH en 
la trobe TASCoi3-RQ 

C) 



Figura 23 Comportamiento de V cfl segun el angulo de acartelamiento: a) a=0°, b) a=6.12°, 

y c) a=9.13° 




Figura 24 Contribucion del concreto en las TACR 

La contribucion del concreto calculada con las ecuaciones 8 y 9, tiene buena 
correspondencia con los involucrados en propuestas anteriores (Debaiky y El-Niema, 1982; 
MacLeod y Houmsi, 1994). La cota superior impuesta al peralte critico tiene la finalidad de 
asegurar que este mantenga su significado fisico, y evitar que numericamente sea mayor que el 
peralte efectivo del apoyo, tal y como podria ocurrir con la propuesta de MacLeod y Houmsi. 

5.2 Contribucion del acero transversal 

La contribucion a cortante del acero de refuerzo transversal, V sR , se determino al 
suponer que los n estribos cruzados por todo el agrietamiento, en la cartela donde ocurrio la 
falla, desarrollaron el esfuerzo de fluencia, es decir: 

v* = nAJ y 



( 10 ) 
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En el calculo de V sR se consideraron todos los estribos cruzados por el agrietamiento, a 
exception de los estribos ubicados a cada lado del estribo de vertice, ya que se supuso que 
estos eran solamente parte del detallado, y controlan exclusivamente la fuerza vertical del acero 
longitudinal inclinado. 



Una manera de comprobar los valores obtenidos de V sR , es compararlo con el termino 
TV, de la ecuacion 7. Lo anterior es posible ya que la contribucion del concreto se puede 
conocer de las ecuaciones 8 y 9, asi como la fuerza vertical que introduce el acero longitudinal 
inclinado se puede estimar de la instrumentation interna. La comparacion entre estos dos 
valores se realizo con el parametro 



7 



V* 



( 11 ) 



La evaluation de rj mostro que es adecuado suponer que los estribos cruzados por 
todo el agrietamiento experimentan un esfuerzo de fluencia, ya que ?/ en promedio tuvo un valor 
de 0.97 con un coeficiente de variation del 24%. 



La evidencia experimental demostro que el estimar la contribucion del acero de refuerzo 
transversal en las TACR como si se tratase de trabes de section constante, es decir, 
suponiendo un agrietamiento de 45° en una proyeccion de un peralte del apoyo, no es 
adecuado, tal y como lo recomiendan Debaiky y El-Niema (1982). Una manera mas realista de 
considerar la contribucion de los estribos en la resistencia a cortante, es cuantificando el 
numero de estribos, n, espaciados cierta distancia s, dentro de la portion de la longitud 
acartelada yt c en la que se extiende todo el agrietamiento (ecuacion 12). El lector interesado 
puede comprobar que en los elementos probados por Debaiky y El-Niema, con geometria 
similar a los comentados en el presente trabajo, se obtuvieron agrietamientos extendidos como 
los mostrados en las figuras 10, 1 2 y 1 5. 

V s rta=Y A ^ L ( 12 ) 

5 




' Yl c 

Figura 25 Extension del agrietamiento diagonal en las cartelas 

Se determino que el agrietamiento se extiende aproximadamente en un 85% de la 
longitud de las cartelas (^=0.85). Considerando esta extension del agrietamiento y una 
profundidad del bloque de compresion del 20% del peralte rmnimo, razonablemente observada 
en los experimentos, se puede obtener un “agrietamiento diagonal equivalente” de 
aproximadamente 24° (figura 25). 



6 METODOLOGIA PROPUESTA PARA EL CALCULO DE LA RESISTENCIA A 
CORTANTE DE LAS TACR 
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Con base en los resultados experimentales, se propone la ecuacion 13 para calcular la 
resistencia nominal a cortante de las TACR: 

Kf/’/l — V C RTA ^sRTA ~ P^sfy SdlO (X (13) 

La contribucion del concreto, V cRT a, se calcula con las ecuaciones 8 y 9, mientras que la 
contribucion del acero de refuerzo transversal, V sRTA , con la ecuacion 12, donde ^0.85. El 
tercer termino de la ecuacion 13 representa la fuerza vertical que introduce el acero longitudinal 
inclinado (T=A s f y ), el factor j3 indica el nivel de fluencia que desarrolla dicho armado, y para el 
cual se recomiendan los siguientes valores: 

[0.1 Trabes sin refuerzo transversal 
P = \ (14) 

[1 .0 Trabes con refuerzo transversal 

6.1 Evaluacion y comparacion de la metodologia propuesta 

Se comparo la resistencia a cortante de las TACR estimada con la ecuacion 13 contra 
la propuesta de Debaiky y El-Niema (1982). Se utilizaron los resultados experimentales de dos 
grupos de elementos: 1) el que corresponde a las trabes con refuerzo transversal del presente 
trabajo, y 2) las trabes que poseen una geometria similar a las de este trabajo, ensayadas por 
Debaiky y El-Niema (1982). La tabla 3 muestra los resultados obtenidos. 

Segun los resultados mostrados en la tabla 3, se observa que de manera general no 
existen diferencias significativas entre ambos metodos para el calculo del cortante resistente. 
Pero un desglose de las fuerzas involucradas en el mecanismo resistente (figura 26), 
demuestra que si existen diferencias entre las dos propuestas, principalmente en la 
contribucion del acero de refuerzo transversal y en la fuerza que introduce el acero longitudinal 
inclinado. En parrafos anteriores se explico como Debaiky y El-Niema consideraron la 
contribucion de los estribos, la cual no tiene una explicacion satisfactoria a partir de la evidencia 
experimental. Respecto a la fuerza que introduce el acero longitudinal inclinado, existe poca 
transparencia en la forma en que la estimaron Debaiky y El-Niema, ya que en un principio 
sostienen que efectivamente se trata de dicha fuerza vertical, pero al final de su publicacion 
mencionan que se trata de la contribucion del mecanismo de dovela del mismo acero. A juzgar 
por estos resultados, es evidente que aun falta mucho por hacer para entender el 
comportamiento a cortante de las TACR. 



Ton 

25.0 

20.0 

15.0 

10.0 
5.0 
0.0 



Trabe TASC a 2-R1 

a=6.12° 




V c rta V s rta pAgfySenoa V nTA V u 

s Este estudio □ Debayky y El-Niema ■ Experimental 



Ton Trabe El (Debaiky y E-Niema (1982)) 

10.0 n a =9.46° 




h Este estudio □ Debaiky y El-Niema ■ Experimental 



a) b) 

Figura 26 Comparacion detallada entre dos metodos que estiman el cortante nominal resistente en 
TACR: a) trabe ensayada por Archundia (2004), y b) trabe ensayada por Debaiky y El-Niema 
(1982). 
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Tabla 3 Comparacion en la prediccion del cortante nominal de las TACR segun dos metodos 

diferentes 


Trabe 


a 


Vu-Exp. 

(ton) 


VnTA-Este 

estudio. 

(ton) 


VnTA-Debaiky y 
El-Niema. 

(ton) 


^ * ril A- Este estudio 


V nT A- Debaiky et al. 


^ u— Exp. 


^ m- Exp. 


Este estudio 


TASCal-RI 


3.07° 


20.00 


21.55 


17.67 


1.08 


0.88 


TASCa2-R1 


6.12° 


17.00 


16.59 


15.95 


0.98 


0.94 


TASCa3-R1 


9.13° 


12.00 


11.23 


13.82 


0.94 


1.15 


TASCa4-R1 


12.10° 


8.00 


5.18 


10.87 


0.65 


1.36 


Promedio 


- 


- 


- 


- 


0.91 


1.08 


a 


- 


- 


- 


- 


0.16 


0.19 


c.v. 


- 


- 


- 


- 


18% 


18% 


Deba 


iky y El-Niema (1982) 


A4 


4.76° 


5.13 


5.14 


5.09 


1.00 


0.99 


C5 


“ 


5.75 


5.43 


5.74 


0.94 


1.00 


E2 


“ 


7.50 


10.18 


7.66 


1.36 


1.02 


A5 


9.46° 


5.70 


1.69 


4.90 


0.30 


0.86 


C4 


“ 


6.10 


2.20 


5.51 


0.36 


0.90 


D5 


“ 


6.50 


8.03 


8.56 


1.24 


1.32 


D6 


“ 


7.50 


8.84 


9.11 


1.18 


1.21 


El 


“ 


9.50 


8.31 


8.18 


0.87 


0.86 


FI 


“ 


6.70 


10.58 


8.86 


1.58 


1.32 


F2 


“ 


7.05 


9.34 


8.54 


1.33 


1.21 


Promedio 


- 


- 


- 


- 


1.02 


1.07 


a 


- 


- 


- 


- 


0.40 


0.17 


c.v. 


- 


- 


- 


- 


39% 


16% 


TOTALES 


Promedio 










0.99 


1.07 


a 










0.35 


0.18 


c.v. 










35% 


17% 



7 CONCLUSIONES 

Se realizo el ensaye de diez trabes acarteladas de concreto reforzado de tamano casi 
real, disenadas para fallar por cortante. Del ensaye de estos elementos, se obtuvieron 
expresiones que permiten determinar el cortante nominal resistente en trabes cuyas cartelas 
tengan aproximadamente una longitud de un tercio del claro efectivo y cuyo angulo de 
acartelamiento oscile entre 3°y 12°. Las ecuaciones propuestas tienen el mismo formato que 
se ha utilizado en el reglamento mexicano para diseno de estructuras de concreto (NTCC-04). 
Las ecuaciones se propusieron para que tuvieran un significado fisico. 

La evidencia experimental mostro un evidente mecanismo de arco, el cual permitio que 
el dano se distribuyera en las cartelas propiciando una falla menos fragil en comparacion con la 
tipica falla subita por cortante que se presenta en trabes prismaticas. A partir de este resultado 
se estudiara el mecanismo a cortante de estos elementos con la Teona de Puntales y 
Tensores. 



La comparacion entre dos metodos analiticos basados en la teoria elastica para el 
calculo de deflexiones en trabes de seccion variable, mostro que al utilizar metodos que 
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modelan las trabes como elementos viga-columna de seccion variable, se obtienen resultados 
comparables con metodos mas refinados como el de los elementos finitos. Para el modelado 
crudo de la condition agrietada, es suficiente con obtener la menor inercia agrietada y aplicarla 
en todo el elemento. Sin embargo, y como se ha demostrado anteriormente para otros 
elementos estructurales construidos con materiales petreos, cuando se utiliza el metodo de 
“agrietamiento de modulos” generalmente se subestiman las deflexiones cuando se comparan 
con resultados experimentales (Tena-Colunga y Abrams, 1990). 

Se presento una formulation que permite estimar el cortante resistente en trabes 
acarteladas de concreto reforzado, la cual, sin duda, debera calibrarse experimentalmente con 
trabes provistas de menores cuantias de acero longitudinal y sometidas a cargas reversibles. 
Ademas se debe estudiar el comportamiento en elementos donde en una misma seccion se 
presenten los maximos valores de cortante y momento, tal es el caso de los voladizos. 
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RESUMEN 

El Indice de Gestion de Riesgos, IGR, que se describe en este documento, trae 
consigo una serie de indicadores que mide el desempeno y la efectividad de la gestion 
de riesgos. Estos indicadores reflejan la organizacion, el desarrollo y la accion 
institucional para reducir la vulnerabilidad y las perdidas, la preparacion para 
responder en caso de crisis y la capacidad para recuperarse con eficiencia. Este 
indice proporciona una medida cuantitativa de la gestion con base en unos niveles 
preestablecidos -targets- o referentes deseables -benchmarks- cualitativos hacia los 
cuales se debe dirigir la gestion del riesgo. El diseno del IGR implica establecer una 
escala de niveles de desempeno o determinar la distancia de las situaciones actuales 
con respecto a ciertos umbrales objetivo o al desempeno logrado por un pais, una 
region subnacional o una ciudad Ifder considerada como referente. Para la formulacion 
del IGR se tuvieron en cuenta cuatro politicas publicas, cada una de e I las descritas 
por seis indicadores: la identif icacion del riesgo, la reduccion del riesgo, el manejo de 
desastres y la gobernabilidad y proteccion financiera. La identificacion del riesgo 
comprende la percepcion individual, la representacion social y la evaluacion objetiva; 
la reduccion del riesgo corresponde a la prevencion y mitigacion, el manejo de 
desastres tiene como objetivo la respuesta y recuperacion post desastre apropiadas; y 
las politicas de gobernabilidad y proteccion financiera estan relacionadas con la 
institucionalizacion y la transferencia del riesgo. Al final se muestran resultados a 
nivel nacional, subnacional y urbano que ilustran la aplicacion del IGR en esas 
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escalas. Palabras clave Gestion de riesgos, desempeno de la gestion de riesgos, 
fndice de gestion del riesgo, toma de decisiones. 



ABSTRACT 



The Risk Management Index, RMI, described in this paper, brings together a group of 
indicators that measure risk management performance and effectiveness. These indicators 
reflect the organizational, development, capacity and institutional actions taken to reduce 
vulnerability and losses in a given area, to prepare for crisis and to recover efficiently from 
disasters. This index is designed to assess risk management performance. It provides a 
quantitative measure of management based on predefined qualitative targets or benchmarks 
that risk management efforts should aim to achieve. The design of the RMI involved establishing 
a scale of achievement levels or determining the distance between current conditions and an 
objective threshold or conditions in a reference country, subnational region, or city. The 
proposed RMI is constructed by quantifying four public policies, each of which is described by 
six indicators. The mentioned policies include the identification of risk, risk reduction, disaster 
management, and governance and financial protection. Risk identification comprises the 
individual perception, social representation and objective assessment; risk reduction involves the 
prevention and mitigation; disaster management comprises response and recovery; and, 
governance and financial protection policy is related to institutionalization and risk transfer. 
Results at the urban, national and subnational levels, which illustrate the application of the RMI 
in those scales, are finally given. 

Keywords: risk management, performance of risk management, risk management index, 
decision making. 



1 ENFOQUE METODOLOGICO UTILIZANDO INDICADORES 

Se han propuesto varios metodos basados en indicadores y otras variables para la 
evaluacion de la vulnerabilidad y del riesgo de desastre. Las contribuciones de Bates (1992) 
(1994), Tucker et al. (1994), Davidson (1997), Puente (1999), Cardona et al. (2003 a, b), UNDP 
(2004), World Bank (2004) y Carreno et al. (2005, 2006), entre otros, han intentado medir 
aspectos relacionados con la vulnerabilidad y el riesgo utilizando indicadores cuantitativos o 
cualitativos. En estos estudios, la vulnerabilidad o el riesgo de desastre ha sido evaluado desde 
diferentes puntos de vista, usando tecnicas que son claramente similares en su metodologfa 
pero con diferentes propositos y alcances al utilizado en el presente artfculo. Los estudios 
mencionados tienen como objetivo la evaluacion de la vulnerabilidad y el riesgo mediante 
indicadores, pero no la evaluacion de la efectividad o el desempeno de la gestion de riesgos, 
que es el objetivo del presente artfculo. 

Actualmente no existen indicadores especfficos en los pafses, ampliamente aceptados, 
para evaluar directamente el desempeno o performance de la gestion del riesgo u otros 
aspectos relevantes que reflejen lo que se desea medir como gestion del riesgo. Se han tenido 
algunas iniciativas para los niveles regional y subnacional (Mitchell, 2003), sin embargo, este 
tipo de iniciativas han sido consideradas subjetivas y arbitrarias, debido a su caracter normativo. 
Uno de los principales esfuerzos para definir los aspectos que caracterizan la gestion del riesgo 
ha sido el marco de accion o framework liderado por la EIRD (2003), en el cual se proponen 
preliminarmente varias areas tematicas, componentes y los posibles criterios para la valoracion 
del desempeno (Cardona et al. 2003b). En cualquiera de los casos es necesario calificar las 
variables con una escala cualitativa que puede ser de 1 a 5 6 de 1 a 7 (Benson 2003b; 
Briguglio 2003 a/b; Michell 2003) o mediante valoraciones lingufsticas Davis 2003; Masure). 
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2 EL I'NDICE DE GESTION Del RIESGO, IGR 



2.1 DESCRIPCION DEL IGR 



El IGR (Disaster Risk Management Index, DRM, en ingles) ha sido disenado para 
evaluar el desempeno de la gestion de riesgos. Proporciona una medida cuantitativa de la 
gestion con base en unos niveles preestablecidos o referentes deseables cualitativos hacia los 
cuales se debe dirigir la gestion del riesgo. El diseno del IGR implica establecer una escala de 
niveles de desempeno (Davis, 2003; Masure, 2003) o determinar la distancia con respecto a 
ciertos umbrales objetivo o al desempeno obtenido por un pais Ifder considerado como el 
referente (Munda, 2003). Para la formulacion del IGR se tuvieron en cuenta cuatro politicas 
publicas, cada una compuesta por seis indicadores que caracterizan el desempeno de la 
gestion en el pais. El fndice de identificacion del riesgo, IGRir es una medida de las 
percepciones individuales de como esas percepciones son entendidas por la sociedad como un 
todo y la estimacion objetiva del riesgo. El fndice de reduccion del riesgo, IGRrr involucra las 
medidas de prevencion y mitigacion. El fndice de manejo de desastres, IGR M d tiene como 
objetivo la respuesta y recuperacion post desastre apropiadas; y el fndice de gobernabilidad y 
proteccion financiera, el IGR PF mide el grado de institucionalizacion y transferencia del riesgo. 
Las cuatro politicas publicas y sus indicadores se definieron despues de un acuerdo con varios 
expertos y evaluadores. Cualquier pais o ciudad podrfa redefinirlos de acuerdo a sus propias 
especificaciones, mientras que los parametros se mantienen en las distintas evaluaciones en el 
tiempo para hacer un seguimiento consistente de la gestion de riesgos. El IGR se define como 
el promedio de los cuatro indicadores compuestos. 

IGR = ( IGR m + IGR rr + 1GR md + IGR PF ) / 4 ( 1 ) 

Se proponen seis indicadores para cada polftica publica que caracterizan el desempeno 
de la gestion de riesgos de un pais, region o ciudad. Un numero muy alto de indicadores podrfa 
ser redundante e innecesario y harfa muy diffcil la asignacion de factores de importancia o 
pesos a cada indicador. Siguiendo el metodo de evaluacion del desempeno de la gestion de 
riesgos propuesto por Carreno et al. (2004), la valoracion de cada indicador se realiza utilizando 
cinco niveles de desempeno (bajo, incipiente, apreciable, notable y optimo) que corresponden a 
un rango de 1 (bajo) a 5 (optimo). Este enfoque metodologico permite utilizar cada nivel de 
referencia simultaneamente como un objetivo de desempeno y, por lo tanto, permite la 
comparacion e identificacion de resultados o logros hacia los cuales los gobiernos deben dirigir 
sus esfuerzos de formulacion, implementacion y evaluacion de cada polftica. Es posible estimar 
alternativamente el IGR como la suma ponderada (pesos) de valores numericos fijos (1 a 5 por 
ejemplo), en vez de los conjuntos difusos de valoracion lingufstica, sin embargo esa 
simplificacion elimina la no linealidad de la gestion del riesgo, obteniendose resultados menos 
apropiados. 

Los subindices de condiciones de gestion de riesgo para cada tipo de polftica publica 
( IR,RR,MD,PF) se obtienen de la ecuacion: 

N 

X VV < 1 l ‘C 

TC'R 1 — 

1 ^ JIX c(IR,RR,MD,PF) N 



i = 1 




(. IR,RR,MD,PF ) 



( 2 ) 
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Donde, w, es el peso asignado a cada indicador, corresponde a cada indicador para 
la unidad territorial en consideracion c y el periodo t -normalizado u obtenido de la 
desfusificacion tie las valoraciones linguisticas; que representan los niveles de desempeno de la 
gestion de riesgo definidos para cada politica publica respectivamente. Dichas valoraciones 
linguisticas, de acuerdo con la propuesta de Cardona (2001) y Carreno (2001, 2006) equivalen 
a un conjunto difuso que tienen una funcion de pertenencia tipo campana y sigmoidal (en los 
extremos), dadas parametricamente por las ecuaciones 

, . , , 1 ( 3 ) 

campana yx,a,b ,c) = — 



a 



sigmoidal (x\a,c) = F n W 

1 + exp [-a(x-c)] 

Donde el parametro b usualmente es positivo y donde a controla la pendiente en el 
punto de cruce, 0.5 de pertenencia, x = c. La figura 1 .a muestra estas funciones de pertenencia. 



Figura 1 a) Funciones que representan el nivel de calificacion, b) Grado de efectividad de la 

gestion. 

La forma y cobertura de estas funciones de pertenencia siguen un comportamiento no 
lineal, en forma de sigmoideo, como lo proponen Carreno et at. (2004) para caracterizar el 
desempeno -o profundidad de la gestion del riesgo y su grado -o probabilidad de su efectividad. 

La respuesta de un sistema socio-tecnico ante el riesgo corresponde a un nivel de 
adaptacion segun el nivel de efectividad de su estructura tecnica y de su organizacion, los 
cuales producen varios patrones de accion, inaccion, innovacion y determinacion ante el riesgo. 
Segun Contort (1999) se pueden presentar varios tipos de respuesta dependiendo de la 
estructura tecnica, la flexibilidad y la apertura cultural al uso de tecnologia. Estos tipos de 
respuesta son: respuesta no adaptativa (inadecuada para nivel de riesgo que existe, el 
desempeno de la gestion es bajo o inexistente), adaptativa emergente (insuficiente, pero 
incipiente), adaptativa operativa (gestion adecuada con restricciones, apreciable) y auto 
adaptativa (innovadora, creativa y espontanea, es decir notable o optima). 

Las funciones de pertenencia para los conjuntos difusos definidas, representan los 
niveles de calificacion para los indicadores y se utilizan durante el procesamiento de la 
informacion. En el eje x de la figura se representa el valor de los indicadores y en el eje y el 
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grado de pertenencia a cada nivel de calificacion, siendo 1 la total pertenencia y 0 la no 
pertenencia. El desempeno de la gestion de riesgos, tal como aquf se propone, lo definen estas 
funciones, que conforman la curva sigmoideo que se ilustra en la figura l.b, donde se indica el 
grado de efectividad de la gestion del riesgo en funcion del nivel de desempeno calculado con 
los diferentes indicadores. Desde el punto de vista teorico, es importante destacar que la figura 
l.b ilustra que el aumento de la efectividad de la gestion de riesgo no es lineal, como en la 
evolucion de cualquier proceso complejo; en un principio se tiene un menor progreso y en la 
medida que se logra una mayor gestion del riesgo, y se hace sostenible, el desempeno 
aumenta y mejora la efectividad. En un alto grado de desempeno, esfuerzos menores 
adicionales aumentan significativamente la efectividad. Por el contrario, pequenos logros en la 
gestion del riesgo se traducen en un desempeno despreciable y poco sostenible, por lo que sus 
resultados tienen poca o ninguna efectividad. 

Es necesario que los expertos, ademas de calificar los indicadores, tambien asignen 
importancias relativas entre los indicadores de cada politica publica, estas importancias o pesos 
son asignados utilizando el Proceso Analftico Jerarquico (PAJ), que se describe en el Apendice 
2.4-3. Un vez han sido ponderadas y agregadas estas funciones forman un conjunto difuso del 
cual se espera obtener una respuesta o resultado. Para lograr esta trasformacion es necesario 
llevar a cabo la desfusificacion de la funcion de pertenencia que se obtiene y extraer de ella su 
valor “concentrado”; lo que equivale a extraer un “indice”. 

Los pesos asignados suman 1 y son utilizados para ponderar (darle altura) las 
funciones de pertenencia de los conjuntos difusos correspondientes a las calificaciones 
realizadas. 



N 

'L w i= l (5) 

j = i 

Donde N es el numero de indicadores que intervienen en cada caso. La calificacion de 
cada politica publica (IGRir, IGRrr, IGR M d y IGR PF ) es el resultado de la union de los conjuntos 
difusos escalados por los pesos, 

Mrmip = max(w 1 xju c {C l ),...,w N xju c (C N )) (6) 

Donde w, a w N son los pesos de los indicadores de la figura 2, mc(Ci) a ju c (C N ) son las 
funciones de pertenencia para las calificaciones de cada indicador, y jurmh es la funcion de 

pertenencia para la calificacion del IGR para cada politica publica p. El valor del indice de 

gestion de riesgos es obtenido de la desfusificacion de esta funcion de pertenencia, utilizando el 
metodo del centroide de area, COA 

IGR p = [max (w, x /j c (C 1 w N x ju c ( C N ))] Centroide (7) 

Esta tecnica consiste en estimar el area y el centroide de cada conjunto difuso y 
obtener el valor concentrado dividiendo la sumatoria del producto entre ellos por la sumatoria 
total de las areas: 



X = 




( 8 ) 



6 
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COA = 



| /j a ( x)xdx 

X 

| //„ (x)dx 



(9) 



Finalmente, el promedio de los cuatro indices proporciona el indice total de gestion de 
riesgos, IGR. 



Figura 2. Indicadores componentes del IGR 



2.2 Indicadores de identif icacion del riesgo 

Es importante reconocer y entender el riesgo colectivo para disenar las medidas de 
prevencion y mitigacion. La identificacion del riesgo significa comprender como se percibe el 
riesgo desde el punto de vista individual y de la sociedad, y cuales son los enfoques 
metodologicos para evaluarlo. Para poder hacer intervenir el riesgo es necesario reconocerlo, 
dimensionarlo (medirlo) y representarlo mediante modelos, mapas, indices, etc. que tengan 
significado para la sociedad y para los tomadores de decisiones. Metodologicamente involucra 
la valoracion de las amenazas factibles, de los diferentes aspectos de la vulnerabilidad de la 
sociedad ante dichas amenazas y de su estimacion como una situacion de posibles 
consecuencias de diferente indole en un tiempo de exposicion definido como referente. La 
valoracion del riesgo con fines de intervencion, es relevante (Carreno et al, 2005; Carreno etal., 
2006) cuando la poblacion lo reconoce y lo entiende. La figura 2 muestra la composicion del 
IGR ir . 
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2.3 Indicadores de reduccion del riesgo 

La principal accion de la gestion de riesgos es reducir el riesgo. En general, 
corresponde la implementacion de medidas estructurales y no estructurales de prevencion y 
mitigacion. Es la accion de anticiparse con el fin de evitar o disminuir el impacto economico, 
social y ambiental de los fenomenos peligrosos potenciales. Implica procesos de planificacion, 
pero fundamentalmente de ejecucion de medidas que modifiquen las condiciones de riesgo 
mediante la intervencion correctiva y prospectiva de los factores de vulnerabilidad existente o 
potencial, y control de las amenazas cuando eso es factible. La figura 2 muestra la composicion 
del IGRrr. 

2.4 Indicadores de manejo de desastres 

El manejo de desastres, IGRmd, corresponde a la apropiada respuesta y recuperacion 
post desastre, que depende del nivel de preparacion de las instituciones operativas y la 
comunidad. Esta polftica publica de la gestion del riesgo tiene como objetivo responder eficaz y 
eficientemente cuando el riesgo ya se ha materializado y no ha sido posible impedir el impacto 
de los fenomenos peligrosos. Su efectividad implica una real organizacion, capacidad y 
planificacion operativa de instituciones y de los diversos actores sociales que se verfan 
involucrados en casos de desastre. La figura 2 muestra la composicion del IGRmd- 



3 EJEMPLOS DE APLICACION 

En el marco del Programa de Indicadores de Gestion de Riesgos de Desastres en las 
America, Colombia y otros 10 pafses de Latinoamerica y el Caribe fueron evaluados (Cardona, 
2005). Adicionalmente, se calculo el IGR para Bogota y para los 32 departamentos de Colombia 
con la participacion representantes de las instituciones encargadas de la ejecucion de las 
polfticas publicas de gestion de riesgos en cada caso. Las valoraciones de los indicadores para 
cada polftica publica y sus respectivos pesos se establecieron mediante consultas con 
evaluadores locales y asesores reconocidos internacionalmente. En esta seccion se presentan 
algunos casos de estudio para ilustrar la aplicacion del IGR a nivel local (urbano), nacional y 
subnacional. 



Nivel Urbano: Bogota, Colombia 



Para la evaluacion de los indicadores de gestion de riesgos se conto con la 
colaboracion de expertos de la Direccion de Prevencion y Atencion de Emergencias de la 
ciudad de Bogota, Colombia (DPAE), y por academicos de la ciudad. En las tablas 2 a 5 se 
pueden ver las calificaciones realizadas en diferentes periodos entre 1985 y 2003 para los 
indicadores mostrados en la figura 2. Como se menciono en la seccion dos, la escala de estas 
calificaciones es: (1) bajo, (2) incipiente, (3) apreciable, (4) notable, y (5) optimo. 



Tabla 1 Calificaciones para los indicadores de identificacion del riesgo (IR) 



Indicador 


1985 


1990 


1995 


2000 


2003 


IR1 


1 


1 


2 


3 


3 


IR2 


1 


1 


2 


3 


3 


IR3 


1 


2 


3 


4 


5 


IR4 


1 


1 


1 


3 


4 


IR5 


1 


1 


2 


2 


3 


IR6 


1 


1 


1 


2 


4 
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Tabla 2 Calificaciones para los indicadores de reduccion del riesgo (RR) 



Indicador 


1985 


1990 


1995 


2000 


2003 


RR1 


1 


2 


2 


3 


4 


RR2 


1 


1 


1 


1 


2 


RR3 


1 


1 


1 


3 


4 


RR4 


1 


2 


2 


3 


4 


RR5 


2 


2 


2 


4 


4 


RR6 


1 


1 


1 


2 


3 



Tabla 3 Calificaciones para los indicadores de manejo de desastres (MD) 



Indicador 


1985 


1990 


1995 


2000 


2003 


MD1 


1 


2 


2 


3 


3 


MD2 


1 


1 


1 


2 


3 


MD3 


1 


1 


1 


2 


2 


MD4 


1 


1 


1 


1 


3 


MD5 


1 


1 


1 


2 


3 


MD6 


1 


1 


1 


1 


2 



Tabla 4 Calificaciones para los indicadores de protection financiera (PF) 



Indicador 


1985 


1990 


1995 


2000 


2003 


PF1 


1 


2 


2 


3 


3 


PF2 


1 


4 


4 


4 


4 


PF3 


1 


1 


3 


3 


4 


PF4 


1 


1 


1 


1 


1 


PF5 


1 


1 


1 


2 


3 


PF6 


1 


1 


2 


2 


3 



Los expertos, ademas de realizar las calificaciones, tambien asignaron importancias 
relativas entre los indicadores de cada politica, a las que se aplico el Proceso Analitico 
Jerarquico (PAJ) para determinar dichos pesos. Este metodo se describe en el Apendice A. Los 
pesos calculados para los indicadores de cada politica publica se indican en la tabla 6. Aunque 
tambien es factible asignar un peso a cada subindice compuesto que el desempeiio de la 
ciudad en cada una de las cuatro politicas, estos pesos se asumieron iguales. 



Tabla 5 Pesos para el conjunto de indicadores 



Peso 


IR 


RR 


MD 


PF 


W] 


0.05 


0.14 


0.11 


0.21 


W2 


0.22 


0.09 


0.11 


0.46 


w 3 


0.36 


0.07 


0.40 


0.12 


w 4 


0.22 


0.31 


0.22 


0.05 


W5 


0.05 


0.20 


0.05 


0.12 


w 6 


0.12 


0.19 


0.11 


0.04 



La figura 3 muestra un ejemplo del calculo de un fndice, en este caso el IGR M d, la figura 
3. a muestra las funciones de pertenencia de las calificacion de los indicadores de manejo de 
desastres para el ano 2003, que se muestra en la tabla 4. La figura 3b muestra la union de los 
conjuntos difusos correspondientes a las calificaciones ponderadas. El proceso de 
desfusificacion correspondiente a la ecuacion 7 y se ilustra en la figura 3b. La tabla 7 muestra los 
resultados finales para la ciudad de Bogota. 
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Figura 2 Agregacion y desfusificacion para calcular el IGR M d para el 2003 



Tabla 6 Indices de gestion de riesgos para Bogota 



Indice 


1985 


1990 


1995 


2000 


2003 


IGR.r 


4.56 


13.90 


35.57 


56.15 


67.10 


IGRrr 


11.03 


13.90 


13.90 


46.14 


56.72 


IGRmd 


4.56 


8.25 


8.25 


24.00 


32.33 


IGRrf 


4.56 


57.49 


54.80 


57.64 


61.44 


IGR 


6.18 


23.38 


28.13 


45.98 


54.40 



Con estos resultados se ilustra como se ha desarrollado la gestion de riesgos en 
Bogota en los ultimos 25 anos, pero tambien enfatiza que aspectos es necesario mejorar de las 
cuatro polfticas publicas que se estudiaron. La politica publica que ha tenido menor desarrollo 
en la ciudad ha sido la politica de manejo de desastres, mientras que la politica que ha tenido el 
mayor desarrollo es la de identificacion del riesgo seguida por la politica de proteccion 
financiera. Se hizo el mismo estudio detallando cada una de las localidades en las que esta 
dividida la ciudad, siguiendo el mismo procedimiento, utilizando las mismas funciones y con 
calificaciones dadas por expertos de la DPAE para el ano 2003. La figura 4 muestra los 
resultados finales del IGR por localidades. Por los resultados, es claro que a nivel urbano, la 
gestion de riesgos es una tarea que fundamentalmente la debe llevar a cabo la administracion 
central de la ciudad. Las localidades no tienen la posibilidad de desarrollarse en esta materia 
independientemente, ya que son areas muy pequehas que no tienen la suficiente capacidad ni 
autonomia para hacerlo. 



Figura 3 Valores del IGR para las localidades de Bogota, 2003. 
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Nivel Nacional: Colombia 

La evaluacion de los indicadores de la gestion de riesgos y sus pesos fueron asignados 
por personal experto de la Direccion de Prevencion y Atencion de Desastres (DPAD) de 
Colombia y por academicos del Centro de Estudios en Desastres y Riesgos (CEDERI) de la 
Universidad de los Andes. Las calificaciones realizadas para diferentes periodos entre 1985 y 
2003, de acuerdo a la escala dada, para los indicadores de las politicas publicas se presentan 
en las tablas 8 a 1 1 . 



Tabla 7 Calificaciones para los indicadores de identificacion del riesgo (IR) 



Indicador 


1985 


1990 


1995 


2000 


2003 


IR1 


2 


3 


3 


4 


4 


IR2 


1 


2 


3 


3 


3 


IR3 


2 


2 


3 


4 


4 


IR4 


1 


1 


2 


3 


3 


IR5 


1 


1 


3 


2 


2 


IR6 


1 


2 


3 


2 


2 



Tabla 8 Calificaciones para los indicadores de reduccion del riesgo (RR) 



Indicador 


1985 


1990 


1995 


2000 


2003 


RR1 


1 


2 


2 


3 


3 


RR2 


1 


2 


3 


2 


2 


RR3 


1 


1 


2 


2 


2 


RR4 


1 


2 


3 


2 


2 


RR5 


2 


2 


3 


4 


4 


RR6 


1 


1 


2 


3 


3 



Tabla 9 Calificaciones para los indicadores de manejo de desastres (MD) 



Indicador 


1985 


1990 


1995 


2000 


2003 


MD1 


1 


2 


2 


3 


3 


MD2 


1 


1 


2 


2 


2 


MD3 


1 


2 


2 


2 


2 


MD4 


1 


1 


1 


2 


2 


MD5 


1 


1 


2 


1 


1 


MD6 


1 


1 


1 


2 


2 



Tabla 10 Calificaciones para los indicadores de proteccion financiera (PF) 



Indicador 


1985 


1990 


1995 


2000 


2003 


PF1 


1 


2 


3 


2 


2 


PF2 


1 


2 


3 


2 


2 


PF3 


1 


1 


2 


2 


2 


PF4 


1 


1 


2 


2 


2 


PF5 


1 


1 


2 


3 


3 


PF6 


1 


2 


2 


3 


3 



Nuevamente, los pesos fueron obtenidos partiendo de la opinion de las autoridades de 
la DPAD mediante la aplicacion del Proceso Analitico Jerarquico. La tabla 12 muestra los 
resultados finales de los indices de gestion de riesgos obtenidos para Colombia despues de la 
union de pesos y la desfusificacion de las calificaciones de cada polftica publica. 
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Tabla 11 Indices de gestion de riesgos para Colombia 



Indice 


1985 


1990 


1995 


2000 


2003 


IGR ir 


10.54 


25.07 


32.46 


48.41 


48.41 


IGRrr 


10.97 


13.96 


39.28 


44.46 


44.46 


IGRmd 


4.56 


12.49 


12.49 


28.73 


28.73 


IGRrf 


4.56 


12.49 


31.50 


39.64 


39.64 


IGR 


7.66 


16.00 


28.93 


40.31 


40.31 



Las figuras 5 y 6 muestran que la identificacion del riesgo y la reduccion del riesgo han 
sido intensivas en Colombia durante los periodos de analisis. De acuerdo a este analisis, 
actualmente, el gobierno de Colombia intenta dirigir sus esfuerzos a formular, implementar y 
evaluar politicas de gestion de riesgos de acuerdo a estos logros y objetivos de desempeno 
(SIRE, 2005). 



Figura 4 Indices de Gestion de Riesgos para cada polftica publica 



Figura 5 Evolucion del IGR desde 1 985 hasta 2003 

Resultados para Latinoamerica y el Caribe 

La gestion de riesgos se evaluo para los siguientes paises de la region de 
Latinoamerica y el Caribe (LAC): Argentina, Chile, Colombia, Costa Rica, Republica 
Dominicana, Ecuador, El Salvador, Guatemala, Jamaica, Mexico y Peru. 

La evaluacion de los indicadores de gestion de riesgos y sus pesos fue realizada por 
asesores nacionales y oficiales de instituciones relacionadas con el manejo de riesgos de 
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desastres de cada pais (Cardona et al. 2004, 2005). Las figuras 7 a 10 ilustran los valores de 
los componentes del IGR y la figura 11 muestra los resultados finales del IGR para los paises 
cada cinco anos desde 1 985 al ano 2000. 



Figura 6 IGR de Identificacion del Riesgo en LAC, 1 985 - 2000 



Figura 7 IGR de Reduccion del Riesgo en LAC, 1 985 - 2000 



Figura 8 IGR de Manejo de Desastres en LAC, 1985 - 2000 





Indice para medir el desempeno de la gestion de riesgos 



37 



Figura 9 IGR de Proteccion Financiera, 1 985 - 2000 



FiguralO IGR paises, 1985 - 2000 



El analisis muestra que Republica Dominicana, Ecuador y Argentina presentan un 
progreso muy bajo en los ultimos anos. El Salvador y Guatemala presentan un desempeno un 
poco mayor que los paises anteriores. Peru y Colombia representan un nivel de desempeno 
mejor que los anteriores, mientras que Chile, Costa Rica, Jamaica y Mexico obtienen los 
avances mas significativos en la practica de la gestion de riesgos. Desde los alios 80 se ha 
presentado una preocupacion creciente con relacion al tema de los desastres, por lo que la 
tendencia en los paises para mejorar la gestion de riesgos ha ido aumentando. Como resultado, 
los avances han pasado de tener un nivel de desempeno “bajo” a un nivel de desempeno 
“significativo” en la mayoria de los casos. En promedio, el desempeno de la gestion de riesgos 
es algo mejor que “incipiente” y la efectividad factible todavia es baja (0.2 - 0.3). Esto sugiere 
que se requieren esfuerzos considerables para promover una gestion de riesgos efectiva y 
sostenible, incluso en los paises mas avanzados, En general, los mayores avances se han 
presentado en la identificacion de riesgos y el manejo de desastres. La reduccion del riesgo, la 
proteccion financiera y la organizacion institucional han presentado un menor avance. 



Nivel Subnacional: Departamentos de Colombia 



La metodologia fue adaptada para evaluar el desempeno de la gestion de riesgos a 
nivel subnacional. El IGR se evaluo para los 32 departamentos de Colombia. La figura 12 
muestra el mapa del IGR. 
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Figura 11 IGR departamentos de Colombia en el 2004 

La evaluacion de los indicadores de la gestion de riesgos se realizo por funcionarios de 
la DPAD, Colombia y por academicos del CEDERI. El IGR solo se evaluo para el ano 2004. El 
Departamento de Antioquia y Bogota D.C. presentaron los avances mas significativos en la 
gestion del riesgo, seguidos por Valle del Cauca, Risaralda, Quindio, Narino, Magdalena, 
Cundinamarca y Caldas que tienen el mismo nivel en gestion de riesgos. Los menores valores 
del IGR los presentan los departamentos de Vichada, Vaupes, Putumayo, Guajira, Guaviare, 
Guainia, Choco, Cordoba, Cesar y Arauca, los cuales presentan un desarrollo minimo de las 
cuatro politicas publicas de la gestion de riesgos. 



CONCLUSIONES 

El Indice de Gestion de Riesgos, IGR, es el primer mdice internacional sistematico y 
consistente desarrollado para rnedir el desempeno de la gestion de riesgos. Las bases 
conceptuales y tecnicas de este mdice son robustas, a pesar de su subjetividad inherente. El IGR 
permite una realizacion sistematica y cuantitativa de cada pals durante diferentes perlodos as! 
como comparaciones entre palses. Este Indice permite realizar la evaluacion de la gestion de 
riesgos a nivel nacional, subnacional y urbano, estableciendo puntos de referenda o patrones de 
funcionamiento de la gestion del riesgo, para definir objetivos de funcionamiento y mejorar su 
eficiencia. 

El IGR es novedoso y con un mayor alcance que otros intentos si mi I arcs en el pasado. 
Puede mostrar rapidamente los cambios que se presentan debido a mejoramientos en las politicas 
o deterioramiento de la gobernabilidad. El IGR puede mostrar mejoramientos anuales o bienales 
puntuales debido a las decisiones politicas y a medidas de implementacion de la gestion de 
riesgos. Es importante desde el punto de vista del refuerzo positivo que se da a los gobiernos 
nacionales o subnacionales, as! como tambien para mejorar la proteccion social y lograr un 
progreso socio-economico. Esta es la razon por la que el IGR es una herramienta util para 
organizaciones como el Banco Interamericano de Desarrollo (BID), que financio este estudio. 
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Este fndice tiene la ventaja de estar compuesto por medidas que forman directamente 
conjuntos de decisiones/acciones especfficas dentro de grupos con resultados deseables. Aunque 
el metodo pueda ser redefinido o simplificado en el futuro, esta aproximacion es muy novedosa 
debido a que permite la medicion de la gestion de riesgos y su efectividad. El nuevo Plan de 
Accion del BID para la reduccion de desastres en Latinoamerica y el Caribe y los artfculos de las 
nuevas estrategias del pals para el desarrollo de asistencia realizada por los tomadores de 
decisiones del banco y de los paises, han estado basados en resultados y analisis de estas 
evaluaciones. 
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APENDICE A 

El PAJ es una tecnica bastante usada para la toma de decisiones con atributos 
multiples (Saaty 1980, 1987; Saaty y Vargs, 1991). Permite la descomposicion de un problema 
en una jerarquia y asegura que tanto los aspectos cualitativos como cuantitativos de un 
problema sean incorporados en el proceso de evaluacion, durante la cual la opinion es extrafda 
sistematicamente por medio de comparaciones entre pares. El PAJ permite la aplicacion de 
datos, experiencia, conocimiento, e intuicion de una forma logica y profunda dentro de una 
jerarquia como un todo. 

El nucleo del PAJ es una comparacion de atributos entre pares ordinales; 
subindicadores en este contexto en los cuales enunciados de preferencia son identificados. 
Para un objetivo dado, las comparaciones son realizadas por pares de subindicadores, primero 
planteando la pregunta “<j,Cual de los dos es el mas importante?” y segundo “^Por cuanto?”. La 
fortaleza de la preferencia es expresada en una escala semantica de 1 a 9, lo que permite la 
medida dentro del mismo orden de magnitud. La preferencia de 1 indica igualdad entre dos 
subindicadores mientras que una preferencia de 9 indica que un subindicador es 9 veces mas 
grande o mas importante que aquel con el que es comparado. Los pesos relativos de los 
subindicadores son calculados usando una tecnica de vectores propios. Una de las ventajas de 
este metodo es que es capaz de verificar la consistencia de la matriz de comparacion a traves 
del calculo de los valores propios (eigenvalues). 

Las matrices permiten la comparacion de la importancia relativa asignada, con sus 
respectivos valores propios, fndice y relacion de consistencia y el calculo de pesos. Los pesos 
calculados fueron utilizados para la evaluacion de los indicadores de cada polftica para el 
ejemplo de Bogota (vertablaA.1 a A. 8). 
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Tabla A. 1 Matriz de comparaciones para la identificacion del riesgo 





IR1 


IR2 


IR3 


IR4 


IR5 


IR6 


IR1 


1 


0.2 


0.2 


0.2 


1 


0.33 


IR2 


5 


1 


0.5 


1 


5 


2 


IR3 


5 


2 


1 


2 


5 


4 


IR4 


5 


1 


0.5 


1 


5 


2 


IR5 


1 


0.2 


0.2 


0.2 


1 


0.33 


IR6 


3 


0.5 


0.25 


0.5 


3 


1 



Valor propio = 6.0877 Cl = 0.018 CR = 0.014 

Tabla A. 2 Importancia para los indicadores de identificacion del riesgo 



Indicador 


Vector propio Principal 


Vector prioritario 


IR1 


0.0982 


0.05 


IR2 


0.4441 


0.22 


IR3 


0.7280 


0.36 


IR4 


0.4441 


0.22 


IR5 


0.0969 


0.05 


IR6 


0.2381 


0.12 



Tabla A. 3 Matriz de comparaciones para la reduccion del riesgo 





RR1 


RR2 


RR3 


RR4 


RR5 


RR6 


RR1 


1 


1 


0.25 


0.5 


3 


1 


RR2 


1 


1 


0.25 


0.50 


3 


1 


RR3 


4 


4 


1 


2 


5 


4 


RR4 


2 


2 


0.5 


1 


5 


2 


RR5 


0.33 


0.33 


0.2 


0.2 


1 


0.33 


RR6 


1 


1 


0.25 


0.5 


3.0 


1 



Valor propio = 6.1343 Cl = 0.027 CR = 0.022 

Tabla A.4 Importancia para los indicadores de reduccion del riesgo 



Indicador 


Vector propio Principal 


Vector prioritario 


RR1 


0.3172 


0.14 


RR2 


0.1896 


0.09 


RR3 


0.1597 


0.07 


RR4 


0.6900 


0.31 


RR5 


0.4382 


0.20 


RR6 


0.4122 


0.19 
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Tabla A.5 Matriz de comparaciones para el manejo de desastres 





MD1 


MD2 


MD3 


MD4 


MD5 


MD6 


MD1 


1 


2 


2 


5 


4 


5 


MD2 


0.5 


1 


1 


5 


2 


5 


MD3 


0.5 


1 


1 


5 


2 


5 


MD4 


0.2 


0.2 


0.2 


1 


0.33 


1 


MD5 


0.25 


0.5 


0.5 


3 


1 


3 


MD6 


0.2 


0.2 


0.2 


1 


0.33 


1 



Valor propio = 6.0684 Cl = 0,014 CR = 0.011 

Tabla A.6 Importancia para los indicadores de Manejo de Desastres 



Indicador 


Vector propio Principal 


Vector prioritario 


MD1 


0.2272 


0.11 


MD2 


0.2272 


0.11 


MD3 


0.8023 


0.40 


MD4 


0.4392 


0.22 


MD5 


0.0923 


0.05 


MD6 


0.2272 


0.11 



Tabla A.7 Matriz de comparaciones para proteccion financiera 





PF1 


PF2 


PF3 


PF4 


PF5 


PF6 


PF1 


1 


0.33 


2 


5 


2 


5 


PF2 


3 


1 


5 


6 


5 


6 


PF3 


0.5 


0.2 


1 


3 


1 


3 


PF4 


0.2 


0.167 


0.33 


1 


0.33 


1 


PF5 


0.5 


0.2 


1 


3 


1 


3 


PF6 


0.2 


0.167 


0.167 


1 


0.33 


1 



Valor propio = 6.0909 Cl = 0.018 CR = 0.,015 

Tabla A.8 Importancia para los indicadores de proteccion financiera 



Indicador 


Vector propio Principal 


Vector prioritario 


PF1 


-0.3942 


0.21 


PF2 


-0.8583 


0.46 


PF3 


-0.2159 


0.12 


PF4 


-0.0887 


0.05 


PF5 


-0.2159 


0.12 


PF6 


-0.0828 


0.04 
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RESUMEN 

El riesgo ha sido definido, para efectos de su gestion, como las consecuencias 
economicas, sociales y ambientales potenciales de eventos peligrosos que pueden ocurrir en un 
periodo de tiempo especffico. No obstante, en el pasado, el concepto de riesgo ha sido definido, 
en muchos casos, de una forma fragmentada, de acuerdo con cada disciplina cientffica 
involucrada en este estudio. Desde la perspectiva de este artfculo, el riesgo requiere de una 
evaluacion multidisciplinar que tenga en cuenta, no solo el dano ffsico esperado y el numero o 
tipo de vfctimas o perdidas economicas, sino tambien las condiciones relacionadas con la 
fragilidad social y la falta de resiliencia que favorecen los efectos de segundo orden (efectos 
indirectos) cuando un evento peligroso impactagolpea un centro urbano. El metodo general 
propuesto de evaluacion del riesgo urbano es multi-amenaza y holfstico: es decir, es un enfoque 
integrado y comprensivo para guiar la toma de decisiones. La evaluacion del dano ffsico 
potencial (enfoque duro) como resultado de la convolucion de la amenaza y la vulnerabilidad 
ffsica de edificios e infraestructura es el primer paso de este metodo. Subsecuentemente, tambien 
se considera un conjunto de condiciones sociales de contexto que agravan los efectos ffsicos 
(enfoque blando). En este metodo, la evaluacion holfstica del riesgo esta basada en la utilizacion 
de indicadores de riesgo urbano. De acuerdo con este procedimiento, el fndice de riesgo ffsico se 
obtiene para cada unidad de analisis de los escenarios de perdida existentes. El fndice de riesgo 
total se obtiene de afectar el primer fndice por un factor de impacto que involucra un coeficiente 

©2006 ESPE, Quito, Ecuador ISSN: 1390-0315 Recibido: Julio de 2006 

Aprobado: Septiembre de2006 



46 



M.L Carreno, O.D Cardona, M.C Marulanda, A.H.Barbat 



de agravamiento el cual esta basado en variables asociadas con las condiciones socio- 
economicas de cada unidad de analisis. Finalmente, el metodo propuesto es aplicado en su forma 
de amenaza individual a la evaluacion holfstica del riesgo sfsmico para las ciudades de Bogota 
(Colombia) y Barcelona (Espana). 

Palabras claves: Enfoque holfstico, evaluacion del riesgo, riesgo sfsmico, vulnerabilidad socio- 
economica. 



ABSTRACT 



Risk has been defined, for management puiposes, as the potential economic, social and 
environmental consequences of hazardous events that may occur in a specified period of time. 
However, in the past, the concept of risk has been defined in a fragmentary way in many cases, 
according to each scientific discipline involved in its appraisal. From the perspective of this 
article, risk requires a multidisciplinary evaluation that takes into account not only the expected 
physical damage, the number and type of casualties or economic losses, but also the conditions 
related to social fragility and lack of resilience conditions, which favour the second order effects 
(indirect effects) when a hazard event strike an urban centre. The proposed general method of 
urban risk evaluation is multi hazard and holistic, that is, an integrated and comprehensive 
approach to guide decision-making. The evaluation of the potential physical damage (hard 
approach) as the result of the convolution of hazard and physical vulnerability of buildings and 
infrastructure is the first step of this method. Subsequently, a set of social context conditions that 
aggravate the physical effects are also considered (soft approach). In the method here proposed, 
the holistic risk evaluation is based on urban risk indicators. According to this procedure, a 
physical risk index is obtained, for each unit of analysis, from existing loss scenarios, whereas 
the total risk index is obtained by factoring the former index by an impact factor or aggravating 
coefficient, based on variables associated with the socio-economic conditions of each unit of 
analysis. Finally, the proposed method is applied in its single hazard form to the holistic seismic 
risk evaluation for the cities of Bogota (Colombia) and Barcelona (Spain). 

Keywords: Holistic approach, risk evaluation, seismic risk, socio-economic vulnerability. 



1 LA NOCION DE RIESGO 

Muchos enfoques conceptuales de riesgo tienen su origen en los estudios de amenaza 
tecnologica y algunos de ellos son extrapolados al campo del riesgo de desastres naturales. Tal 
vez, las primeras investigaciones especializadas en el tema de desastres naturales comenzaron 
a principios de lo sesenta basadas en las contribuciones pioneras de Gilbert White (1964) desde 
el punto de vista de la ecologia y la geografia. Sociologitas como Enrico Quarantelli (1988) y 
Russell Dynes (1994) desde 1963 hicieron esfuerzos para explicar la respuesta social de los 
desastres siguiendo analogfas como la respuesta en caso de ataques nucleares. Geografos como 
Robert Kates (197) y Roger Kasperson (1988) y el ffsico Christopher Hohenemser enfocaron su 
investigacion en riesgos naturales y nucleares. Desde el punto de vista de la ingenierfa civil estas 
investigaciones han sido materializadas en desarrollos realizados en el campo del riesgo ffsico. 
Asf, comenzando del trabajo en evaluacion del dano de Whitman (1973), innumerables 
metodologfas en la evaluacion del riesgo sfsmico ffsico han sido desarrolladas en todo el mundo. 
Mas tarde, este proceso evoluciono hacia una vision mas integral del riesgo sfsmico 
incorporando otros de sus aspectos (Coburnm y Spence 1992) hasta alcanzar la metodologfa 
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generalizada HAZUS (1999), ahora disponible para la evaluacion del riesgo multiamenaza. Durante los 
anos 90, muchos investigadores, motivados por la Decada lnternacional de Reduccion de los Desastres 
Naturales, D1RDN, trabajaron con riesgos y desastres en el mundo. El tema gano importancia y se ha 
venido aumentando el reconocimiento que los terminos amenaza, vulnerabilidad y riesgo han tenido 
diferentes significados e implicaciones desde los diferentes angulos metodologicos y practicos (Cardona, 
2004). 

Un ejemplo de un modelo sistematico de riesgo fue proporcionado por Kates (1971) desde el pensamiento 
de la escuela ecologista. El describe la nocion de “ajuste” de las amenazas naturales 
considerando las interacciones entre naturaleza, hombre y tecnologia. Palmlund (1989) propuso 
un modelo analogico con la estructura clasica de la tragedia Griega (con actores, escenario, 
drama y roles) para explicar el desastre ambiental desde la perspectiva politica y social. Un 
contexto clasico o explicacion estructural, donde el riesgo es visto como un atributo de las 
estructuras sociales, es el propuesto por Douglas y Wildavsky (1982). Una teoria cultural del 
riesgo es propuesta por Rayner (1992), mientras que los enfoques de la escuela de la 
Econorma Politica son propuestos por Wesgate y O’Keefe (1976), Wijkman y Timbarlae (1984), 
Susmal et al. (1984) y Chambers (1989). Las contribuciones de Wisner (1993), Cannon (1994), 
Blaikie et al. (1996) y los miembros de la Red de Estudios Sociales de Prevencion de Desastres 
en Latinoamerica (La RED) (Maskrey, 1994; Lavell, 1996; Cardona, 1996; Mansilla, 1996) 
tambien deben ser consideradas constructivistas, enfatizando la construccion social de la 
vulnerabilidad y el riesgo. Una de las contribuciones conceptuales que se derivan en un enfoque 
multidisciplinario fue hecha por Wilches-Chaux (1989); el propuso diferentes clases de 
vulnerabilidad (cultural, ambiental, social, economica, ffsica, etc.). 

El reporte Analisis de los Desastres Naturales y Vulnerabilidad (UNDRO, 1980), basado 
en la reunion de expertos realizada en 1979, propuso la unificacion de desastre relacionada con 
las definiciones de amenaza (A), vulnerabilidad (V), elementos expuestos (E) y riesgo (R) y 
sugirio una expresion para asociarlas, considerada un estandar en el presente, 

R = EH V (1) 

Muchas metodologfas de evaluacion del riesgo, basadas en esta formulation, han sido 
desarrolladas desde diferentes perspectivas en las ultimas decadas, y, recientemente se desarrollo un 
enfoque holistico o multidisciplinario para el caso de centros urbanos (Cardona y Hurtado 2000; Masure, 
2003). 

Cardona (2001) desarrollo un marco conceptual y un modelo para el analisis de riesgo sismico de una 
ciudad desde una perspectiva hollstica. Este modelo considera las variables de riesgo “dura” y “blando” 
del centra urbano, teniendo en cuenta la exposition, las caracterlsticas socioeconomicas de las diferentes 
localidades (unidades) de la ciudad y su capacidad de resistir desastres o su grado de resiliencia. El 
modelo se construyo para guiar la toma de decisiones en la gestion de riesgos, ayudar a identificar zonas 
crlticas de la ciudad y su vulnerabilidad desde diferentes disciplinas profesionales. 



2 METODOLOGIA DE EVALUACION 

Este artlculo presenta un rnetodo alternative de evaluacion de riesgo urbano basado en el 
modelo de Cardona (Cardona, 2001; Barbat y Cardona, 2003), utilizando un enfoque holistico y 
describiendo el riesgo sismico mediante Indices. El dano esperado en edificios y perdidas en la 
infraestructura, obtenidas de futuros escenarios de perdidas, es information basica para la 
evaluacion del riesgo flsico en cada unidad de analisis. El Indice de dano flsico se obtiene a 
partir de estos datos. 
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El metodo propuesto se desarrolla para una evaluacion multiamenaza y por lo tanto es necesario 
disponer de estimaciones de dano ffsico para todas las amenazas significativas. Frecuentemente, cuando se 
dispone de information historica, la principal amenaza puede ser identificada y de esta forma, su situation 
potencial mas crftica. 

La evaluacion holistica del riesgo mediante indices se logra afectando el riesgo fisico con un factor de 
impacto obtenido de las condiciones contextuales como la fragilidad socioeconomica o la falta de 
resiliencia que agravan el escenario fisico de perdida inicial. Se necesitan datos disponibles de estas 
condiciones a nivel urbano para aplicar el metodo. Mas adelante se hace una explication del modelo y 
tambien se describen algunos ejemplos de aplicacion para las ciudades de Bogota, Colombia, y Barcelona, 
Espana, para ilustrar los beneficios de este enfoque que contribuye a la efectividad del manejo del riesgo, 
invitando a la action, identificando las debilidades duras (hard) y blandas (soft) del centra urbano. La 
fFigura 1 muestra el marco teorico del modelo alternativo. 




Figura 1 Marco teorico y modelo para el enfoque holistico de riesgo de desastres 
(adaptado de Cardona y Barbat, 2000) 



Desde una perspectiva holistica del riesgo, R es una funcion del dano fisico potencial, D ; , y un 
factor de impacto. If. El primero se obtiene de la susceptibilidad de los elementos expuestos Yu„ a 
amenazas, //,, con respecto a sus intensidades potenciales, /, de eventos en un periodo de tiempo f, y el 
ultimo depende de las fragilidades sociales y Fi , y los temas relacionados con la falta de resiliencia, y Ri , de 
que tan propenso esta el sistema socio-tecnico o contexto a un desastre. Usando los meta COnceptOS de la 
teoria de control y sistemas dinamicos complejos para reducir el riesgo, es necesario intervenir de forma 
correctiva y prospectiva los factores de vulnerabilidad y cuando sea posible, intervenir directamente las 
amenazas. El manejo de riesgos requiere de un sistema de control (estructura institucional) y un sistema de 
actuation (politicas publicas y acciones) para implementar los cambios necesarios en los elementos 
expuestos o sistemas complejos donde el riesgo es un proceso social. 
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La evaluacion holfstica del riesgo propuesta en este artfculo se realiza utilizando un conjunto de 
variables de entrada, denominadas descriptores. Estos reflejan el riesgo ffsico y las condiciones de 
agravamiento que contribuyen al impacto potencial. Los descriptores, que sera discutido mas adelante, son 
obtenidos de los escenarios de perdidas y de la informacion de la capacidad socioeconomica y de 
resistencia de un contexto expuesto (Carreno el al. 2005). 

El modelo de evaluacion holfstica del riesgo urbano propuesto en este artfculo mejora los 
aspectos conceptuales y metodologicos de la primera propuesta de Cardona (2001), refinando las tecnicas 
numericas aplicadas y llegando a una herramienta mas versatil. El mejoramiento conceptual proporciona 
un soporte teorico y analftico mas solido al nuevo modelo, eliminando aspectos innecesarios y dudosos del 
metodo anterior y dando mayor transparencia y aplicabilidad en algunos casos. El modelo de Cardona 
permite la evaluacion del riesgo sfsmico en un centra urbano teniendo en cuenta las caracterfsticas del 
riesgo ffsico, amenaza sfsmica, exposicion ffsica, fragilidad socioeconomica y falta de resiliencia, que 
identifican las caracterfsticas de la ciudad que incrementan el nivel de riesgo asf como tambien las areas 
crfticas. Este modelo analiza diferentes tipos de informacion mediante indicadores y utiliza un proceso de 
normalizacion de los resultados basados en la media y la desviacion estandar aplicados a cada indicador. 
Como consecuencia, los resultados obtenidos con el modelo de Cardona permiten una comparacion 
relativa del riesgo sfsmico holistico entre diferentes areas de una ciudad. pero no una comparacion, en 
terminos absolutos, con otras areas urbanas. El modelo de Cardona utiliza un sistema neuro-difuso, con 
conjuntos difusos que identifican las calificaciones lingiifsticas de los descriptores, pero no existe la 
informacion necesaria de la calibracion de este sistema. 

El nuevo modelo aquf propuesto conserva el enfoque basado en indicadores, pero mejora el 
procedimiento de normalizacion y calcula los indices finales de una manera absoluta (no relativa). Esta 
caracterfstica facilita la comparacion del riesgo entre centres urbanos. En este metodo, la exposicion y la 
amenaza sfsmica han sido eliminadas debido a que se considera que estan contempladas en el calculo del 
riesgo ffsico. El descriptor de densidad de poblacion, un componente de exposicion en el modelo de 
Cardona, aquf se considera, como un descriptor de fragilidad social. El nuevo enfoque conserva el uso de 
indicadores y de los conjuntos difusos o funciones de pertenencia originalmente propuestos por Cardona, 
pero de una forma diferente. Otras cambios en el modelo propuesto son las unidades en las que se 
expresan algunos de los descriptores, en ciertos casos es mas importante normalizar los valores de 
entrada con respecto a la poblacion que con respecto al area de la zona de estudio. Por ejemplo, en el caso 
del numero de camas en hospitales existentes en el area urbana estudiada. 

La fragilidad socioeconomica y la falta de resiliencia son un conjunto de factores (relacionados 
con los efectos indirectos o intangibles) que agravan el riesgo ffsico (efectos directos potenciales). Asf, el 
riesgo total depende de los efectos directos, o el riesgo ffsico y los efectos indirectos expresados como un 
factor de efectos directos. Por lo tanto, el riesgo total puede estar expresado de la siguiente forma: 



R t =R f {1 + f) 0 (2) 

Esta expresion se conoce como la ecuacion de Moncho, Rj es el fndice de riesgo total, Rf 6S 6 1 

indice de riesgo ffsico y (1 +F) es el factor de impacto, en el cual F es el coeficiente de 
agravamiento que depende de la fragilidad socio-economica, FS, y de la falta de resiliencia del 
contexto, FR. El riesgo total, para el caso sfsmico tambien se le identifica como Urban Seismic 
Risk Index, USRi, por su denominacion en ingles. 



El riesgo ffsico Rf se obtiene mediante la suma ponderada de factores de riesgo ffsico, 

p 

r f — y>„ x F RFi 

i=i 



( 3 ) 
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Donde F rf , son los factores de riesgo fisico, w RF , es el peso para cada factor y p es el 
numero total de descriptores de riesgo fisico. Los pesos suman 1 y se obtienen utilizando el 
Proceso Analttico Jerarquico (PAJ), el cual es usado para obtener la jerarquta de las 
comparaciones de pares discretos y continuos (Saaty, 2001). Este proceso, explicado en el 
Apendice A, se ha llevado a cabo partiendo de la opinion de expertos, recolectada en las tablas 
A.1 a A. 3, mediante el metodo Delphi. Esta es la forma mas adecuada de calcular la 
importancia relativa de las variables de diferente naturaleza y sus pesos relativos. 

La tabla 1 presenta los descriptores propuestos para el riesgo fisico y las unidades en 
las cuales se utilizan. Es importante mencionar que el calculo de escenarios de riesgo fisico no 
es el objetivo de la metodologfa desarrollada en este artfculo, pero el fndice de riesgo fisico es 
obtenido partiendo de las evaluaciones de perdidas existentes. 



Tabla 1 Descriptores de riesgo fisico, sus unidades e identificadores 



Descriptor 


Unidades 


Xrfi 


Area destruida 


% (area destruida / area construida) 


XrF2 


Muertos 


Numero de muertos cada 1,000 
habitantes 


Xrf3 


Heridos 


Numero de heridos cada 1,000 
habitantes 


Xrf4 


Roturas red de acueducto 


Numero de roturas / Km 2 


Xrfs 


Roturas red de gas 


Numero de roturas / Km 2 


Xrf6 


Longitud de redes electricas cafdas 


m de longitud cafda / Km 2 


Xrf7 


Vulnerabilidad de centrales 

telefonicas 


fndice de vulnerabilidad 


Xrfs 


Vulnerabilidad subestaciones 

electricas 


fndice de vulnerabilidad 


Xrf9 


Dario en la red vial 


fndice de dano 



Los factores de riesgo fisico se obtienen utilizando funciones de transformacion 
mostradas en la figura 2, con las que conocido el valor bruto de cada descriptor (numero 
de fallecidos, area destruida, etc.) se obtiene el valor de cada factor de riesgo fisico, el 
cual toma un valor entre 0 y 1. Por ejemplo, la funcion de transformacion para el area 
construida danada, de acuerdo con la opinion de los expertos, es 0 cuando este 
descriptor tiene el valor rmnimo y es 1 cuando el area danada puede llegar a ser del 20% 
del area construida. 

En forma similar se calcula el coeficiente de agravamiento, F, de la suma 
ponderada de los factores de agravamiento por fragilidad social, F FS /, y por falta de 
resiliencia, F FRy 



( 4 ) 

Donde i/ir FS/ y w FRj son pesos que toman en consideracion su importancia relativa, 
suman 1 y son calculados utilizando el Proceso Analftico Jerarquico (PAJ), m y n indican 
el numero total de descriptores para fragilidad socio-economica y falta de resiliencia 
respectivamente. 



( 5 ) 
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La tabla 2 presenta los descriptores propuestos para conformar el coeficiente de 
agravamiento por fragilidad social y falta de resiliencia. En la tabla 3 se observan los 
descriptores con las unidades utilizadas. Los factores de agravamiento son calculados 
mediante las funciones de transformacion de la figura 3 y la figura 4, que relacionan los 
valores brutos de las variables o descriptores que representan la fragilidad social y la 
resiliencia con los valores correspondientes del agravamiento. El peso w FS o tv FR de cada 
factor F fs/ o F FR; se calcula utilizando el Proceso Analttico Jerarquico (PAJ). 






Figura 2 Funciones de transformacion utilizadas para normalizar los factores de riesgo ffsico. 
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Tabla 2 Descriptores utilizados para estimar el coeficiente F 



Aspecto 


Descriptor 


Fragilidad social 


Area de barrios marginales 


Tasa de mortalidad 


Tasa de delincuencia 


Indice de disparidad social 


Densidad de la poblacion 


Falta de resiliencia 


Camas hospitalarias 


Recursos humanos en salud 


Espacio publico 


Personal de socorro 


Nivel de desarrollo de la localidad 


Operatividad en emergencias 



Tabla 3 Descriptores de agravamiento, sus unidades e identificadores 



Descriptor 


Unidades 


Xfsi 


Area de barrios marginales 


Area barrios marginales / Area localidad 


XfS2 


Tasa de mortalidad 


Numero de muertos cada 10,000 
habitantes 


XfS3 


Tasa de delincuencia 


Numero de delitos cada 100,000 
habitantes 


XfS4 


Indice de disparidad social 


Indice entre 0 y 1 


XfS5 


Densidad de poblacion 


Habitantes / Km 2 de area construida 


Xfri 


Camas hospitalarias 


Numero de camas cada 1 ,000 hab 


XfR2 


Recurso humano en salud 


Recurso humano en salud cada 1 ,000 hab 


XfR3 


Espacio publico 


Area de espacio publico/ Area total 


XfR4 


Personal de socorro 


Personal de socorro cada 10,000 hab 


XfR5 


Nivel de desarrollo de la 
localidad 


Calificacion de 1 a 4 


XfR6 


Operatividad en emergencias 


Calificacion de 0 a 1 * 



• Este indice se definio por Carreno et al. (2005) 



Para determinar las funciones de transformacion, en la mayoria de los casos, se 
utilizaron funciones sigmoidales. Los valores de los descriptores se encuentran en el eje x y los 
factores correspondientes, o los valores escalados, se encuentran en el eje y. Una vez definida 
la forma de estas funciones, se fijaron todos sus valores maximos (correspondientes a los 
valores 1 6 0 para cada factor) utilizando informacion existente de desastres pasados asi como 
tambien apreciaciones de expertos, tanto americanos como europeos. En los casos de falta de 
resiliencia, de los descriptores del nivel de desarrollo de la comunidad y del plan de 
emergencias o preparacion se utilizo una relaciongresion lineal. 
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Area barrios marginales / Area total P[0.05 0.75] 






Figura 3 Funciones de transformacion utilizadas para normalizar los factores de fragilidad social. 

Se considera que los efectos indirectos de los eventos peligrosos, medidos por el factor 
F en la ecuacion 2, pueden ser del mismo orden que los efectos directos. De acuerdo con la 
Comision Economica de Naciones Unidas para Latinoamerica y el Caribe, CEPAL, se estima 
que los efectos economicos indirectos de un desastre natural dependen del tipo de fenomeno 
(Zapata, 2004). El orden de magnitud de los efectos economicos indirectos para un desastre 
‘humedo’ (por ejemplo, causado por una inundacion) podrian llegar a ser del orden de 0.50 a 
0.75 de los efectos directos. En el caso de un desastre ‘seco’ (por ejemplo, un sismo), los 
efectos indirectos podrian llegar a ser del orden de 0.75 a 1 .00 de los efectos directos. La 
diferencia radica en el tipo de dano que ocasionan (destruccion de cultivos, infraestructura, 
viviendas, etc). Esto significa que el riesgo total R T , podrfa ser del orden de 1 .5 y 2 veces el f? F . 
En este metodo se asume que el valor maximo es 2. Por esta razon, en este caso, el factor de 
impacto, F, toma valores entre 0 y 1 en la ecuacion 2. 

La Figura 5 muestra el proceso de calculo del indice de riesgo total para las unidades 
de analisis, que pueden ser distritos, municipios, comunas o localidades. 
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Figura 4 Funciones de transformacion utilizadas para normalizar los factores de falta de resiliencia. 

Expresar el resultado del tndice R F y de Fcomo una combinacion lineal de indicadores 
relativos implica que no existe interaccion entre los mismos o entre dichas variables y los pesos 
que se utilizan en la ponderacion. Esto no es muy realista pero en ocasiones se ha considerado 
aceptable teniendo en cuenta las incertidumbres e imprecisiones inherentes de los datos y para 
efectos de simplificacion. Una vez los descriptores se han escalado y se han convertido en 
factores (conmensurables) se debe realizar su ponderacion segun sea su importancia relativa al 
interior del tndice del cual hacen parte. La participacion de cada factor ilustra que tan importante 
es frente a los demas en la composicion del tndice de riesgo ftsico o en el coeficiente de 
agravamiento. Dado que el riesgo no puede ser medido directamente y los indicadores no estan 
bien correlacionados, la unica opcion que queda para la ponderacion es la evaluacion subjetiva 
por parte de expertos de las importancias relativas de los factores que componen un tndice. 
Aunque esta tecnica parece ser indeseable por algunos por la falta de una base expltcitamente 
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replicable, en evaluaciones caracteristicas dentro de lo que se conoce como ciencia pos-normal 
(Funtowicz y Ravetz, 1992), como es el caso del riesgo de desastre, esta es la unica opcion 
factible y razonable para capturar el criterio, la experiencia y el juicio de expertos. Como 
alternativa se ha propuesto el Proceso Analitico Jerarquico (PAJ) que facilita el analisis 
multicriterio basado en importancias relativas. Una description resumida de esta tecnica se 
presenta en el Apendice A. 



Frfi 


Area destruida 


Wr 


Frf2 


Muertos 


Wr 


Frf3 


Heridos 


Wr 


Frf4 


Roturas red de acueducto 


Wr 


FrF5 


Roturas red de gas 


Wr 


F RF6 


Longitud redes electricas cafdas 


Wr 


F RF7 


Vulnerabilidad centrales telefonicas 


Wr 


Frf8 


Vulnerabilidad subestaciones elect. 


Wr 



Ffsi 


Area de barrios marginales 


Wfs 


F FS2 


Tasa de mortalidad 


Wrs 


FfS3 


Tasa de delincuencia 


Wfs 


F FS4 


Indice de disparidad social 


Wfs 


F FS5 


Densidad poblacional 


Wfs 


Ffr, 


Camas hospitalarias 


Wf 


Ffr2 


Recursos humanos en salud 


Wf 


FfR3 


Espacio publico 


Wf 


Ffr4 


Personal de Socorro 


Wf 


F FR5 


Nivel de desarrollo 


Wf 


F FR6 


Operatividad en emergencias 


Wf 







r t 



Riesgo total 



Figura 5 Factores de riesgo ffsico, fragilidad social y falta de resiliencia y sus pesos 



3 EJEMPLOS DE APLICACION 
3.1 Riesgo sismico de Bogota 

En Bogota, la capital de Colombia, las localidades o alcaldias menores, son 
subdivisiones politico-administrativas del territorio urbano con competencias claras y criterios de 
financiacion y aplicacion de recursos. Se crearon con el objetivo de atender de una forma mas 
eficaz las necesidades de la poblacion de cada territorio. Desde 1992, Bogota esta dividida en 
20 localidades, las cuales se pueden ver en la Figura 6: Usaquen, Chapinero, Santafe, San 
Cristobal, Usme, Tunjuelito, Bosa, Ciudad Kennedy, Fontibon, Engativa, Suba, Barrios Unidos, 
Teusaquillo, Martires, Antonio Narino, Puente Aranda, Candelaria, Rafael Uribe, Ciudad Bolivar 
y Sumapaz. En este estudio, se consideraron solo 19 de estas localidades ya que la localidad 
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de Sumapaz corresponde a area rural. Estas localidades estan subdivididas en 117 unidades 
territoriales (UPZ). 




Figura 6 Division politico-administrativa de Bogota, Colombia. 

Como es bien conocido, la amenaza sismica es la amenaza mas significativa para 
Bogota. El escenario de riesgo sismico fisico ilustrado en la Figura 7 (Universidad de los Andes, 
2005) se utilizo como punto de partida para la aplicacion del modelo. Este mapa representa el 
area destruida en celdas predefinidas considerando que un terremoto con magnitud Ms de 7. ,4 
y un periodo de retorno de 500 anos ocurre en la falla frontal de la cordillera oriental 
(Universidad de Los Andes, 2005). El escenario de riesgo sismico se calculo mediante 
simulaciones edificio por edificio, obteniendo los descriptores del riesgo fisico para cada UPZ. 
Sin embargo, la informacion relacionada con los factores de agravamiento ha sido calculada 
para cada localidad y no para cada UPZ, como se vera mas adelante. 




Figura 7 Un escenario de riesgo sismico fisico (Universidad de los Andes, 2005) 
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Las tablas 4 y 5 muestran los pesos calculados usando el PAJ como se describe en el 
Apendice A, para los componentes del riesgo fisico y para los factores de agravamiento, 
respectivamente. 



Tabla 1 Pesos para los factores de riesgo fisico 



Factor 


Peso 


Valor 


F R FI 


Wrfi 


0.31 


F RFP 


Wrfp 


0.10 


F RFP 


Wrfp 


0.10 


F R FA 


Wrfa 


0.19 


Frff 


Wrfp 


0.11 


Frfr 


Wrfp 


0.11 


F RF7 


WrF7 


0.04 


F RFP 


Wrfp 


0.04 



Tabla 2 Pesos para los factores de las condiciones de agravamiento 



Factor 


Peso 


Valor 


Ffri 


WpRI 


0.18 


F FRP 


Wfrp 


0.04 


F FRP 


Wfrp 


0.04 


F FRA 


WfRA 


0.18 


F FRF 


Wfrp 


0.18 


F FR 1 


Wfri 


0.06 


Ffrp 


Wfrp 


0.06 


Ffrp 


Wfrp 


0.04 


F FRA 


Wfra 


0.03 


Ffrf 


Wfrp 


0.09 


Ffrp 


Webb 


0.09 



Las tablas 6 a 8 muestran los valores de los descriptores utilizados en esta aplicacion, 
los cuales representan el riesgo fisico y la fragilidad social y la falta de resiliencia de la ciudad. 
La tabla 7 muestra los valores de los factores de riesgo fisico obtenidos de la aplicacion de las 
funciones de la figura 4. La tabla 9 muestra los factores de agravamiento de los efectos 
indirectos debidos a la fragilidad social y la falta de resiliencia que se obtienen de la aplicacion 
de las funciones de la figura 2 y 3. El indice de riesgo fisico, fi F , y el coeficiente de 
agravamiento, F, tambien estan indicados en estas tablas. Igualmente, se muestran los valores 
promedio de la ciudad. Estos son calculados normalizando la densidad de poblacion. La tabla 
10 muestra los resultados del riesgo fisico, el factor de impacto, el riesgo total de cada localidad 
y los valores promedios para la ciudad. 

Las figuras 8 a 12 muestran graficamente los resultados de la evaluacion holfstica del 
riesgo sismico de Bogota utilizando el modelo propuesto. Los valores promedio del riesgo fisico 
y el riesgo total por localidad se muestran en las figuras 10 y 12. En estas figuras se puede 
observar que la localidad de La Candelaria presenta la situacion mas crftica desde el punto de 
vista del riesgo sismico fisico y total, ya que su factor de impacto es significativo, aunque no es 
el mayor en la ciudad. Las localidades con mayor factor de impacto son Usme, Ciudad Bolivar, 
Ciudad Kennedy y Bosa, donde los valores mas bajos son los de Barrios Unidos, Teusaquillo y 
Chapinero. El mayor riesgo fisico lo presentan, ademas de La Candelaria, las localidades de 
Santa Fe, Chapinero y Los Martires, mientras que el menor riesgo fisico es el de Ciudad 
Kennedy y Bosa. Los mayores valores de riesgo total aparece en las localidades de La 
Candelaria, Santafe y Los Martires, y los valores menores son los de Ciudad Kennedy, Barrios 
Unidos y Bosa. 
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Indice de riesgo fisico, K, 

□□ o-o.i 
□□ 0 . 1 - 0.2 
[ZZI 0.2 -0.3 
L_! 0.3 -0.4 
_ 0.4 -0.5 

LJ 0.5 -0.6 
0.6 -0.7 
0.7 -0.8 
0.8 -0.9 

19 0.9 -1.0 



Figura 8 Indice de riesgo fisico, ftp, para las UPZ de Bogota 



Tabla 3 Descriptores de riesgo fisico, Rf 



UPZ 


Nombre UPZ 


Local idad 


Xrfi 


XrF2 


XrF3 


Xrf4 


XrF5 


XrF6 


XrF7 


XrF8 


1 


Paseo los 

Libertadores 


San Cristobal 


8.5 


0 


0 


0.96 


0.17 


33.69 


0.68 


0.90 


2 


La Academia 


Suba 


5.6 


0 


1 


0.81 


0.17 


19.14 


0.66 


0.77 


3 


Guaymaral 


Suba 


5.7 


1 


5 


0.93 


0.18 


19.14 


0.66 


0.77 


9 


Verbenal 


Usaquen 


8.4 


1 


3 


0.65 


0.16 


23.87 


0.7 


0.83 


10 


La Uribe 


Usaquen 


12.6 


1 


4 


0.85 


0.16 


23.87 


0.7 


0.83 


11 


San Cristobal 
Norte 


Usaquen 


12.3 


1 


4 


0.67 


0.14 


23.87 


0.7 


0.83 


12 


Toberin 


Usaquen 


13 


6 


19 


1.08 


0.23 


23.87 


0.7 


0.83 


13 


Los Cedros 


Usaquen 


12.8 


8 


17 


0.83 


0.20 


23.87 


0.7 


0.83 
























110 


Ciudad 

Salitre 

Occidental 


Fontibon 


10 


11 


21 


0.89 


0.16 


5.49 


0.64 


0.7 


111 


Puente 

Aranda 


Puente 

Aranda 


21.6 


55 


181 


1.11 


0.27 


20.19 


0.69 


0.7 


112 


Granjas 

Techo 


Fontibon 


30.4 


102 


337 


0.94 


0.22 


5.49 


0.64 


0.7 


113 


Bavaria 


Kennedy 


15.5 


22 


57 


1.14 


0.24 


10.80 


0.54 


0.7 


114 


Modelia 


Fontibon 


11.7 


2 


7 


1.12 


0.29 


5.49 


0.64 


0.7 


115 


Capellania 


Fontibon 


27.2 


23 


74 


1.24 


0.31 


5.49 


0.64 


0.7 


116 


Alamos 


Negativa 


36.3 


21 


63 


1.23 


0.33 


2.82 


0.66 


0.8 


117 


Aeropuerto 
El Dorado 


Fontibon 


16.4 


2 


6 


1.27 


0.34 


5.49 


0.64 


0.7 
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Tabla 4 Factores, F rf , y el indice de riesgo ffsico, Rf 



UPZ 


Nombre 

UPZ 


Localidad 


F RF1 


F RF2 


FrF3 


FrF4 


FrF5 


Frf 6 


FrF7 


FrF8 


Rf 


1 


Paseo los 
Libertadores 


San 

Cristobal 


0.361 


0 


0 


0.019 


0.00243 


0.058 


0.68 


0.9 


0.188 


2 


La 

Academia 


Suba 


0.157 


0 


0.000356 


0.0143 


0.0024 


0.0181 


0.66 


0.77 


0.113 


3 


Guaymaral 


Suba 


0.162 


0.0008 


0.00889 


0.018 


0.00249 


0.0181 


0.66 


0.77 


0.116 


9 


Verbenal 


Usaquen 


0.353 


0.0008 


0.0032 


0.00985 


0.00201 


0.0288 


0.7 


0.83 


0.179 


10 


La Uribe 


Usaquen 


0.726 


0.0008 


0.00569 


0.0157 


0.00209 


0.0288 


0.7 


0.83 


0.298 


11 


San 

Cristobal 

Norte 


Usaquen 


0.704 


0.0008 


0.00569 


0.0105 


0.00158 


0.0288 


0.7 


0.83 


0.290 


12 


Toberin 


Usaquen 


0.755 


0.0288 


0.128 


0.0229 


0.00414 


0.0288 


0.7 


0.83 


0.324 


13 


Los Cedros 


Usaquen 


0.741 


0.0512 


0.103 


0.0148 


0.00332 


0.0288 


0.7 


0.83 


0.318 


























110 


Ciudad 

Salitre 

Occital. 


Fontibon 


0.5 


0.0968 


0.157 


0.0168 


0.00215 


0.00125 


0.64 


0.7 


0.242 


111 


Puente 

Aranda 


Pte.Aranda 


1 


1 


1 


0.0239 


0.00571 


0.02 


0.69 


0.7 


0.584 


112 


Granjas 

Techo 


Fontibon 


1 


1 


1 


0.0182 


0.00371 


0.00125 


0.64 


0.7 


0.579 


113 


Bavaria 


Kennedy 


0.899 


0.387 


0.885 


0.0252 


0.00444 


0.00605 


0.54 


0.7 


0.470 


114 


Modelia 


Fontibon 


0.656 


0.0032 


0.0174 


0.0182 


0.00655 


0.00125 


0.64 


0.7 


0.268 


115 


Capellania 


Fontibon 


1 


0.423 


1 


0.0297 


0.00788 


0.00125 


0.64 


0.7 


0.522 


116 


Alamos 


Negativa 


1 


0.353 


0.949 


0.0292 


0.00886 


0.00045 


0.66 


0.8 


0.515 


117 


Aeropuerto 
El Dorado 


Fontibon 


0.935 


0.0032 


0.0128 


0.031 


0.00902 


0.00125 


0.64 


0.7 


0.358 




Figura 9 indice de riesgo total, Rt, para las UPZ de Bogota 
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Tabla 5 Valores de los descriptores de agravamiento para los factores de fragilidad social 
y falta de resiliencia de Bogota 



Locaiidad 


Xfsi 


XfS2 


XfS3 


XfS4 


XfS5 


Xfri 


XfR2 


XfR3 


Xf=R4 


XfR5 


XfR6 


Antonio 

Narino 


0.015 


398 


120 


0.2 


28338 


4.871 


5 


0.016 


13 


2 


0.0456 


Barrios 

Unidos 


0.002 


917 


130 


0.29 


25920 


1.254 


33 


0.150 


9 


2 


0.0456 


Bosa 


0.249 


1366 


92 


0.51 


44458 


0.338 


3 


0.038 


4 


2 


0.0702 


Chapinero 


0.131 


612 


249 


0 


9255 


16.265 


89 


0.032 


19 


3 


0.0724 


Ciudad 

Bolivar 


0.239 


1866 


78 


0.92 


48968 


0.503 


3 


0.035 


7 


2 


0.1884 


Engativa 


0.061 


2747 


53 


0.41 


34958 


0.183 


7 


0.142 


11 


2 


0.1715 


Fontibon 


0.129 


1028 


34 


0.39 


37558 


0.206 


4 


0.109 


2 


2 


0.0456 


Kennedy 


0.164 


2546 


68 


0.44 


41451 


0.756 


8 


0.084 


1 


1 


0.0456 


La 

Candelaria 


0.248 


100 


86 


0.34 


13074 


2.509 


0 


0.000 


13 


1 


0.0456 


Los Martires 


0.036 


621 


84 


0.33 


32227 


2.846 


103 


0.033 


48 


2 


0.0596 


Puente 

Aranda 


0.001 


1323 


32 


0.37 


37211 


0.616 


4 


0.062 


8 


2 


0.1715 


Rafael Uribe 


0.237 


1618 


109 


0.5 


36759 


1.605 


11 


0.084 


6 


1 


0.1715 


San 

Cristobal 


0.265 


1648 


55 


0.82 


41875 


3.490 


19 


0.108 


11 


1 


0.3383 


Santa Fe 


0.771 


694 


182 


0.36 


17764 


6.176 


143 


0.150 


46 


1 


0.0904 


Suba 


0.110 


2621 


115 


0.41 


25886 


1.181 


15 


0.045 


4 


2 


0.177 


Teusaquillo 


0.002 


681 


90 


0.05 


14437 


5.556 


20 


0.112 


8 


3 


0.1315 


Tunjuelito 


0.161 


941 


63 


0.45 


37702 


2.540 


13 


0.084 


9 


1 


0.0552 


Usaquen 


0.050 


1473 


74 


0.33 


20836 


3.972 


28 


0.129 


3 


2 


0.0456 


Usme 


0.168 


850 


129 


1 


32863 


0.199 


2 


0.021 


12 


2 


0.0456 



Indice de riesgo fisico, R F 



La Candelaria 




— 


Santa Fe 
Chapinero 




























Teusaquillo 














Kalael Unbe 
































San Cristobal 
Engativa 


























Bctiiios Unidos 










Suba 































10,62 



0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 

Figura 10 Indice de riesgo fisico para las localidades de Bogota, en orden descendente 
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Tabla 6 Coeficiente de agravamiento, F, calculado con los factores de agravamiento de 
fragilidad social y falta de resiliencia 



Locatidad 


Fj FS1 


F FS2 


F F S3 


F FS4 


F FS5 


F FR1 


F FR2 


F FR3 


F FR4 


F FR5 


F FR6 


F 


Antonio 

Narino 


0 


0.0155 


0.00538 


0.2 


1 


0.947 


0.778 


0.996 


0 


0.6 


0.95 


0.50 


Barrios 

Unidos 


0 


0.0964 


0.00702 


0.29 


1 


0.997 


0 


0 


0 


0.6 


0.95 


0.44 


Bosa 


0.162 


0.222 


0.00194 


0.51 


1 


1 


0.92 


0.92 


0.367 


0.6 


0.93 


0.61 


Chapinero 


0.0268 


0.0405 


0.0435 


0 


0.125 


0.419 


0 


0.951 


0 


0.3 


0.93 


0.20 


Ciudad 

Bolivar 


0.146 


0.423 


0.00086 


0.92 


1 


0.999 


0.92 


0.936 


0 


0.6 


0.81 


0.67 


Engativa 


0.000494 


0.799 


0.00000988 


0.41 


1 


1 


0.564 


0.00653 


0 


0.6 


0.83 


0.51 


Fontibon 


0.0255 


0.123 


0 


0.39 


1 


1 


0.858 


0.172 


0.837 


0.6 


0.95 


0.55 


Kennedy 


0.053 


0.729 


0.000356 


0.44 


1 


0.99 


0.436 


0.444 


0.959 


1 


0.95 


0.62 


La 

Candelaria 


0.16 


0.00032 


0.00142 


0.34 


0.373 


0.986 


1 


1 


0 


1 


0.95 


0.49 


Los 

Martires 


0 


0.0418 


0.00127 


0.33 


1 


0.982 


0 


0.946 


0 


0.6 


0.94 


0.48 


Puente 

Aranda 


0 


0.208 


0 


0.37 


1 


0.999 


0.858 


0.724 


0 


0.6 


0.83 


0.52 


Rafael 

Uribe 


0.143 


0.315 


0.00382 


0.5 


1 


0.994 


0.142 


0.444 


0.0408 


1 


0.83 


0.56 


San 

Cristobal 


0.189 


0.327 


0.0000274 


0.82 


1 


0.973 


0 


0.18 


0 


1 


0.66 


0.59 


Santa Fe 


1 


0.0532 


0.0191 


0.36 


0.763 


0.915 


0 


0 


0 


1 


0.91 


0.61 


Suba 


0.0147 






















0.52 


Teusaquillo 


0 


0.051 


0.00176 


0.05 


0.494 


0.931 


0 


0.147 


0 


0.3 


0.87 


0.27 


Tunjuelito 


0.0503 


0.102 


0.000185 


0.45 


1 


0.986 


0.0356 


0.444 


0 


1 


0.94 


0.53 


Usaquen 


0 


0.26 


0.000632 


0.33 


0.921 


0.965 


0 


0.045 


0.633 


0.6 


0.95 


0.46 


Usme 


0.0568 


0.082 


0.00685 


1 


1 


1 


0.964 


0.988 


0 


0.6 


0.95 


0.67 



Tabla 7 Riesgo sismico de Bogota 
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Locatidad 


Rf 


F 


Rt 


Antonio Narino 


0,.41 


0,-50 


0,-62 


Barrios Unidos 


0„27 


0,-44 


0,-38 


Bosa 






0,-28 


Chapinero 


0..47 


0,-20 


0,-57 


Ciudad Bolivar 


0,-39 


0,-75 


0,-65 


Engativa 


0,-36 


0,.51 


0,.54 


Fontibon 


0,-41 


0,-55 


0,-64 


Kennedy 


0,.24 


0,-62 


0,.38 


La Candelaria 


0,.62 


0,-49 


0,.93 


Los Martires 


0,.47 


0,-48 


0,-69 


Puente Aranda 


0,.43 


0,.52 


0,-65 


Rafael Uribe 


0,.39 


0,.56 


0,-61 


San Cristobal 


0,.37 


0,.59 


0,.58 


Santa Fe 


0,.48 


0,-61 


0,-77 


Suba 


0,.25 


0,-52 


0,.38 


Teusaquillo 


0,.40 


0,-27 


0,.51 


Tunjuelito 


0,.39 


0,.53 


0,.60 


Usaquen 


0,.30 


0,-46 


0,-44 


Usme 


0,.37 


0,-67 


0,-62 


Bogota 


0.32 


0.55 


0.50 



Coeficiente de agravamiento, F 



Usme _ 
Ciudad Bolivar 
Ciudad Kennedy 
Bosa 
Santa Fe _ 
San Cristobal 
Rafael Uribe 
Fontibon 
Tunjuelito 
Suba 

Puente Aranda 
Engativa 
Antonio Narino. 
La Candelaria 
Los Martires 

Usaquen 
Barrios Unidos 
Teusaquillo 
Chapinero 
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Figura 11 Coeficiente de agravamiento para las localidades de Bogota, en orden descendente 
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Indice de riesgo total, R T 




Figura 12 Indice de riesgo total para las localidades de Bogota, en orden descendente 



Bogota fue estudiada con el primer modelo (modelo Cardona), por lo que es posible 
comparar los resultados obtenidos. La figura 13 muestra los resultados obtenidos (Cardona, 
2001). Este fue el primer analisis integrado del riesgo sismico de la ciudad. Los resultados 
obtenidos con el modelo, solo permiten ordenar las localidades en funcion de su riesgo total 
relativo. Aunque los valores del indice son diferentes, la clasificacion usando ambos modelos, 
es similar. Por ejemplo, la localidad de Tunjuelito presenta el menor valor de riesgo total, y la 
localidad de La Candelaria presenta el mayor riesgo total con el modelo alternativo propuesto y 
serfa el segundo con el modelo anterior. Las otras localidades mantienen un orden similar 
(figuras 12 y 13). 



Riesgo total, I RT 



USME 
CANDELARIA 
SAN CRISTOBAL 
USAQUEN 
BARRIOS UNIDOS 
TEUSAQUILLO 
CHAPINERO 
MARTIRES 
CIUDAD BOLIVAR 
SANTAFE 
RAFAEL URIBE URIBE 
BOSA 
SUBA 

PUENTE ARANDA 
ANTONIO NARINO 
ENGATIVA 
FONTIBON 
CIUDAD KENNEDY 
TUNJUELITO 



Figura 13 Indice de riesgo total para las localidades de Bogota, obtenidas con el modelo de Cardona. 
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3.2 Riesgo Sismico de Barcelona 

La ciudad de Barcelona, Espana, esta subdividida en 10 distritos (ver figura 14), las 
cuales estan dirigidas por un Alcalde. Los distritos tienen competencias en temas de urbanismo, 
espacio publico, mantenimiento de infraestructuras, etc. Estas localidades son: Ciutat Vella, 
Eixample, Sants-Montjui'c, Les Corts, Sarria-Sant Gervasi, Gracia, Horta-Guinardo, Nou Barris, 
Sant Andreu and Sant Marti. Los distritos a su vez se subdividen en 38 barrios o zonas 
estadisticas grandes. Barcelona tambien esta subdividida en 248 zonas estadisticas pequenas 
(ZRP). El indice de riesgo fisico se calculo de un escenario de riesgo probabilista en el marco 
del proyecto RISK-. UE (ICC/CIMNE, 2004). La figura 15 presenta el escenario de riesgo fisico 
calculado, considerando las 248 zonas ZRP. El factor de impacto se calculo por distrito debido a 
la disponibilidad de datos a este nivel solamente. 




Figura 14 Division territorial de Barcelona 



sn 




Grado de dano 



0-0.5 No dano 
0.5 - 0.7 Leve 
0.7- 0.9 
0.9- 1.1 

1.1 - 1.3 

1.3- 1.5 

1.5 - 1.7 Moderado 

1.9- 2.1 

1.7- 1.9 

2.1 - 2.3 

2.3- 2.5 

2.5 - 2.7 Severo 

2.7- 2.9 

2.9- 3.1 

3.1 - 3.3 

3.3- 3.5 

3.5-4 Completo 
Sin datos 



Figura 15 Escenario de riesgo fisico para Barcelona, usando los 248 zonas estadisticas menores (ZRP) 
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La tabla 11 muestra ejemplos de los descriptores de riesgo ffsico para algunos de los 
248 ZRP. La tabla 12 presenta ejemplos de los factores de riesgo ffsico. La tabla 13 muestra los 
valores de los descriptores de fragilidad social y falta de resiliencia, y la tabla 14 muestra los 
factores de agravamiento obtenidos mediante la aplicacion de funciones de transformacion 
(Figuras 2 a 4). Adicionalmente, la tabla 14, abajo, muestra los valores promedio de los factores 
para la ciudad, normalizados usando la densidad de poblacion. La tabla 15 presenta algunos 
ejemplos de los resultados finales del fndice de riesgo ffsico, el factor de impacto y el fndice de 
riesgo total para cada ZRP. Los pesos son los mismos que los utilizados para Bogota (tablas 4 y 
5). Las figura 16 a 18 muestran los resultados del fndice de riesgo ffsico, el factor de impacto y 
el fndice de riesgo total para Barcelona usando el modelo arriba propuesto. 



Tabla 8 Ejemplos de valores de descriptores de riesgo ffsico 




Figura 16 fndice de riesgo ffsico para Barcelona, usando los 248 zonas estadfsticas menores (ZRP) 
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Tabla 9 Factores e fndice de riesgo ffsico, Rf, para Barcelona 



1 ZRP 






F Ft Ft 


F RFd 


F Rf=q 


Frfr 


F RF7 


F HF« 


F RFQ 


mm 




0.952 


0.0288 


0.0512 


0.02 


0.02 


0.02 


0.17 


0.32 


0.025 


0.306 




0.999 


0.08 


0.157 


0.02 


0.02 


0.02 


0.17 


0.32 


0 


0.331 




1 


0.0648 


0.128 


0.02 


0.02 


0.02 


0.17 


0.32 


0 






1 


0.0288 


0.0512 


0.02 


0.02 


0.02 


0.17 


0.32 


0.2 






1 


0.0392 


0.08 


0.02 


0.02 


0.02 


0.17 


0.32 


0.2 


0.340 




1 


[ilikBi 


[IllEift 


0.02 


0.02 


0.02 


0.17 


0.32 


0 






1 


0.0392 


0.0697 


0.02 


0.02 


0.02 


0.17 


0.32 


0.2 






[tW 1 


0.0072 


0.0128 


0.02 


0.02 


0.02 


0.17 


0.32 


0 






0.396 


0.0032 


0.00889 


0.02 


0.02 


0.02 


0.17 


0.32 


0 






0.344 


0.0128 


0.0228 


0.02 


0.02 


0.02 


0.17 


0.32 


0 




























0.0761 


0.0072 


0.0128 


0.02 


0.02 


0.02 


0.17 


0.32 


0.05 


0.053 


241 


0.0128 


0.0128 


0.0228 


0.02 


0.02 


0.02 


0.17 


0.32 


0 


mm 


242 


0.0288 


0.02 


0.0356 


0.02 


0.02 


0.02 


0.17 


0.32 


0.025 


0.041 


243 


0.604 


0.289 


0.564 


0.02 


0.02 


0.02 


0.17 


0.32 


0.025 


0.279 


244 


0.042 


0.0392 


0.08 


0.02 


0.02 


0.02 


0.17 


0.32 


0 




245 


0.353 


0.205 


0.411 


0.02 


0.02 


0.02 


0.17 


0.32 


0 




246 


0.0544 


0.0392 


0.08 


0.02 


0.02 


0.02 


0.17 


0.32 


0.025 


0.054 


247 


0.0544 


0.0512 


0.115 


0.02 


0.02 


0.02 


0.17 


0.32 


0.025 


0.058 


248 


EKB 


p vm 




PXSP9 






[tin 


PIS 


PJI#: 


mm 


Barcelona 


1 0.152 


\0.017 


1 0.033 


\ 0.020 






\0.170 


0.320 


\0.03l\ 


0.076\ 



Tabla 10 Valores para los descriptores de agravamiento de fragilidad social y falta de 

resiliencia de Barcelona. 



Distrito 


Xfri 


•XecJ 


XpFt 


Xffh 


Xfrf 




Xfr? 


Xfri 


Xffm 


Xfrf 


Xfrk 


Ciutat Vella 


PS 


EEEI 


252.87 


0.8 


12690 


4.9650 


11 


0.0828 


15 


1 


1 


EixarriDle 


0 


119 


60.04 


0.3 


14186 


6.1475 


14 


0.0180 


18 


4 


1 


Sant - Montiuic 


0 


Big 


73.61 


0.3 


6834 


0 


0 


0.1219 


15 


3 


1 


Les Corts 


0 


81 


30.99 


0.1 


14080 


10.6864 


24 


0.0424 


18 


4 


1 


Sarria-Sant Gervasi 


0 


95 


30.99 


0 


11647 


10.8704 


24 


0.0194 


8 


4 


1 




0 


115 


42.66 


0.2 


16570 


7.1269 


16 


0.0324 


8 


4 


1 


Horta-Guinardo 


pn 


95 


36.00 


0.5 


21573 


16.1716 


36 


0.0369 


8 


2 


1 


Nou Barris 


0.1 


95 


31.54 


0.8 


28256 


0 


0 


0.0430 


10 


1 


1 


Sant Andreu 


0.1 


91 


31.54 


0.5 


19890 


1.1325 


3 


0.0198 


10 


2 


1 


Sant Marti 


0. 3 


93 


42.44 


0.8 


19069 


0 


0 


0.0337 


3 


1 


1 



Tabla 11 Coeficiente de agravamiento, F, calculado con los factores de agravamiento de fragilidad 
social, F F s, y falta de resiliencia, F F r para Barcelona 



Distrito 


F FF1 


F FF7 


Fffi 


F FRd 


Ffrr 


Ffri 


FpR? 


Ffri 


F FRd 


Ff rf 


F 


Ciutat 


0.0918 


0.00061 


0.0452 


0.8 


0.342 


0.964 


0,142 


0.461 


0 


0.5 


0.4437 


EixamDle 


0 


0.00061 


0.000111 


0.3 


0.471 


0.925 


0,00889 


0.993 


0 


0.5 


0.2796 


Sant 


0 


0.000347 


0.000612 


0.3 


0.0364 


1 


1 


0.0806 


0 


0.5 


0.2558 


Les Corts 


0 


0.000123 


0 


0.1 


0.461 


0.755 


0 


0.893 


0 


0.5 


0.2270 


Sarria- 


0 


0.00026 


0 


0 


0.265 


0.769 


0 


0.991 


0 


0.5 


0.1785 


Gracia 


0 


0.000542 


0 


0.2 


0.678 


0.894 


0 


0.949 


0 


0.5 


0.2947 


Horta- 


0.0102 


0.00026 


0 


0.5 


0.947 


0.436 


0 


0 


0 


0.5 


0.3875 


Nou 


0.0102 


0.00026 


0 


0.8 


1 


1 


1 


0.889 


0 


0.5 


0.6164 


Sant 


0.0102 


0.000215 


0 


0.5 


0.882 


0.997 


0,92 


0.99 


0 


0.5 


0.5042 


Sant 


0.255 


0.000237 


0 


0.8 


0.84 


1 


1 


0.943 


0.633 


0.5 


0.6591 


Barcelona 


0.04 


0.00033 


0.00352 


0.48 


0.69 


0.87 


0.45 


0.75 


0.54 


0.5 


0.42 



Tabla 12 Riesgo sfsmico de Barcelona 
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ZRP 


Rf 


F 


Rt 




001 


0.306 


0.444 


0,442 


002 


0.331 


0.444 


0,479 


003 


0.328 


0.444 


0,473 


004 


0.336 


0.444 


0,485 


005 


0.340 


0.444 


0,491 


006 


0.316 


0.444 


0,456 


007 


0.339 


0.444 


0,489 


008 


0.172 


0.444 


0,248 
















010 


0.126 


0.444 


0,182 












240 


0.053 


0.659 


0,088 


241 


0.032 


0.659 


0,053 


242 


0.041 


0.659 


0,068 


243 


0.279 


0.659 


0,462 


244 


0.048 


0.659 


0,080 


245 


0.183 


0.659 


0,303 


246 


0.054 


0.659 


0,089 


247 


0.058 


0.659 


0,097 


248 


0.078 


0.659 


0,130 


Barcelona 


0.0759 


0.42 


0.1102 



Figura 17 Coeficiente de agravamiento para los distritos de Barcelona 
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Figura 18 Indice de riesgo total, usando las 248 zonas estadfstica menores (ZRP) 

3.3 Comparacion y Discusion de los Resultados 

Los resultados obtenidos para Bogota han sido comparados con los obtenidos para 
Barcelona. La tabla 16 muestra los valores promedios del riesgo para ambas ciudades. La 
ciudad de Bogota esta localizada en una zona con amenaza sismica intermedia, mientras que 
Barcelona esta localizada en una zona con amenaza sismica baja a moderada. Los valores 
promedios obtenidos para el indice de riesgo fisico, R F , reflejan claramente esta situacion. Cabe 
mencionar que los resultados obtenidos del coeficiente de agravamiento, F, en las dos ciudades 
no son muy diferentes. Los valores mas bajos de este factor de impacto son similares (0.20 
para la localidad de Chapinero en Bogota y 0.18 para el distrito de Sarria en Barcelona). La 
diferencia es mas notable para los valores mayores en las dos ciudades (0.67 para la localidad 
de Usme en Bogota y 0.71 para el distrito de Sant Marti en Barcelona). Aunque el valor mas alto 
para Barcelona es mayor que el valor mas alto para Bogota, el valor promedio de Barcelona es 
menor que el de Bogota, lo que refleja la gran diferencia del riesgo sismico holistico entre las 
dos ciudades. La metodologia propuesta, permite una evaluacion holistica del riesgo unificada, 
lo que permite realizar comparaciones entre otras ciudades diferentes del mundo. 

Tabla 1 3 Comparacion de los valores promedio entre Bogota y Barcelona 



Indice 


Bogota 


Barcelona 


Riesgo fisico, R F 


0.23 


0.08 


Factor de Impacto, F 


0.66 


0.42 


Riesgo total, R T 


0.37 


0.11 



4 CONCLUSION ES 
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La estimacion del riesgo requiere un enfoque multidisciplinar que tenga en cuenta no 
solo el dano fisico esperado, el numero y tipo de vfctimas o perdidas economicas sino tambien 
los aspectos sociales, organizacionales e institucionales relacionados con el desarrollo de 
comunidades, que contribuyen con la creacion del riesgo. A nivel urbano, por ejemplo, la 
vulnerabilidad vista como un factor de riesgo inferno, puede ser relacionada no solo con el nivel 
de exposicion o susceptibilidad ffsica de los edificios y la infraestructura potencialmente 
afectados, sino tambien con la fragilidad social y la falta de resiliencia de la comunidad expuesta. 
La ausencia de organizaciones institucionales y comunitarias, una debil preparacion para la 
respuesta a la emergencia, la inestabilidad polftica y una economfa debil en un area geografica, 
contribuyen a aumentar el riesgo. Por todo esto, las consecuencias potenciales negativas no 
solo estan relacionadas con los efectos de los eventos peligrosos, sino tambien con la 
capacidad para absorber los efectos y el control de sus implicaciones en un area geografica 
dada. 



Se ha sugerido una representacion simplificada pero multidisciplinar del riesgo sismico 
urbano para la modelacion, basado en el uso parametrico de las variables que reflejan 
diferentes aspectos del riesgo. Dicho enfoque permite formular un modelo lo mas realista 
posible, en donde mejoras o variables alternativas pueden ser realizadas continuamente. La 
consideracion de aspectos ffsicos permite la construccion de un fndice de riesgo fisico. Tambien, 
la variables contextuales (social, economica, etc.) permitieron la construccion de un coeficiente 
de agravamiento y mediante este un factor de impacto. El primero es construido con la 
informacion de escenarios sfsmicos de dano fisico (efectos directos) y el ultimo es el resultado 
de la estimacion de condiciones de agravamiento (efectos indirectos) basados en descriptores y 
factores relacionados con la fragilidad social y la falta de resiliencia de los elementos expuestos. 

Este nuevo modelo para la evaluacion holfstica del riesgo facilita la gestion integral del 
riesgo por parte de los diferentes actores implicados en la toma de decisiones sobre la 
reduccion del riesgo. Su aplicacion permite hacer un seguimiento de la situacion de riesgo y la 
eficacia de las medidas de prevencion y mitigacion. Los resultados pueden ser analizados y 
permiten establecer prioridades para la mitigacion y acciones de prevencion y planeamiento con 
el fin de modificar las condiciones que tienen una mayor influencia sobre el riesgo de la ciudad. 
Una vez que los resultados han sido expresados en graficas para cada localidad o distrito, es 
facil identificar los aspectos mas relevantes del fndice de riesgo total sin necesidad de 
realizarlos nuevamente para futuros analisis e interpretacion de resultados. Finalmente, este 
metodo permite comparar el riesgo entre diferentes ciudades del mundo y realizar un analisis de 
riesgo multiamenaza. 



APENDICE A CALCULO DE LOS FACTORES DE PONDERACION. 

El Proceso Analftico Jerarquico - (PAJ) es una tecnica usada para la toma de 
decisiones con atributos multiples (Saaty, 1987, 2001; Saaty y Vargas, 1991). Permite la 
descomposicion de un problema en una jerarqufa y asegura que tanto los aspectos cualitativos 
y cuantitativos de un problema sean incorporados en el proceso de evaluacion durante la cual la 
opinion es extrafda sistematicamente por medio de comparaciones entre pares (pairwise). El 
PAJ permite la aplicacion de datos, experiencia, conocimiento e intuicion de una forma logica y 
profunda dentro de una jerarqufa como un todo. 

El PAJ es una metodologfa de decision compensatoria porque las alternativas que son 
eficientes con respecto a uno o mas objetivos pueden compensarse mediante su desempeno 
con respecto a otros objetivos. En particular, el PAJ como metodo de ponderacion permite a los 
tomadores de decisiones derivar pesos contrastantes a los asignados en forma arbitraria (JRC- 
EC, 2002, 2003). 





70 



M.L Carreno, O.D Cardona, M.C Marulanda, A.H.Barbat 



El nucleo del PAJ es una comparacion de atributos entre pares ordinales; 
subindicadores en este contexto, en los cuales los enunciados de preferencia son identificados. 
La fortaleza de la preferencia es expresada en una escala semantica de 1 a 9, lo que permite la 
medida dentro del mismo orden de magnitud. Una preferencia de 1 indica igualdad entre dos 
subindicadores mientras que una preferencia de 9 indica que un subindicador es 9 veces mas 
grande o mas importante que aquel con el que es comparado. Estas comparaciones resultan en 
una matriz (tabla A.1), en la cual, por ejemplo, el factor F RF7 es cinco veces mas importante que 
el factor F RF i. 

Los pesos relativos de los subindicadores son calculados utilizando una tecnica de 
vectores propios. Una de las ventajas de este metodo es que permite verificar la consistencia de 
la matriz de comparacion a traves del calculo de los valores propios y de un tndice consistente. 

El PAJ tolera la inconsistencia a traves de la cantidad de redundancia. Para una matriz 
de tamano n x n, solo comparaciones n-1 son requeridas para establecer pesos para los 
indicadores n. El numero real de comparaciones realizadas en el PAJ es n x n(n-1)/2. Esta 
redundancia es una caractertstica util como lo es, en forma analoga, la estimacion de un 
numero como resultados de obtener el promedio de repetidas observaciones. Esto resulta en un 
conjunto de pesos que son menos sensitivos a errores de juicio. Ademas, esta redundancia 
permite la medicion de los errores de juicio al dar medios para calcular una relacion de 
inconsistencia (Saaty 1980; Karlsson 1998). 

Cl 

CR = (A.1) 

Cl 

random 

Obtenida como la relacion entre un tndice de consistencia 

Cl = ^ max ~ H (A.2) 

n -1 



y el valor del mismo tndice de consistencia Cl random obtenido para una matriz de comparaciones 
generada aleatoriamente, donde X max es el valor propio principal de la matriz de comparacion 
por pares. De acuerdo con Saaty, se asegura una buena exactitud para relaciones de 
consistencia pequenos (CR menor que 0.1 sugiere como regia de medida gruesa, aunque 
incluso 0.2 es citado). Si no se consigue esta condicion, el problema tiene que ser estudiado 
nuevamente y la matriz de comparacion revisada. Una vez que se logra una buena consistencia, 
el vector propio principal es calculado y normalizado. Esta normalizacion se realiza mediante la 
division de cada elemento del valor propio por la suma de los valores de sus elementos. Los 
elementos de este valor propio son los valores de los factores de peso. La tabla A.2 muestra los 
factores de peso obtenidos comenzando por la matriz de comparacion de pares de la tabla A.1 . 
La tabla A. 3 contiene la matriz de comparacion de pares para los factores de agravamiento 
obtenidos de la opinion de expertos y los pesos calculados mediante la aplicacion de PAJ se 
dan en la Tabla A. 4. 



Tabla A 1 Matriz de comparaciones para el riesgo ftsico 





Frfi 


FrF2 


FrF3 


FrF4 


FrF5 


FrF6 


FrF7 


FrF8 
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F RF1 


1 


4 


4 


2 


3 


3 


5 


5 


FrF2 


0.25 


1 


1 


0.5 


1 


1 


3 


3 


FrF3 


0.25 


1 


1 


0.5 


1 


1 


3 


3 


FrF4 


0.50 


2 


2 


1 


2 


2 


4 


4 


FrF5 


0.33 


1 


1 


0.5 


1 


1 


3 


3 


FrF6 


0.33 


1 


1 


0.5 


1 


1 


3 


3 


FrF7 


0.20 


0.33 


0.33 


0.25 


0.33 


0.33 


1 


1 


FrF8 


0.20 


0.33 


0.33 


0.25 


0.33 


0.33 


1 


1 



Valor propio = 8.1 1 
Cl= 0.0152 
CR = 0.01 08 



Tabla A2 Importancia para el riesgo fisico 





Vector 

propio 

principal 


Vector 

prioritario 


Frfi 


0.7410 


0.31 


Frf? 


0.2420 


0.10 


Frfi 


0.2420 


0.10 


FrF4 


0.4368 


0.19 


Frfr 


0.2496 


0.11 


Frfr 


0.2496 


0.11 


FrF7 


0.0958 


0.04 


Frfr 


0.0958 


0.04 



Tabla A3 Matriz de comparaciones para el coeficiente de agravamiento 





Ffsi 


FfS2 


FfS3 


FfS4 


FfS5 


Ffri 


FfR2 


FfR3 


FfR4 


FfR5 


Ffrb 


Ffsi 


1 


4 


4 


1 


1 


3 


3 


4 


4 


3 


3 


FfS2 


0.25 


1 


1 


0.25 


0.25 


0.5 


0.5 


1 


1 


0.5 


0.5 


FfS3 


0.25 


1 


1 


0.25 


0.25 


0.5 


0.5 


1 


1 


0.5 


0.5 


FfS4 


1 


4 


4 


1 


1 


3 


3 


4 


4 


3 


3 


FfS5 


1 


4 


4 


1 


1 


3 


3 


4 


4 


3 


3 


Ffri 


0.33 


2 


2 


0.33 


0.33 


1 


1 


2 


2 


0.5 


0.5 


FfR2 


0.33 


2 


2 


0.33 


0.33 


1 


1 


2 


2 


0.5 


0.5 


FfR3 


0.25 


1 


1 


0.25 


0.25 


0.5 


0.5 


1 


2 


0.33 


0.33 


FfR4 


0.25 


1 


1 


0.25 


0.25 


0.5 


0.5 


0.5 


1 


0.33 


0.33 


FfR5 


0.33 


2 


2 


0.33 


0.33 


2 


2 


3 


3 


1 


1 


FfR6 


0.33 


2 


2 


0.33 


0.33 


2 


2 


3 


3 


1 


1 



Valor propio = 1 1 .24 
C/ = 0.024 
CR = 0.01 6 
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Tabla A4 Importancia para el coeficiente de agravamiento 





Vector 

propio 

principal 


Vector 

prioritario 


Ffsi 


1.0000 


0.18 


F FS2 


0.2136 


0.04 


FfS3 


0.2136 


0.04 


FfS4 


1 


0.18 


FfS5 


1 


0.18 


Ffri 


0.33928 


0.06 


F FR2 


0.33928 


0.06 


FfR3 


0.21601 


0.04 


F FR4 


0.1895 


0.04 


FfR5 


0.47833 


0.09 


FfR6 


0.47833 


0.09 
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RESUMEN 

En los procedimientos de analisis y diseno sismorresistentes de estructuras se hace 
necesario introducir simplificaciones o aproximaciones que permitan estudiar la 
equivalencia entre el analisis elastico a partir del cual se efectuan los calculos, con el 
comportamiento elasto plastico, que es el comportamiento real de las estructuras cuando 
son sometidas a las la accion de terremotos fuertes. Una de estas simplificaciones, consiste 
en proyectar las estructuras de manera que estas dispongan de una resistencia lateral 
inferior a la resistencia elastica, lo cual se consigue mediante la aplicacion de cargas 
laterales calculadas a partir de espectros de diseno inelasticos, que no son mas que los 
espectros de diseno elasticos afectados por factores de reduccion de respuesta tambien 
conocidos como factores R. En este artfculo se realiza una revision exhaustiva del estado 
del arte de los factores de reduccion de respuesta, que comprende desde su propuesta 
original, pasando por su implementacion en los codigos de diseno sismorresistente, hasta 
llegar a las nuevas propuestas de racionalizacion de su aplicacion. Para poder alcanzar una 
mejor comprension sobre el tema, es necesario realizar una comparativa entre la aplicacion 
de los factores de reduccion de respuesta en los diferentes codigos de diseno 
sismorresistente a nivel mundial. 

ABSTRACT 

In the procedures of seismic analysis and design of structures, it is necessary to 
introduce simplifications that allow achieving equivalence between results of elastic analysis, 
on which the seismic codes are based, with results considering the nonlinear, plastic, 
behavior of structures, which correspond to their real behavior when subjected to strong 
earthquake actions. One of these simplifications consists of designing the structures in such 
a way to assure lower lateral strength than the elastic one. This lateral strength is obtained 
by applying calculated equivalent seimic loads using inelastic design spectra, obtained from 
the elastic ones by dividing them by the response reduction factors, also known as R factors. 
In this study, an exhaustive revision of the state-of-the-art of these response reduction 
factors is made, which starts from its original proposal, goes through its implementation into 
the seismic design codes, until arriving to the new proposals of application. Finally, a 
comparative study among the applications of the response reduction factors in different 
seismic design codes existing in the world is made. 
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1 FACTORES DE REDUCCION DE RESPUESTA 



El concepto de factor de reduccion de respuesta se basa en la premisa de que un 
sistema estructural bien detallado es capaz de sostener grandes deformaciones sin llegar a 
colapsar. Al aplicar en el proyecto de estructuras sismorresistentes factores de reduccion 
mayores que la unidad, el proyectista acepta una simplificacion importante: la de que con 
las herramientas de calculo lineal se pueden obtener unas cuantificaciones razonables de la 
respuesta real de las estructura. La segunda simplificacion asumida es que si se acepta un 
comportamiento significativamente no lineal, es logico esperar que ocurra un dafio global 
importante en la estructura. 



2 ANTECEDENTES 



La adopcion de los factores de reduccion para determinar fuerzas sismicas de 
diseno a partir de las correspondientes a la respuesta elastica se ha realizado partiendo de 
los valores sugeridos inicialmente por Veletsos y Newmark (1960), que vincularon los 
factores de reduccion de respuesta con los valores esperados de ductilidad global de 
desplazamiento, de acuerdo con los rangos de periodos dentro del espectro de 
aceleraciones. Los valores adoptados obedecen a los criterios derivados de la comparacion 
de los espectros elasticos e inelastico de respuesta, los cuales mostraban iguales fuerzas 
para el rango de frecuencias altas, igual energia absorbida para el rango de frecuencias 
intermedias e iguales desplazamientos para las frecuencias bajas. Esto permitio plantear los 
siguientes valores para los factores de reduccion de la Tabla 1 . 



Tabla 1 Valores del factor de reduccion de respuesta propuestos por Newmark y 

Veletsos, 1960. 



R = 1 


Para T = 0 


R = p 


Para el rango de los periodos cortos. 


R = 


Para el rango de los periodos moderados. 


R = p 


Para el rango de los periodos intermedios y largos. 



En esta tabla, R es el factor de reduccion, p es la ductilidad global y T el periodo de la 
estructura. El procedimiento se ha mantenido a lo largo de los anos en los codigos de 
diseno sismorresistente debido al buen desempeno exhibido por las estructuras disehadas 
con este metodo, que han mantenido su integridad y preservado vidas al haber sido 
sometidas a sismos fuertes. La dependencia del factor de reduccion respecto a la ductilidad 
global y del periodo estructural ha prevalecido en los codigos de diseno, a pesar de que 
recientes investigaciones han demostrado que aplicar dichos factores redunda en diseno 
poco seguro para el rango de los periodos bajos y excesivamente conservador para los 
periodos intermedios y largos, (Ordaz y Perez- Rocha 1998). 






Factores de reduccion de respuesta estado del arte y estudio comparative entre codigos 79 



Al nivel normativo, los factores de reduccion son considerados por primera vez en el 
articulado del Blue Book (SEAOC, 1959) de forma explicita en el calculo del esfuerzo 
cortante minimo considerado en el analisis de edificios, segun la expresion: 

V b = C • K.W (1) 



En esta expresion V b es el cortante en la base rmnimo, C es un factor a 

determinar que depende del periodo fundamental de la estructura, K es un factor de 
reduccion que depende del tipo estructural considerado y cuyo valor es el inverso del factor 
de reduccion de respuesta propiamente dicho y W es el peso sismico total del edificio. 
Desde esta primera aproximacion se realiza la recomendacion de que los factores de 
reduccion deberian ser aplicados de acuerdo con el juicio del proyectista, pero requiriendo 
profundizar en la determinacion mas precisa de los valores de estos. 

Posteriormente, el UBC (ICBO, 1961) incluye el factor de zonificacion Z en la 
expresion de calculo del cortante rmnimo de diseno, manteniendo el factor de reduccion en 
el formato de K de la expresion anterior. 

V b = Z • C • K.W (2) 



Los factores de reduccion R , tal como se les conoce actualmente son presentados 
en la publicacion del ATC-3-06 (ATC 1978), que modernizo todo el procedimiento general 
aplicado en la ingenieria sismica. El procedimiento propuesto, que permite la aplicacion del 
metodo de las fuerzas sismicas equivalentes mediante la reduccion del cortante en la base 
calculado a partir del espectro de diseno elastico con un amortiguamiento igual al 5% del 
amortiguamiento critico y para el correspondiente peso sismico de la edificacion: 



V S fl5 xW 

y e <0 



R 



R 



( 3 ) 



El cortante en la base de diseno, V b , se calcula escalando el cortante en la base 
elastico, V e , por el factor de reduccion R correspondiente al diseno. En el calculo del 
cortante en la base elastico V e se considera la aceleracion del espectro elastico de diseno 
calculado por los metodos dinamicos conocidos, para un amortiguamiento del 5% del 
amortiguamiento critico, lo que permite obtener la aceleracion de diseno S a 5 . Esta 

aceleracion afecta al peso sismico W que, generalmente, corresponde al 100% de las 
cargas permanentes mas un porcentaje de las sobrecargas de uso. 

Posteriormente, la adecuacion del las normativas de diseno sismorresistente al analisis 
de secciones basado en resistencia que sustituia al analisis por esfuerzos admisibles, 
propicio la aparicion del factor de reduccion R w , que en terminos generates tenia valores 
mayores que su predecesor. Este nuevo factor de reduccion se encuentra en el articulado 
del UBC-97 y del IBC-2003, sin embargo en las ultimas revisiones de dichos codigos han 
sufrido importantes reducciones respecto a las versiones anteriores. 



3 COMPONENTES DEL FACTOR DE REDUCCION DE RESPUESTA R 

Seguidamente se presentan tres de las propuestas que han procurado establecer 
los componentes del factor de reduccion de respuesta. 
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3.1 Propuesta de Bertero (1986) 

Esta propuesta ha sido posteriormente apoyada por Witthaker et al. (1999), e 
incluye la redundancia estructural como un factor independiente en la determinacion de los 
factores de reduccion de respuesta. Sin embargo es necesario senalar que hasta la fecha 
se carece de investigaciones que confirmen su grado de participacion o rango de valores 
aplicables a las distintas tipologias estructurales. La propuesta se muestra en la expresion 
siguiente: 



R - Rs R h R vg r sr 



(4) 



De acuerdo con la expresion anterior, el factor de reduccion R queda definido como el 
producto de cuatro factores: el factor de reduccion por resistencia R s , el factor de reduccion 
por ductilidad R^ , el factor de reduccion por multiples grados de libertad R VG y el factor de 
reduccion por redundancia R sr . 

3.2 Propuesta de Witthaker et al. (1987) 

En esta propuesta se estima que el factor de reduccion de respuesta es el producto 
de tres factores interrelacionados entre si: 



R “ R S -R n -R i; 



(5) 



Los factores son: el factor de reduccion por resistencia R s que contiene la reduccion debida 
a las consideraciones conservadoras asumidas en el diseno en general, el factor de 
reduccion por ductilidad R^ que presenta la influencia de la ductilidad y el factor de 

reduccion por amortiguamiento R«= , que incluye el efecto del amortiguamiento propio de los 

materiales que constituyen la estructura. De estos tres factores, los dos primeros pueden 
determinarse directamente de la curva de capacidad, en base a la relacion de las fuerzas 
cortantes en la base, asi el factor de resistencia R s es: 



R s = 



Vk 



( 6 ) 



Siendo V u el cortante en la base ultimo y V b el cortante en la base de diseno. 

El factor de reduccion por ductilidad R^ , viene dado por la relacion entre los cortantes en la 
base elastico V e y ultimo V u de la curva de capacidad de la estructura. 



R =Xi 
^ V 



(7) 



Witthaker et al. (1987) proponen el calculo de los factores de reduccion por amortiguamiento 
R;, a partir de las caracteristicas dinamicas de la estructura, representadas por el periodo 

fundamental, y de acuerdo con el amortiguamiento de los materiales que constituyen la 
estructura. 



3.3 Propuesta ATC-19 (1995) 

En la propuesta del ATC-19 (1995), el factor de amortiguamiento ha sido excluido. Se 
considera que el factor de reduccion de respuesta R es el producto de tres factores: el de 
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reduccion por resistencia R s , el de reduccion por ductilidad R^ y otro factor, denominado 
factor de redundancia R r , segun la expresion siguiente: 

R = (R s .R^).R r (8) 

De estos tres factores los dos primeros son dependientes del periodo, y se definen 
de forma similar a la expuesta en el apartado anterior, mientras que los valores del factor de 
redundancia suelen presentarse como valores fijos sugeridos para distintos tipos 
estructurales, sin que exista claridad en cuanto a su determinacion ya que suele incluirse 
dentro del factor de sobre resistencia. 



4 RECIENTES INVESTIGACIONES SOBRE LOS COMPONENTES DE R 

Seguidamente, se presenta una resena de recientes investigaciones efectuadas 
sobre cada uno de los componentes del factor de reduccion de respuesta. 



4.1 Factores de reduccion basados en ductilidad 

Dada la evidente importancia de los factores de reduccion basados en ductilidad 
R^ en el calculo de factores de reduccion de respuesta, se incluye en este apartado una 

revision del estado del arte que contempla trabajos realizados aplicando diferentes 
enfoques, asi como modelos estructurales. 

Nassar y Krawinkler (1992) 

Utilizaron el modelo elasto plastico bilineal en el calculo de la respuesta inelastica de 
sistemas de un grado de libertad, para diferentes niveles de ductilidad global. Los espectros 
de R^ se expresan en funcion del periodo estructural T , de la ductilidad global p y del 

endurecimiento post-cedente a : 



=[c(pi-l) + l]/^ 



c(T,a) 



T a _b 
1 + T a + T 



( 9 ) 



Donde los parametros de ajuste coeficiente a y bdependen de la rigidez post-cedente a . 

Vidic et al. (1994). 

Plantearon espectros de R tl obtenidos como la idealizacion de espectros obtenidos a partir 

de modelo histeretico poligonal, cuyos parametros de ajuste se calculan a partir de las 
caracteristicas del material de la estructura. El espectro presenta dos ramas lineales de 
pendientes del periodo: 
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T < T 0 

R,=c>-i) c «-^+i 

T 0 (10) 

T>T 0 

Rn = Cj (|x i) Cr + 1 

Donde el periodo caracteristico de la estructura T 0 depende de la ductilidad global, la 
relacion entre la demanda-capacidad de la estructura y el tipo de material. Los parametros 
C, y C R dependen del comportamiento histeretico y del amortiguamiento respectivamente. 

Miranda y Bertero (1994) 

Esta investigacion introdujo el efecto de sitio en el calculo de factores de reduccion por 
ductilidad, para lo que utilizo el analisis de la respuesta no lineal de los suelos para 
determinar la respuesta en superficie y el modelo elasto plastico perfecto para determinar la 
respuesta inelastica de estructuras modeladas como sistemas de un grado de libertad. La 
expresion utilizada es: 

R m =-^7“ + 1 (11) 

<P 

El factor ()) contiene la dependencia del periodo y la informacion relativa a la influencia del 
tipo de suelo: 

<(> = 1 + — — • e L5(T 0 6) Roca (suelo tipo SI ) 

10-T-p-T 2T v K ' 



, , 1 2 

(b = 1 + e 

12 • T - p • T 5T 



-1.5(T-0.2) 2 



Depositos aluvionales (suelo tipo S2) 



(|) = 1 + 




3 -T 




Suelos blandos (suelo tipo S3) 



Siendo T„ el periodo caracteristico del suelo, SI corresponde a suelos rocosos, S2 a 
depositos aluvionales y S3 a suelos blandos. 

Ordaz y Perez Rocha (1998) 

Proponen espectros de calculados a partir de espectros elasticos de 

desplazamiento y de velocidades escalados mediante el desplazamiento y velocidad 
maxima del terreno alcanzada. La novedad de este enfoque es que permite incorporar 
directamente el efecto de sitio sin tener expresiones que dependan explicitamente del 
periodo estructural. La expresion propuesta es la siguiente: 





P 



( 13 ) 





Factores de reduccion de respuesta estado del arte y estudio comparative entre codigos 83 



Es de hacer notar sin embargo, que el periodo queda incluido de forma implicita en los 
espectros de velocidad V(t) y de desplazamiento D(t). 



Lee, Han y Oh (1999) 

Utilizaron diferentes modelos histereticos poligonales que les permitieron establecer 
el grado de influencia del tipo de degradacion en los valores de R g de sistemas de un 

grado de libertad. La expresion depende del periodo de la estructura y de la ductilidad global 
esperada: 

R^A 0 [l-e- B « T ] 

A 0 = 0.99-p + 0.15 (14) 

B 0 = 23.69 • pT 0 ' 83 



Chakraborti y Gupta (2004) 

Este trabajo presenta un enfoque similar al trabajo de Ordaz y Perez Rocha, pero, 
ademas de aplicar modelos que incorporan degradacion de la estructura, incorporan las 
caracteristicas de sitio mediante al escalar los espectros de diseno (aceleraciones) 
previamente normalizados respecto a la aceleracion pico alcanzada en el terreno. En este 
caso los espectros de R^ si poseen una dependencia explicita respecto al periodo 
estructural: 



R u =P( T ) 



psa(t) 

PGA 



«(T) 



+ W 



V l0y 



(15) 



Vielma et al. (2004) 

Consideran el comportamiento histeretico de suelo y estructura, utilizando modelos 
no lineales con degradacion de rigidez, deterioro de resistencia y estrangulamiento propios 
de tipos estructurales prescritos en codigos sismorresistentes. Los factores R g propuestos 

se calculan a partir de espectros promedio para cada combinacion de casos especificos, 
obteniendo espectros idealizados bilineales dependientes de la ductilidad global, el periodo 
estructural y de parametros de ajuste que incluyen la respuesta degradada y el efecto de 
sitio para suelos que van desde los mas rigidos hasta los mas blandos. 

T <T g 

R g =l + p.T/(a.T g ) 

m m ( 16 ) 

g 

R .. = 1 + + bT 

r* /a 



Dolsek y Fajfar (2004) 

Investigan sobre la influencia de la degradacion de estructuras de hormigon armado y 
mampostena, realizando un estudio parametrico aplicando un modelo poligonal 
representative de la degradacion combinada de ambos materiales. Obtienen expresiones 
para los factores de reduccion haciendo enfasis en dos aspectos, el primero que la 
influencia de la degradaciones importante en la determinacion de factores de reduccion, 
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segundo, que los espectros inelasticos de diseno adquieren una forma diferente a la clasica 
en los codigos de diseno, consistente en tres ramas donde la rama de aceleraciones se 
mantiene constante, al tener valores de factores de reduccion variables esta rama no es 
horizontal sino que presenta pendiente negativa y tercero, resaltan la importancia de 
compatibilizar los espectros de diseno con los espectros de reduccion para obtener rangos 
de perfodos para el cambio de las expresiones de reduccion. La expresion basica propuesta 
por Dolsek y Fajfar es: 



R ^ =c(p-p 0 ) + R o ( 17 ) 

Donde R 0 es la reduccion que depende del perfodo y de la (reducida) ductilidad aportada 
por la mamposterfa, c es un coeficiente que se obtiene mediante el perfodo esquina, el 
grado de deterioro de resistencia esperado y el perfodo fundamental de la estructura, p es 

la ductilidad global y p 0 es la ductilidad complementaria que es tambien funcion del aporte 
de ductilidad de la mampostena. 

4.2 Factores de reduccion por sobre resistencia 

La gran cantidad de factores que intervienen en la sobre resistencia hacen que esta sea 
diffcil de cuantificar a nivel de cada elemento estructural evitando que esta pueda ser 
incorporada de forma racional en los codigos de diseno sismorresistentes a pesar de que se 
le considera muy importante en el mejoramiento de la respuesta de los edificios; por esto es 
mas conveniente considerar el efecto de la sobre resistencia sobre la respuesta global de la 
estructura mediante la utilizacion de las relaciones de fuerza-deformacion obtenidas del 
analisis estatico inelastico (pushover). 

La incorporacion del factor de sobre resistencia no se realiza de forma explfcita en los 
codigos de diseno debido a esta gran cantidad de factores; sin embargo un mejor 
conocimiento del fenomeno permitirfa su aplicacion desacoplada del factor de reduccion de 
respuesta R . 

El factor de sobrerresistencia R s se origina por la diferencia que existe entre el 
cortante en la base de diseno, V b , y el cortante en la base ultimo, V u . Esta diferencia 
resulta del hecho de proyectar para grupos de fuerzas mayoradas de forma estadfstica con 
base en la dispersion que las estas presentan (mayores a las resistidas por la estructura 
durante su vida util) conjuntamente con propiedades minoradas de los materiales. 
Adicionalmente, el diseno de los elementos sismorresistentes implica satisfacer condiciones 
de desplazamientos maximos (derivas de piso) y de detallado de partes de la estructura 
para mejorar su capacidad de disipacion ductil. 

En la Figura 1 se puede observar como se establecen las relaciones de resistencia. 
Primero, la curva de comportamiento permite definir un comportamiento elasto plastico 
idealizado, considerando el punto para el que se alcanza la maxima fuerza cortante en la 
base, claramente esta fuerza cortante maxima se alcanza en un punto en el que la 
estructura ha alcanzado el comportamiento plastico, permitiendo ubicar el segmento 
horizontal B-C. Seguidamente se define la recta de comportamiento elastico O-B uniendo el 
origen con el punto sobre la curva de comportamiento que corresponde al 75% del valor del 
cortante maximo idealizado, puede notarse que este punto presenta un cortante en la base 
superior al que se alcanza al ocurrir la primera plastificacion o formacion de la primera rotula 
plastica en algun punto de la estructura. Ademas este punto permite definir un 
desplazamiento ‘elastico’. 



Figura 1 Resistencias y desplazamientos caracterfsticos para el calculo de la 
sobrerresistencia y la ductilidad de desplazamiento. 





Factores de reduccion de respuesta estado del arte y estudio comparative entre codigos 85 



Tabla 2 Factores de reduccion por sobrerresistencia. 



Los aspectos anteriores conducen finalmente a secciones con resistencia superior a 
la calculada por simple resistencia. Investigaciones recientes han permitido determinar 
factores de sobrerresistencia (Bertero 1986) que demuestran la clara dispersion de los 
valores de R s para edificios porticados de hormigon armado, por lo que es necesario 
profundizar en este campo, (ver Tabla 2). De la tabla anterior se tiene que los valores de 
sobrerresistencia mas altos corresponden a los porticos en los que las cargas gravitatorias 
poseen una mayor influencia sobre el diseno de las secciones, por esta razon la 
sobrerresistencia de los porticos resistentes a momento alcanza valores mas altos que los 
porticos perimetrales, para los que las cargas gravitatorias son menores pero que 
mantienen las mismas fuerzas sismicas de diseno que los porticos resistentes a momento. 

Al establecer la relacion entre el diseno actual de edificios de hormigon armado, se 
Mega a alcanzar valores de los factores de resistencia R s que en todo caso presentan la 
misma dispersion de valores (ver la Tabla 3). 

Tabla 3 Factores de sobre resistencia, segun diferentes autores. 



La estimacion de los factores de reduccion por resistencia R s se complica por 
efectos como la resistencia actualizada de los materiales (una vez que se ha incursionado 
en comportamiento plastico o ha ocurrido dano), los efectos del confinamiento, presencia de 
elementos estructurales tales como losas de hormigon y elementos no estructurales no 
considerados de forma explicita en el analisis, haciendo necesario reconocer que no todos 
estos efectos contribuyen favorablemente al incremento del factor de sobrerresistencia. 
Entre los efectos que pueden tener contribucion desfavorable se incluyen el incremento de 
la sobrerresistencia en vigas puede dar origen al fallo fragil por cizallamiento, elementos de 
cerramiento pueden causar efectos de columna corta y las irregulares distribuciones de las 
alturas de las columnas y la discontinuidad en las lineas resistentes verticales pueden 
producir el efecto de piso blando (Park 1996). En la Tabla 4 se presenta una lista que 
contiene los factores anteriormente citados y la etapa del proyecto de la edificacion en la 
cual son incorporados. 
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Tabla 4 Factores que influyen en la sobrerresistencia. 



Causa de la sobrerresistencia 


Etapa del proyecto en la cual se 


Modelado 


Analisis 


Diseno 


Diferencia entre resistencia actual y nominal de los materiales 


X 






Diseno conservador y requerimientos por dutilidad 






X 


Factores de carga y multiples casos de carga combinados 




X 




Consideration de la torsion accidental 




X 




Provisiones para cumplir el estado limite de servicio (flechas) 






X 


Participation de los elementos no estructurales (losas) 


X 






Cumplimiento de las cuantias minimas de acero 






X 


Redundancia estructural 


X 






Endurecimiento resistente 


X 






Efectos del confinamiento 


X 






Analisis con base en el periodo elastico 




X 




Diseno controlado por rigidez 






X 



Jain y Navin (1995), estudiaron los valores de sobrerresistencia en edificios de 
hormigon armado, considerando la influencia de la intensidad sismica representado por la 
aceleracion de diseno basica correspondiente a diferentes zonas sismicas de la India y del 
numero de niveles del edificio. Los resultados de la Figura 2, en la que las zonas sismicas 
van desde la de menor amenaza (I) a la de mayor amenaza (V), muestran que la 
sobrerresistencia es mayor para edificios ubicados en zonas de baja sismicidad en los que 
predomina el diseno por cargas gravitacionales. Igualmente, debe resaltarse el hecho que 
para edificios de baja altura en los que las fuerzas de piso no adquieren valores tan altos 
como en los edificios de mayores masas sismicas, los factores de sobrerresistencia son 
mayores. 




0123456789 10 

Numero de niveles 

Figura 2 Factores de reduccion por sobrerresistencia determinados por Jain y Navin 

(1995). 



Mas recientemente, Massumi et al. (2004), retoman el estudio de la influencia del 
numero de niveles e incluyen el numero de vanos en el calculo de factores de 
sobrerresistencia obtenidos mediante los resultados de las relaciones de fuerza deformacion 
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resultantes del analisis estatico y dinamico no lineal. Los resultados obtenidos muestran que 
el numero de vanos o de lineas resistentes verticales son inversamente proporcionales a los 
valores de sobrerresistencia, lo que contradice los resultados obtenidos de forma 
probabilista por Bertero y Bertero (1999). 




0123456789 10 

Numero de niveles 

Figura 3 Factor de reduccion por sobrerresistencia, determinado por 
Massumi et al. (2004). 



En la Figura 3 se muestran factores de sobrerresistencia calculados en este trabajo, 
apreciando que los valores senalados son semejantes a los de la figura anterior, para el 
caso de edificios disenados para zonas de alta sismicidad. En esta misma figura es posible 
observar que los valores de sobrerresistencia varian muy poco con el metodo de analisis no 
lineal aplicado, sea este dinamico o estatico, conclusion similar a la de otros trabajos 
(Elnashai y Mwafi 2004). 

Otro de los factores a considerar al evaluar la sobrerresistencia global es el efecto 
producido por la ductilidad global de diseno. Las secciones y el armado resultantes del 
diseno de una estructura para un nivel alto de ductilidad son mayores que en el caso que la 
misma estructura se disenase para un nivel inferior de ductilidad. De esta forma es posible 
obtener valores de sobrerresistencia mayores para edificios de alta ductilidad que para 
edificios de baja ductilidad (Mwafi y Elnashai, 2002). Sin embargo el analisis de diversos 
niveles de ductilidad escapa de los alcances de este trabajo. 

Mwafi y Elnashai (2002), concluyen que los valores de la sobrerresistencia se deben 
fundamentalmente a tres fuentes principales: las caracteristicas resistentes reales de los 
materiales, la utilizacion del periodo elastico en el calculo de las demandas sismicas y el 
diseno de secciones de pilares utilizando la flexion biaxial. De acuerdo con estos factores, 
consideran que adoptar valores rmnimos de sobrerresistencia igual a 2, conduce a 
resultados conservadores. 



4.3 Factores de reduccion por redundancia estructural 



La redundancia estructural ha cobrado valor a partir del comportamiento de edificios 
de hormigon armado durante los sismos de Northridge (1994) y Kobe (1995). El termino 
‘redundancia’ dentro del contexto de la ingenieria sismorresistente no debe ser interpretado 
como superfluo o sobrante, Bertero y Bertero (1999). 
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De todos los factores el que menos ha sido estudiado es el factor de reduccion por 
redundancia, encontrando entre las investigaciones resultados que en algunos casos son 
contradictorios, a pesar de que todos coinciden en el efecto beneficioso que tiene la 
redundancia en la adecuada respuesta e incluso en la supervivencia de las estructuras 
sometidas a acciones sismicas severas. 

La mayoria de las referencias se fundamentan en la investigacion de Moses 1974, 
que realizo un estudio estadistico de la respuesta de edificios porticados frente a las 
acciones de viento, variando el numero de lineas resistentes verticales. Posteriormente, 
Bertero y Bertero 1999, retoman el tema adicionando al planteamiento estadistico la 
influencia del tipo de sismo, pulsativo o armonico, en la posible respuesta de edificios de 
multiples grados de libertad. Como ejemplo ilustrativo, se propuso el analisis de la respuesta 
de una estructura de un solo portico con otra de dos porticos, para este ultimo se establece 
la correlacion entre la calidad constructiva de cada linea resistente vertical independiente 
(linea de pilares) mediante un coeficiente de correlacion: 

(18) 

Pi -\/l+ Pv 

En la relacion anterior p, es el indice de seguridad para una linea resistente, P 2 es 
el indice de seguridad para dos lineas resistentes y p v es el coeficiente de correlacion de 
calidad constructiva, que resume la relacion entre la calidad de construccion existente entre 
dos lineas resistentes verticales de la misma estructura. 

Queda claro que para estructuras con un numero de lineas resistentes mayor, el 
coeficiente de correlacion de la calidad constructiva se reduce, incrementando la relacion de 
indice de seguridad de un edificio con varias lineas resistentes, respecto al edifico con una 
sol linea resistente. Valores considerados por Bertero y Bertero (1999) indican que, para 
que una estructura disponga de un factor de reduccion por redundancia al menos igual que 
la unidad, debera disponer de al menos 4 lineas resistentes verticales (ver la Tabla 4). 

Tabla 5 Factores de reduccion por redundancia (Bertero 1999). 



Las causas de esta ultima consideracion son las siguientes: 

■ Capacidad de rotacion en secciones en las cuales ya ha ocurrido la plastificacion 

■ Se evita el conocido efecto de tamano por concentracion de solicitaciones de 

grandes magnitudes en un numero limitado de secciones o elementos resistentes. 

La redundancia de edificios duales (porticados y con muros de corte) de hormigon 
armado, analizada siguiendo las directrices del trabajo anterior, concluye que el incremento 
de lineas resistentes verticales combinando las bondades del sistema dual, permite esperar 
incrementos del factor de resistencia del orden de 24%. 



5 LOS FACTORES DE REDUCCION DE RESPUESTA EN LAS NORMAS 
SISMORRESISTENTES 
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5.1 Factores de reduccion de respuesta en el pr EN-1 998-1 :2003(E) (EC-8) 



Al igual que en la mayoria de las normas de diseno sismorresistente, en el EC-8 se 
admite el analisis elastico de las estructuras, partiendo de la reduccion de las ordenadas del 
espectro de diseno elastico mediante la aplicacion de factores de reduccion. En la Figura 4 
se muestran los espectros elasticos de diseno correspondientes a la version 2003 del EC-8, 
para riesgo sismico alto (tipo 1). 

Los factores de reduccion de respuesta son identificados en el EC-8 mediante la 
letra q. En lo sucesivo, en este trabajo se utilizara la letra R para identificar los factores de 
reduccion, inclusive a los del EC-8. Estos factores de reduccion son aplicados a los 
espectros elasticos de diseno para obtener las aceleraciones de diseno y se expresados 
mediante la siguiente ecuacion: 



R = q 0 k w >l,5 (19) 

En esta ecuacion, q 0 es el factor basico de reduccion de respuesta que depende de 
la tipologia estructural predominante, de su regularidad en elevacion, asi como del nivel de 
ductilidad que se espera que la estructura sea capaz de desarrollar. Para dicha ductilidad 
existen dos niveles: DCM (clase de ductilidad media) y DCH (clase de ductilidad alta) y sus 
valores se muestran en la Tabla 6. El otro termino de la Ecuacion (19) es el factor k w que 
considera el modo de fallo predominante en estructuras con muros de cortante. 




Figura 4 Espectros elasticos de diseno, riesgo sismico alto (EC-8). 



Tabla 6 Valores sugeridos para el factor basico de reduccion de respuesta q 0 . 
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En la Tabla 6, el valor a,j es el coeficiente sismico correspondiente a la formacion 
de la primera rotula plastica en cualquiera de los elementos de la estructura. El valor a u es 
el coeficiente sismico requerido para que se produzca la inestabilidad de la estructura, en 
otras palabras, para que se haya alcanzado el fallo de la misma. La relacion de estos 
valores, a u /ai, describe el efecto de la redundancia estructural sobre el factor de reduccion 
de respuesta R . De acuerdo con las especificaciones de la norma EC-8, los valores y 
a u de la Tabla 6 se pueden obtener directamente de la curva de capacidad de la 
estructura, determinada mediante el analisis inelastico estatico. En este caso, (Xj es el valor 
de la ordenada del cortante en la base de disefio y a u es el cortante en la base inelastico, 

correspondiente al desplome ultimo, de colapso. En consecuencia, la relacion ajai equivale 
a un factor de reduccion por redundancia. Sin embargo, para los casos en los que no se 
realiza un analisis inelastico estatico, se recomienda aplicar los valores de la Tabla 7 
cuando la estructura es predominantemente de porticos o de porticos acoplados a muros de 
cortante. Cabe senalar que este factor multiplicador tiene en cuenta el efecto de la 
redundancia estructural, al considerarse valores mas altos a medida que el numero de 
vanos y el de niveles son mayores. Si el sistema estructural es predominantemente de 
muros de cortante, los factores de reduccion por resistencia R s , se obtienen de la Tabla 8. 

Los multiplicadores, adquieren valores entre 1,1 y 1,30, pero presentan una gran 
dispersion, por lo que se recomienda la aplicacion del analisis estatico no lineal para poder 
determinar de forma precisa los valores para cada edificio. 

Si el sistema estructural es predominantemente de muros de cortante, los 
multiplicadores se obtienen de la Tabla 8. 

Tabla 7 Valores sugeridos para el factor basico de reduccion de respuesta q 0 . 



Tabla 8 Valores sugeridos para los multiplicadores. 



En la Figura 5 se muestra la curva de capacidad obtenida del analisis estatico no 
lineal, a partir de la cual se obtienen las relaciones de coeficientes para calcular los factores 
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multiplicadores. En el eje de las abscisas se representan los desplomes del nivel de cubierta 
expresados en %, normalizados respecto a la altura total del edificio y en el eje de las 
ordenadas se representa el cortante en la base normalizado respecto al peso sismico del 
edificio. Los multiplicadores obtenidos aplicando el procedimiento en base al analisis 
inelastico estatico no superaran el valor de 1 ,5. 



El otro termino del factor de reduccion de respuesta tiene que ver con la respuesta 
de estructuras con predominio de muros de cortante. En este caso se aplica un factor de 
predominio del modo de fallo k w , que es menor o igual que la unidad, al factor basico q 0 y, 
por tanto penaliza el predominio del fallo por cortante. Las estructuras con muros de 
cortante tienen menor capacidad de desarrollar respuesta ductil en comparacion con las 
edificaciones porticadas. De esta manera, en la medida en que los muros son mas esbeltos, 
el valor de k w se aproxima a su cota inferior que es 1/3, sin llegar a alcanzarla ya que el 
EC-8 lo limita expresamente a 0,5. 




Efecto de la irregularidad estructural 

La irregularidad en planta es penalizada con factores multiplicadores iguales a 1, lo 
que equivale a no aplicarle ningun incremento al factor basico de reduccion por resistencia y 
redundancia. Sin embargo, en el articulado del EC-8 no se distingue entre los diferentes 
tipos de irregularidad en planta, pero que determinan de modo diferente los efectos 
torsionales. Ademas, los valores de q 0 pueden reducirse hasta un 20% si la estructura 
analizada presenta irregularidades en elevacion; el porcentaje de reduccion queda a juicio 
del proyectista, dependiendo del grado de irregularidad. De igual forma, se establece un 
incremento de los valores de q 0 cuando se disponga de un plan especial de calidad, que 
nunca superara el 20%. 

Factores de reduccion R . Edificios porticados 
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En las Tablas 9 y 10 se muestran los factores de reduccion para edificios porticados 
teniendo en cuenta los factores basicos de calculo asignados de acuerdo con el nivel de 
ductilidad esperado (alto o medio), la regularidad de ia estructura, la existencia de 
irregularidades en planta, en elevacion y una combinacion de ambas. Se aplican los valores 
extremos sugeridos en el EC-8, entendiendose que entre los maximos y los mmirnos se 
encuentra un rango de valores que, tal como se ha indicado, quedan a juicio del proyectista. 
Finalmente, se aplican incrementos a los factores de reduccion en virtud de la existencia de 
planes especiales de control de calidad en la elaboracion del proyecto. En las Tablas 9 a la 
13, las siglas CPC y SPC significan que el proyecto cuenta o no, respectivamente, con un 
plan de control de calidad. 



Tabla 9 Factores de reduccion para estructuras porticadas con nivel de ductilidad alto. 



Irregulares 











regulares 


en planta 


en elevacion 








Factor basico 
qO 


Factor de reduccion 
por redundancia 

aval 


Factor por modo 
predominante de 
fallo Kw 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


un solo nivel 


4,50 


1,10 


1,00 


4,95 


5,94 


4,50 


5,40 


4,13 


4,95 


3,75 


4,50 


multiples niveles 
un solo vano 


4,50 


1,20 


1,00 


5,40 


6,48 


4,50 


5,40 


4,50 


5,40 


3,75 


4,50 


multiples niveles 
multiples vanos 


4,50 


1,30 


1,00 


5,85 


7,02 


4,50 


5,40 


4,88 


5,85 


3,75 


4,50 



SPC: Edificio sin proyecto de control de calidad 
CPC: Edificio con proyecto de control de calidad 



Tabla 10 Factores de reduccion para estructuras porticadas con nivel de ductilidad bajo. 



Irregulares 











regulares 


en planta 


en elevacion 








Factor basico 
qO 


Factor de reduccion 
por redundancia 
ctu/al 


Factor por modo 
predominante de 
fallo Kw 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


un solo nivel 


3,00 


1,10 


1,00 


3,30 


3,96 


3,00 


3,60 


2,75 


3,30 


2,50 


3,00 


multiples niveles 
un solo vano 


3,00 


1,20 


1,00 


3,60 


4,32 


3,00 


3,60 


3,00 


3,60 


2,50 


3,00 


multiples niveles 
multiples vanos 


3,00 


1,30 


1,00 


3,90 


4,68 


3,00 


3,60 


3,25 


3,90 


2,50 


3,00 



SPC: Edificio sin proyecto de control de calidad 
CPC: Edificio con proyecto de control de calidad 



Factores de reduccion R . Edificaciones con muros de cortante 

Para este tipo estructural, que comprende una amplia gama de configuraciones, 
inclusive combinaciones con estructuras porticadas, se presenta una variante que consiste, 
en la aplicacion de un factor corrector que varia en funcion del modo de fallo predominante. 
En las Tablas 11 a la 13, se muestran los valores correspondientes a estructuras con muros 
de cortante para ductilidades bajas y altas, adoptandose el valor de k w = 1 para los casos 
con predominio de fallo por flexion y de k w = 0,5 que es la cota minima norma para 
edificios con fallos por cortante. 
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Tabla 11 Factores de reduccion para estructuras con muros de cortante con nivel de ductilidad 

bajo. 



Irregulares 











regulares 


en planta 


en elevacion 








Factor basico 
qO 


Factor de reduccion 
por redundancia 
ctu/cd 


Factor por modo 
predominante de 
fallo Kw 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


solo dos muros 
desacoplados 


4,00 


1,00 


1,00 


4,00 


4,80 


4,00 


4,80 


3,33 


4,00 


3,33 


4,00 


otros sistemas 
desacoplados de 
muros 


4,00 


1,10 


1,00 


4,40 


5,28 


4,00 


4,80 


3,67 


4,40 


3,33 


4,00 


sistemas duales o 
acoplados a muros 


4,00 


1,20 


1,00 


4,80 


5,76 


4,00 


4,80 


4,00 


4,80 


3,33 


4,00 



SPC: Edificio sin proyecto de control de calidad 
CPC: Edificio con proyecto de control de calidad 



Tabla 12 Factores de reduccion para estructuras con muros de cortante, con nivel de ductilidad 

alto, modo de fallo predominante por cortante. 



Irregulares 











regulares 


en planta 


en elevacion 








Factor basico 
qO 


Factor de reduccion 
por redundancia 
ctu/al 


Factor por modo 
predominante de 
fallo Kw 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


solo dos muros 
desacoplados 


4,00 


1,00 


0,50 


2,00 


2,40 


2,00 


2,40 


1,67 


2,00 


1,67 


2,00 


otros sistemas 
desacoplados de 
muros 


4,00 


1,10 


0,50 


2,20 


2,64 


2,00 


2,40 


1,83 


2,20 


1,67 


2,00 


sistemas duales o 
acoplados a muros 


4,00 


1,20 


0,50 


2,40 


2,88 


2,00 


2,40 


2,00 


2,40 


1,67 


2,00 



SPC: Edificio sin proyecto de control de calidad 
CPC: Edificio con proyecto de control de calidad 



Tabla 13 Factores de reduccion para estructuras con muros de cortante, con nivel de ductilidad 

alto, modo de fallo predominante por flexion. 



Irregulares 











regulares 


en planta 


en elevacion 








Factor basico 
qO 


Factor de reduccion 
por redundancia 
ctu/al 


Factor por modo 
predominante de 
fallo Kw 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


q 

(SPC) 


q 

(CPC) 


solo dos muros 
desacoplados 


3,00 


1,00 


1,00 


3,00 


3,60 


3,00 


3,60 


2,50 


3,00 


2,50 


3,00 


otros sistemas 
desacoplados de 
muros 


3,00 


1,10 


1,00 


3,30 


3,96 


3,00 


3,60 


2,75 


3,30 


2,50 


3,00 


sistemas duales o 
acoplados a muros 


3,00 


1,20 


1,00 


3,60 


4,32 


3,00 


3,60 


3,00 


3,60 


2,50 


3,00 



SPC: Edificio sin proyecto de control de calidad 
CPC: Edificio con proyecto de control de calidad 



La aplicacion de los factores de reduccion del EC-8 no conduce a valores unicos 
para una ductilidad esperada, sino que proporciona un rango de valores que dependen, por 
un lado, de la capacidad para disipar energia que poseen las estructuras y de su 
sobrerresistencia y, por el otro, de los controles de disefio por parte del proyectista. A la 
larga, dichos controles aseguran el cumplimiento de las prescripciones de la norma, lo que 
conduce al incremento de la capacidad disipativa de la estructura. 
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Los factores tipicos de reduccion contemplados en el EC-8 se muestran en la Figura 
6, para estructuras porticadas y una combinacion de suelo rocoso (tipo A). Se observa que 
un espectro de factores de reduccion presenta dos ramas: una rama creciente y constante, 
y una que se reduce una vez alcanzado el periodo maximo de aceleracion constante del 
espectro elastico de diseno. Ademas, el EC-8 define un valor minimo de aceleracion de 
diseno para los periodos altos, razon por la que se observa un cambio de pendiente en la 
rama descendente del espectro de factores de reduccion para un periodo de 2 segundos. 

5.2 Factores de reduccion de respuesta en la Norma NCSE-02 

La norma sismorresistente espanola vigente NCSE-02 permite el analisis elastico de 
estructuras mediante la reduccion de los valores espectrales al aplicar coeficientes de 
respuesta |3 . El coeficiente de respuesta es funcion de la ductilidad esperada y del 
porcentaje de amortiguamiento critico que la misma norma prescribe para los diferentes 
materiales que constituyen el sistema estructural. 

P = - (20) 

Donde v es el factor de modificacion en funcion del amortiguamiento, calculado 
conforme a la expresion (20) y p es el coeficiente de comportamiento por ductilidad. El 
factor de modificacion por amortiguamiento considera la reduccion directa del espectro 
elastico en funcion de la capacidad disipativa que posee la estructura, segun la relacion 
entre el amortiguamiento de 5%, que usualmente se encuentra prescrito en otros codigos de 
diseno y el amortiguamiento real de la estructura, D. . 




Figura 6 Espectros de factores de reduccion. 

La norma establece los valores del amortiguamiento tipicos para estructuras 
usuales, los cuales no solo son aplicados a la reduccion de espectros elasticos, sino 
tambien al incremento de sus ordenadas. Los valores del coeficiente de comportamiento por 
ductilidad se presentan en la Tabla 14. 
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Si bien la norma NCSE-02 no define directamente el uso de factores de reduccion, 
estos se expresan aqui como el inverso de los coeficientes de respuesta [3 , para permitir la 

comparacion con otras normas. De esta forma, se obtienen los valores de R para cada uno 
de los niveles de ductilidad, que dependen fundamentalmente de la capacidad de las 
estructuras de disipar energia frente a la accion de cargas ciclicas, con elementos 
estructurales expresamente disenados para tal fin y acoplados en las dos direcciones 
principales de la estructura. Los valores de los factores de reduccion de respuesta se 
muestran en las Tablas 15 a la 17, para ductilidad alta, intermedia y baja, respectivamente. 



Tabla 14 Parametros del factor de reduccion previstos en la norma NCSE-02. 



Tabla 15 Factores de reduccion norma NCSE-02, ductilidad alta. 



Tipo de estructura 


Compartimentacion 


Ductilidad 


Amortiguamiento 

V 


Multiplicador 

P 


Factor reductor 

R 


Estructura de 
hormigon armado o 
acero laminado 


Diafana 


4,00 


1,09 


0,27 


3,66 


Compartimentada 


4,00 


1,00 


0,25 


4,00 



Tabla 16 Factor de reduccion norma NCSE-02, ductilidad intermedia. 



Tipo de estructura 


Compartimentacion 


Ductilidad 

F 


Amortiguamiento 

V 


Multiplicador 

P 


Factor reductor 

R 


Estructura de 
hormigon armado o 
acero laminado 


Diafana 


3,00 


1,09 


0,36 


2,74 


Compartimentada 


3,00 


1,00 


0,33 


3,00 



Tabla 17 Factor de reduccion norma NCSE-02, ductilidad baja. 








































So/a 
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Tipo de estructura 


Compartimentacion 


Ductilidad 

M 


Amortiguamiento 

V 


Multiplicador 

P 


Factor reductor 

R 


Estructura de 
hormigon armado o 
acero laminado 


Diafana 


2,00 


1,09 


0,55 


1,83 


Compartimentada 


2,00 


1,00 


0,50 


2,00 


Estructura con 
muros de corte 


compartimentada 


2,00 


0,93 


0,46 


2,15 



En la Figura 7 se muestran los espectros correspondientes a un suelo granular de 
compacidad media (clasificado como del tipo III), reducidos por el factor de estructura 
porticada compartimentada, para un nivel de ductilidad alto. Cabe indicar que tanto los 
espectros elasticos como los inelasticos muestran ramas de aceleracion constante para 
distintos intervalos de periodos. Esto se debe a la aplicacion del factor de contribucion k 
contemplado en la norma NCSE-02. 




T(s) 

Figura 7 Espectros de diseno, elastico e inelasticos, para sistemas porticados con alta 

ductilidad, sobre suelo tipo III. 



Finalmente, en la Figura 8 se ha representado el factor de reduccion de respuesta 
tipico aplicado en la norma NCSE-02. Se distinguen las dos ramas: la primera variable 
segun el periodo T y la segunda tiene valores constantes e iguales al valor de la ductilidad 

p. 

5.3 Factores de reduccion de respuesta del UBC 97 



En el Uniform Building Code (UBC 97) se prescriben valores de los factores de 
reduccion de respuesta que dependen del sistema resistente a fuerzas laterales 
predominante y de las caracteristicas del material utilizado en el diseno. El nivel de diseno 
queda implicitamente definido al hacerse referenda a la categoria de “porticos especiales”, 
que corresponde al cumplimiento de las recomendaciones de dimensionado y los detalles 
de armado adecuados para zonas sismicas. La 0 muestra aquellos valores de R que 
interesan en virtud de las estructuras que se consideran en la presente monografia. 
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En la Tabla 18 se muestran tambien los valores de los factores de sobrerresistencia 
estructural,Q 0 , que en el UBC-97 se asignan a cada tipo especffico de estructura. Este 
factor es aplicado para amplificar las componentes horizontales de las cargas sismicas, E h , 
determinada a partir del cortante en la base, mediante la expresion: 

E m =^ 0 E h (21) 

Donde E m es la maxima fuerza sismica estimada. 

Puede observarse que el UBC 97 considera valores de reduccion evidentemente superiores 
a los del resto de las normas de diseno sismorresistente analizadas. Sin embargo, esta 
diferencia es solo aparente, al considerarse en la norma la posibilidad de incrementar las 
fuerzas de diseno teniendo en cuenta la redundancia estructural mediante la siguiente 
expresion: 



1 < p = 2 - 




<1.5 



( 22 ) 




Figura 8 Espectro tipico de factor de reduccion norma NCSE-02, suelos tipo III. 



Tabla 18 Factores de reduccion de respuesta y de sobre rresistencia, UBC-97. 





R 


Qo 


1. Porticos especiales resistentes a momentos 
a. Acero 


8,5 


2,8 


b. Hormigon 


8,5 


2,8 


2. Porticos de mampostena resistentes a momentos 


6,5 


2,8 


3. Portico intermedio resistente a momenta 


5,5 


2,8 


4. Porticos ordinarios resistentes a momentos 
a. Acero 


4,5 


2,8 


b. Hormigon 


3,5 


2,8 


5. Porticos de celosias de acero especiales resitentes a momentos 


6,5 


2,8 



En la Ecuacion (22), A B representa el area de la planta expresada en metros 
cuadrados, r max es la relacion entre el esfuerzo cortante maximo en los pilares y el cortante 
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total de piso de la estructura entera. Esta expresion engloba lo que se puede interpretar 
como una especie de densidad de los pilares que, en el caso de estructuras porticadas 
resistentes a momentos, son los encargados de resistir las fuerzas cortantes inducidas por 
las cargas horizontales de los sismos. La forma de aplicar la ecuacion anterior en el calculo 
de las fuerzas sismicas consiste en calcular las fuerzas horizontales de diseno E h 
mediante la ecuacion. 

E = pE h +E v (23) 



Donde E v representa la componente vertical de la accion sismica. 



5.4 Factores de reduccion norma Covenin 1756-98 

Los factores de reduccion son prescritos de forma directa en la normativa 
venezolana Covenin 1756-98 y sus valores dependen de dos caracteristicas, la primera 
vincula los niveles de reduccion de ordenadas espectrales con el nivel de diseno aplicado, 
obligando esto al cumplimiento de ciertas condiciones que garanticen la ductilidad de el las 
estructuras. La segunda caracteristica tiene que ver con la tipologia estructural utilizada, 
correspondiendo los niveles mas altos de reduccion de la respuesta a las estructuras 
porticadas (Tipo I), hasta las estructuras incapaces de tener un comportamiento asimilable 
al de las que responden como diafragma rigido. 



Los valores maximos de los factores de reduccion son los que en anteriores 
versiones de la normativa correspondian a la ductilidad, sin que intervengan en su 
determinacion factores adicionales, tales como la sobrerresistencia o el amortiguamiento, 
considerado en la determinacion de los espectros elasticos de diseno. 

La presentacion de los valores maximos se realiza para tres familias de estructuras 
que dependen del material: estructuras de hormigon armado, de acero y estructuras mixtas 
de hormigon armado y acero. En cuanto a la influencia de las posibles irregularidades tanto 
en planta como en elevacion, la normativa aconseja la aplicacion global de factores de 
reduccion iguales al 75% de los senalados para estructuras regulares, sin hacer ninguna 
valoracion sobre el tipo de irregularidad y sus posibles consecuencias sobre la respuesta 
ductil de la estructura. En las Tablas 19 a la 21, se incluyen los valores para estos tres tipo 
de estructuras. 



Tabla 19 Factor de reduccion normativa Covenin 1756-98, estructuras de hormigon armado. 



Tabla 20 Factor de reduccion normativa Covenin 1756-98, estructuras de acero. 



En estas tablas aparecen los valores de los factores de reduccion R por tipo de 
estructura, variando de las mas ductiles (tipo I) para las que la estructura predominante 
porticada, a las menos ductiles (tipo IV) dentro de las que se encuentran las estructuras 
sustentadas por una sola columna, sin diafragma rigido o con forjados con vigas planas. 
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Tambien dependen del nivel de diseno, asignando valores de R superiores para el nivel de 
diseno ND3 que impone un detallado mas exigente de secciones y uniones, y valores 
inferiores para el nivel de diseno sin consideraciones sismorresistentes (ND1). 



Tabla 21 Factor de reduccion normativa Covenin 1756-98, estructuras mixtas hormigon 

armado- acero. 



5.5 Factores de reduccion de respuesta del IBC 2003 

El Internacional Building Code (IBC 2003), presenta pocas diferencias respecto al 
UBC 97 en lo referente al calculo de valores de los factores de reduccion de respuesta; 
incluso el formato de presentacion es similar en cuanto a los casos de diseno, sistema 
estructural y tipo de material utilizado. Los valores mencionados se incluyen en la Tabla 22, 
en la que se observa que los valores mas altos corresponden a las estructuras cuya 
tipologia permite desarrollar un comportamiento ductil. Los factores de reduccion pueden 
ser modificados mediante el factor de importancia del edificio, que incrementan las fuerzas 
de diseno en el caso de los edificios considerados como esenciales. 



Tabla 22 Factores de reduccion de respuesta y de sobrerresis tencia IBC 2000. 





R 


Qo 


Porticos especiales de acero resistentes a momentos 


8 


3 


Porticos de celosias de acero especiales resistentes a momentos 


7 


3 


Porticos intermedios de acero resistentes a momentos 


4,5 


3 


Porticos especiales de hormigon resistentes a momentos 


3,5 


3 


Porticos intermedios de hormigon armado resistentes a momentos 


5 


5 


Porticos ordinarios de hormigon armado resistentes a momentos 


3 


3 


Porticos especiales de compuestos resistentes a momentos 


8 


8 


Porticos intermedios de compuestos resistentes a momentos 


5 


5 


Porticos de compuestos parcialmente restringidos resistentes a momentos 


6 


6 


Porticos ordinarios de compuestos resistentes a momentos 


3 


3 


Porticos de muros de mamposteria 


5,5 


5,5 



Tanto en la norma UBC-97 como en la IBC-2003, la sobrerresistencia y la redundancia se 
consideran en forma explicita como factores de amplificacion de las fuerzas sismicas de 
diseno. 



5.6 Factores de reduccion de respuesta en la NZS 1170.0:2002 

La norma neozelandesa NZS 1170.0:2002 utiliza factores de reduccion referidos a 
valores de ductilidad, siendo esta ultima fijada en funcion del tipo de estructura resistente 
predominante y del tipo de material que la compone. Los valores de ductilidad son 
directamente aplicados en la reduccion de las fuerzas de piso. En la Tabla 23 se muestra el 
resumen de los factores de reduccion para tres niveles de ductilidad: estructuras que se 
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espera que tengan una respuesta elastica, estructuras con una ductilidad intermedia y 
estructuras que, bajo condiciones especiales de proyecto, son capaces de desarrollar 
grandes deformaciones evitando el colapso por cargas verticales. 



Tabla 23 Factores de reduccion de respuesta NZS 1170.0:2002. 





Material predominante en la estructura 


Acero 

estructural 


Hormigon 

armado 


Hormigon 

pretensado 


Mamposteria 

reforzada 


Estructuras con respuesta elastica 


1,25 


1,25 


1,00 


1,25 


Estructuras con ductilidad limitada 






Porticos diagonalizados 




Cedencia en tension y compresion 


3,00 


- 


- 


- 


Cedencia en compresion solamente 


3,00 


- 


- 


- 


Porticos resistentes a momentos 


3,00 


3,00 


2,00 


2,00 


Muros 


3,00 


3,00 


- 


2,00 


Muros cargados en voladizo 


- 


2,00 


- 


2,00 


| Estructuras ductiles 






Porticos diagonalizados 


6,00 


- 


- 


- 


Porticos resistentes a momentos 


6,00 


6,00 


5,00 


4,00 


Muros 


- 




- 


4,00 


Porticos con diagonales excentricas 


6,00 




- 


- 



5.7 Factores de reduccion de respuesta del Japanese Building Standard Law 



La norma japonesa presenta como novedad que los factores de reduccion de 
respuesta se aplican como los inversos de los factores de ductilidad. Ademas es posible 
aplicar diferentes valores dependiendo del grado de ductilidad que sea preciso desarrollar 
en un nivel dado. 

La norma permite la aplicacion de valores que dependen del material que constituye 
la estructura, asignandole mayores valores de reduccion a las estructuras de acero, en 
comparacion con las estructuras de hormigon armado. En las Tablas 23 y 24 se muestran 
los valores sugeridos, que tambien dependen de la tipologia estructural predominante en el 
edificio, asi como tambien del nivel de ductilidad que el proyectista dispone para los 
elementos del nivel. 

Tabla 24 Factores de reduccion de respuesta para edificios de acero. 



Tabla 25 Factores de reduccion edificios de hormigon armado . 



Edificios de hormigon armado 


Ductilidad 


Portico resistente 
a momentos 


Otro tipo de 
portico 


Portico con barras 
a compresion 


Excelente 


0,30 


0,35 


0,40 


Buena 


0,35 


0,40 


0,45 


Regular 


0,40 


0,45 


0,50 


Pobre 


0,45 


0,50 


0,55 



Los factores mostrados indican que la norma japonesa considera el diseno para 
ductilidades bajas, lo que a la larga termina condicionando a que el edificio tenga un periodo 
mas bajo que los edificios disenados mediante otras normas sismorresistentes. Tambien es 













































Factores de reduccion de respuesta estado del arte y estudio comparative entre codigos 101 



necesario indicar que la norma contempla la penalizacion de la reduccion mediante la 
aplicacion de factores que valoran la regularidad de la estructura. 



6 COMPARACION DE LOS FACTORES DE DIFERENTES NORMAS 
SISMORRESISTENTES 



El primer aspecto a resaltar es que solo en el EC-8 existe una referenda explicita a 
un factor basico de ductilidad, que depende de la ductilidad global que puedan desarrollar 
las estructuras correctamente disenadas. En las normas NCSE-02 y Covenin 1756-98 los 
factores de reduccion tienen un componente de reduccion que depende exclusivamente de 
la ductilidad R^, para la cual se fijan niveles de diseno que no son mas que un marco de 

especificaciones minimas que se deben cumplir en la etapa de proyecto. Los dos factores 
que completan los factores de reduccion R , que son dependientes de la sobrerresistencia y 
del amortiguamiento, se incorporan en el EC-8 al considerar coeficientes multiplicadores 
que varian de acuerdo con la redundancia estructural R sr y de los materiales que integran 
la estructura R;- , esto ultimo al diferenciar entre los coeficientes multiplicadores asignados a 

estructuras de acero, hormigon armado y muros de corte, aunque es necesario aclarar que 
no se realiza una referencia directa al factor R;. en la norma. 

La NCSE-02 es la unica de las normas revisadas que si presenta una referencia 
directa a la influencia del amortiguamiento en el calculo de los factores de reduccion. Sin 
embargo recomienda el valor promedio para el desempeno de las diferentes estructuras, 
clasificadas de acuerdo con sus materiales, sin hacer referencia a la consideracion de 
valores que dependan de la influencia de elementos no estructurales tales como la densidad 
y distribucion de la mamposteria que se ha demostrado altera los valores esperados del 
amortiguamiento de las estructuras. En la Norma Covenin 1756-98 no se senala la 
participacion del amortiguamiento estructural en la determinacion de los factores de 
reduccion R ; solo es posible seleccionar valores de R que son ligeramente superiores 
para el caso de las estructuras porticadas de acero. 



La consideracion de la amplificacion que experimenta el movimiento sismico cuando 
atraviesa distintos tipos de suelos se encuentra de forma explicita en el EC-8 y en la Norma 
Covenin 1756-98, al adoptar ambas valores de factores de amplificacion dinamica |3 que 
dependen de las caracteristicas geotecnicas, generalmente representadas por la velocidad 
de ondas de cizallamiento y por los espesores de los estratos. La NCSE-02 presenta 
valores maximos como consecuencia de la amplificacion, similares para los suelos 
contemplados en dicha norma, ajustando por esta causa el rango de valores de los periodos 
para los que se tiene aceleracion constante, que se ensancha a medida que el suelo es 
menos rigido, (ver Figura 7). 

En cuanto a los valores de los factores R , es interesante notar que solo el EC-8 
presenta limitaciones a las aceleraciones para las estructuras con periodos largos 
(superiores a los 2 s) limitando explicitamente la aceleracion de diseno inelastico a un valor 
minimo de 0,2 veces la aceleracion espectral normalizada. Este ultimo aspecto hace que el 
espectro de factores de reduccion tenga cuatro ramas bien diferenciadas: la de los periodos 
cortos, conforme a la regia de igual energia; la de los periodos intermedios con 
aceleraciones constantes; una tercera rama decreciente, proporcional a la rama decreciente 
del espectro elastico; y una ultima rama que no decrece de la misma forma que la tercera, 
consecuencia de la cota inferior de las aceleraciones de diseno. 
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El perfodo caracteristico para el cual se produce la transicion entre la rama 
creciente dependiente del perfodo a la rama horizontal dependiente de la ductilidad, 
presenta cierta uniformidad en los tres codigos analizados. Estos valores se encuentran 
entre 0,2 s y 0,4 s, que corresponden a perfodos tfpicos obtenidos de los espectros de 
respuesta obtenidos con registros sobre suelos rocosos o muy duros. Sin embargo, 
recientes trabajos han demostrado la dependencia del perfodo caracteristico respecto al 
perfodo del suelo. 

Finalmente, se reconoce la importancia que se da en las normativas analizadas a la 
relacion que existe entre el valor de la ductilidad de desplazamiento y el nivel de diseno que 
se aplica a la estructura. Por ejemplo, la norma NCSE-02, por permitir ciertas 
irregularidades de la estructura y detalles de armado menos rigurosos para las secciones y 
los nudos, establece niveles maximos de ductilidad que se encuentran por debajo de los 
especificados en el EC-8, similar caracterfstica se encuentra en la norma japonesa, en la 
que se proponen factores de reduccion que dependen de la ductilidad de los elementos 
estructurales del nivel analizado del edificio. 

En casi todas las normas revisadas se ha podido verificar que los factores de 
reduccion de respuesta correspondientes a estructuras de acero, son mayores que los 
correspondientes a las estructuras de hormigon armado. Algunas normas incluso llegan a 
contemplar valores para los sistemas estructurales usuales, tal es el caso del EC-8, 
Covenin-1 756-98, NZBC, IBC-97 y UBC-2003. 

Se comienzan a incorporar aspectos sugeridos en estudios recientes, tales como la 
incorporacion de los factores de redundancia estructural. Ejemplo de esto son los factores 
que incrementan de los factores de reduccion del EC-8, cuando las estructuras cumplen con 
ciertas caracterfsticas de redundancia de elementos resistentes verticales. Igual 
caracterfstica se observa en el IBC-97 y el UBC-2003, sin embargo su aplicacion es de 
caracter aproximado. 

La regularidad estructural es considerada directamente en el EC-8 y la norma 
japonesa, en el resto de las normas el beneficio de la regularidad solo se considera de 
forma indirecta. 



7 CONCLUSIONES 

Se verifica la incorporacion en los codigos de diseno sismorresistente de recientes 
avances en la determinacion de los factores de reduccion de respuesta. Dentro de las mas 
relevantes se encuentran la incorporacion de factores de reduccion dependientes de la 
ductilidad, factores de reduccion dependientes de la sobrerresistencia y de la redundancia 
estructural. Estos avances recientes requieren una mejora en su aplicacion, como por 
ejemplo, determinar valores apropiados para tipos estructurales, tipos de materiales o 
rangos de perfodos. 

Se ha constatado en la mayorfa de los codigos la aplicacion de factores de 
reduccion de respuesta mayores para edificios porticados de acero que para edificios 
porticados de hormigon armado. Sin embargo, se nota la ausencia de factores de reduccion 
recomendados para edificios mixtos de porticos y muros de cortante. En otros casos no se 
suministran factores de reduccion para edificios cuya estructura predominante sea 
unicamente de muros de cortante. 

Las modificaciones normativas dan indicios de que los factores de reduccion de 
respuesta han sido disminuidos progresivamente. Los factores de reduccion no superan el 
valor de 8,5 para el caso de estructuras porticadas. 

La racionalizacion de los factores de reduccion de respuesta, propuesta en el ATC- 
95, implica que el proyectista estructural se involucre en la determinacion del factor de 
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reduccion al aplicar al espectro de diseno, al tener la posibilidad de escoger los factores que 
se adapten a las caracteristicas del edificio que proyecta. 

Se incluyen penalizaciones a las irregularidades estructurales, ignoradas en las 
anteriores versiones de las normativas sismorresistentes en la determinacion de los factores 
de reduccion de respuesta. 

A pesar de la incorporacion de los factores de reduccion por sobrerresistencia y 
redundancia, es necesario hacer mas transparente la relacion entre estos y aspectos 
relevantes del proyecto tales como la regularidad, la tipologia estructural y el grado de 
hiperestatismo. 
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RESUMEN 

La relacion entre el desplazamiento maximo inelastico con relacion al desplazamiento 
maximo elastico en sistemas de un grado de libertad, es un parametro que tiene gran utilidad 
en la evaluacion rapida de la deriva maxima de piso, en el calculo del factor del factor de 
reduccion de las fuerzas sismicas con el cual se pasa del espectro elastico al inelastico y 
tambien para evaluar la respuesta dinamica en sistemas de multiples grados de libertad. 

En el presente trabajo se obtienen relaciones entre los desplazamientos maximos 
inelasticos con respecto a los desplazamientos maximos elasticos, considerando el tipo de 
suelo; en base 80 sismos registrados en el Ecuador y 112 sismos artificiales, los mismos que 
son compatibles con los espectros de diseno elastico del Codigo Ecuatoriano de la 
Construccion. 

Previamente se presentan dos trabajos anteriores realizados en el Centro de 
Investigaciones Cientificas de la Politecnica del Ejercito sobre esta tematica. En uno de ellos no 
se clasificaron los registros sismicos de acuerdo al tipo de suelo y en el otro la base de datos 
fue pequena. 

Se detalla con bastante detenimiento la forma como se obtienen las curvas de ajuste y 
se comparan los resultados con los dos trabajos desarrollados anteriormente. 



ABSTRACT 

The relationship between the maximum inelastic displacement and the maximum elastic 
displacement in single-degree of freedom systems is a very useful parameter for fast evaluation 
of story drift, to calculate the reduction factor for seismic forces switching from an elastic to an 
inelastic spectra, and to evaluate the dynamic response in multiple-degree of freedom systems. 

Relationships between maximum inelastic and maximum elastic displacements are 
obtained in this work, taking into account different soil types. Eighty seisms registered in 
Ecuador and 112 artificial seisms were used, being compatible to the elastic spectra design 
provided by the Ecuadorean Code for Construction. 

The results of two previous works on the same subject, done at the Research Center of 
the Polytechnic of the Army are also reported. In one of them the seisms were not classified 
according to the soil type, and in the other one the seism database was small. 
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The procedure to obtain curve fitting is thoroughly detailed, and the results obtained are 
compared to the two works developed previously. 

1 INTRODUCCION 

Existen tres aplicaciones muy importantes de la relacion entre el desplazamiento 
maximo inelastico con respecto al desplazamiento maximo elastico; la una es en la evaluacion 
rapida de la deriva maxima de pisos en edificios de hormigon armado; otra en el calculo del 
factor de reduccion de las fuerzas sismicas con el cual se pasa del espectro elastico al 
inelastico por comportamiento no lineal de la estructura y una tercera que no se comenta en 
este articulo para encontrar la respuesta no lineal en sistemas de multiples grados de libertad, 
empleando espectros. 

Existen algunas formulaciones para evaluar en forma rapida la deriva maxima de pisos, 
Gupta y Krawinkler (2000), Aguiar (2006) pero una de las primeras fue propuesta por Miranda 
(1 997, 1 999, 2001 ) y se la describe a continuacion: 



_ /?, A A A t 2 

r Art 2 N h 

Donde /?, es el valor de paso de un sistema de un grado de libertad a un sistema de 
multiples grados de libertad; /3 2 es un factor de amplificacion que permite determinar la 
distorsion maxima de entrepiso a partir de la distorsion global de la estructura. La distorsion 
global se define como la relacion entre el desplazamiento lateral maximo en el tope del edificio 
dividido para la altura total del edificio; /? 3 es un factor que permite calcular los 
desplazamientos laterales maximos con comportamiento inelastico a partir de los maximos 
desplazamientos laterales con comportamiento elastico ; jB 4 es un factor que sirve para 
determinar el cociente entre la distorsion maxima de entrepiso y la distorsion global pero 
calculado en una estructura con comportamiento elastico lineal con relacion a la misma relacion 
pero calculada con comportamiento inelastico; N es el numero de pisos de la estructura, h es 
la altura de un entrepiso; S a es la aceleracion espectral elastica. El parametro /? 3 es el que 
estamos analizando en este articulo y se define de la siguiente manera: 

^■iNE~Pi^d => ~ Pi 

^ d 

Donde A INE es el desplazamiento maximo inelastico y S d es el desplazamiento 
maximo elastico, en un sistema de un grado de libertad. 

La deriva maxima de piso es uno de los parametros que mejor relaciona el dano que se 
espera en una estructura. Bertero et al (1 991 ), Moehle (1 992 y 1 996), Rodriguez y Aristizabal 
(1999). La evaluacion rapida de la deriva maxima de piso tiene una gran aplicacion cuando se 
trata de definir las polizas de seguro de viviendas. 

La otra aplicacion se habfa indicado es en la evaluacion del factor de reduccion de las 
fuerzas sismicas R con el cual se pasa del espectro de diseno elastico al espectro de diseno 
inelastico. Existen varias formulaciones para este factor, una de ellas es la propuesta por 
Whittaker et al ( 1987, 1999) que se indica a continuacion: 

R = R m R s R c 

Siendo R es el factor de ductilidad, R s es el factor de resistencia y R £ es el factor 
de amortiguamiento. La relacion entre el desplazamiento maximo inelastico A , NE con el 
desplazamiento maximo elastico S d , el factor R y la ductilidad // viene dada por: 
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^ INE ~ .. 

R M 



A 



INE 



JL 

R, 



J d '-ft 

De tal manera que el contar con ecuaciones para la relacion A INE /S d tiene su 
aplicacion en el calculo en forma indirecta del parametro R . Es importante destacar que 
existen otras formulaciones para el factor de reduccion R como la del ATC-1995 en que 
cambian el factor de amortiguamiento R, de la ecuacion ( 3 ) por el factor de redundancia R r , 
quedando: 

R = R/i Rs Rr 



2 TRABAJOS REALIZADOS 

En el Centro de Investigaciones Cientificas de la Politecnica del Ejercito se han 
realizado tres trabajos en los cuales se determina la relacion A INE / S d . El primero de ellos fue 
realizado por Aguiar y Guerrero (2006) trabajando con 63 acelerogramas registrados en Sur 
America los mismos que se indican en la tabla 1, son eventos cuya aceleracion maxima del 
suelo es mayor a 0.1 g., siendo g la aceleracion de la gravedad. En este trabajo a la relacion 
A ink / S d se denomino /? 3 . 

En las figuras 1 y 2 se indican las curvas de valores medios que se hallaron en el 
estudio para ductilidades de dos, tres y cuatro, trabajando con un modelo elasto perfectamente 
plastico y con un modelo bilineal en el cual la rama post fluencia tiene una pendiente a . 
Siendo a la relacion entre la rigidez elastica y rigidez inelastica. 

La figura 1 corresponde al modelo elasto perfectamente plastico en el cual a - 0 y la 
figura 2 para el modelo bilineal con a = 0.05 . 

Para la determinacion de las curvas que mejor se ajustan a los datos indicados se 
considero la forma propuesta por Nassar y Krawinkler ( 1991) ecuaciones ( 6 ) y ( 7 ). 
Logicamente que los coeficientes “a” y “b" de estas ecuaciones se determinaron en el estudio y 
son los que muestran en la tabla 2. 



Pi = 



[c(//-l) + l] 1/c 



c(T,a ) 



T a b_ 
1 + T a T 



Los valores encontrados en el estudio por Aguiar y Guerrero (2006) fueron comparados 
con los obtenidos por Miranda (2000) quien trabajo con 263 acelerogramas registrados en 
suelo firme en los Estados Unidos de Norte America, la ecuacion a la que Mega, luego de 
trabajar con un modelo elasto perfectamente plastico, es: 



Pi = 





fl-.l 




1 + 


exp(-12T ju (l ' x ) 




U 7 






Estimacion de demandas inelasticas de desplazamiento... 



107 



Tambien se la comparo con la ecuacion propuesta por Nassar y Krawinkler (1991) los 
valores de “a” y “b” que ellos encontraron se indican en la tabla 3 y fueron hallados en base a 
15 sismos registrados en California en suelo firme y aluvial. 

Tabla 1 Sismos considerados en el estudio realizado por Aguiar y Guerrero (2006) 



ACELERACION 



CODIGO 

Oil 



ESTACION CHARACATO 



ESTACION CHARACATO 



ALBERTO GIESECKE MATTO 



ALBERTO GIESECKE MATTO 



ESTACION PARQUE DE LA RESERVA 



ESTACION PARQUE DE LA RESERVA 



ESTACION LA MOLINA 



CESAR VIZCARRA VARGAS 



CESAR VIZCARRA VARGAS 



ESTACION PARQUE DE LA RESERVA 



ESTACION PARQUE DE LA RESERVA 



ESTACION PARQUE DE LA RESERVA 



ESTACION SURCO 



ESTACION SURCO 



ESTACION ZARATE 



ESTACION ZARATE 



ILOCA 



ILOCA 



RAPEL 



PAPUDO 



LLAYLLAY 



LLAYLLAY 



SAN FELIPE 



FECHA MAGNITUD COMPONENTE 



13/06/2005 



17/10/1966 



17/10/1966 



09/11/1974 



lEcgiEMililH 



31/05/1970 
03/10/1974 ' 



03/10/1974 



03/10/1974 
03/10/1974 ' 



05/01/1974 



05/01/1974 



04/09/1985 



IliZHiEMESSlI 

gESHsMUi 



raviaagtta 



LLOLLEO 


CHILE 


LLOLLEO 


CHILE 


LLOLLEO (UNIVERSITY OF CHILE) 


CHILE 


LLOLLEO (UNIVERSITY OF CHILE) 


CHILE 


VINA DEL MAR (UNIVERSITY OF CHILE) 


CHILE 


VINA DEL MAR (UNIVERSITY OF CHILE) 


CHILE 


VALPARAISO EL ALMENDRAL, (UNIVERSITY OF CHILE) 


CHILE 


VALPARAISO EL ALMENDRAL, (UNIVERSITY OF CHILE) 


CHILE 


COLBUN 


CHILE 


ILOCA 


CHILE 


ILOCA 


CHILE 


MELIPILLA 


CHILE 


MELIPILLA 


CHILE 


PICHILEMU 


CHILE 


PICHILEMU 


CHILE 


SAN FERNANDO 


CHILE 


SAN FERNANDO 


CHILE 


SAN ISIDRO 


CHILE 


SAN ISIDRO 


CHILE 


ZAPALLAR 


CHILE 


ZAPALLAR 


CHILE 


ESCUELA EDMUNDO DE AMICIS-MENDOZA 


ARGENTINA 


ESCUELA EDMUNDO DE AMICIS-MENDOZA 


ARGENTINA 


MUNICIPALIDAD LAS HERAS - MENDOZA 


ARGENTINA 


MUNICIPALIDAD LAS HERAS - MENDOZA 


ARGENTINA 


MUNICIPALIDAD DE MAIPU - MENDOZA 


ARGENTINA 


ETNA S/N 2673 


ARGENTINA 


ETNA S/N 2673 


ARGENTINA 


SSA2EVT S/N 1428 


ARGENTINA 


SSA2EVT S/N 1428 


ARGENTINA 


ARMENIA (CARME) 


COLOMBIA 


ARMENIA (CARME) 


COLOMBIA 


ESTACION DE FILANDIA (CFLAN) 


COLOMBIA 


ESTACION DE FILANDIA (CFLAN) 


COLOMBIA 


ESTACION DE PEREIRA (CPER) 


COLOMBIA 


ESTACION DE PEREIRA (CPER) 


COLOMBIA 


ESTACION DE RICAURTE (CRICA) 


COLOMBIA 


ESTACION DE RICAURTE (CRICA) 


COLOMBIA 


ESTACION DE PEREIRA(Mazda) RISARALDA(CMAZP) 


COLOMBIA 



■liSHisareigi 

Ieehe^i 

■Liimra 

wawiwwmi 



08/04/1985 

03/03/1985 



IftKiWiliial 



03/03/1985 

03/03/1985 






HMIMIPB1.T 

PMiWiEMII 



EMSamai 



03/03/1 985 7.80 



mtmwiBBMl 

mimmM 



03/03/1985 

03/03/1985 



liMiSHESBlI 

PMiWiEMII 



26/01/1985 



26/01/1985 



26/01/1985 



26/01/1985 



26/01/1985 






H 1 



25/01/1999 
25/01/1999 ' 



25/01/1999 



25/01/1999 
25/01/1999 ' 



25/01/1999 



17/07/1998 



17/07/1998 



19/02/1997 



MAXIMA [gals] 

-138.5 



10° 


20 


o 

o 


43 
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1 

1 

1 

1 

X 

1 

1 


1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

X 

1 

1 






1 

1 

X 

1 

1 


1 

1 

1 


1 

1 

4 

1 

1 






1 

1 

i 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 


1 

1 

h 

1 

1 

1 

1 











1.00 1.50 2.00 

PERIODO (s) 

- DUCTILIDAD 2 DUCTILIDAD 3 — DUCTILIDAD 4 I 



Figura 1 Valores medios de /? 3 para a — 0 . 




1.00 1.50 2.00 

PERIODO (s) 

-DUCTILIDAD 2 DUCTILIDAD 3 DUCTILIDAD 4 I 



Figura 2 Valores medios de J3 3 para a — 0.05 



Tabla 2 Valores de “a” y “b” encontrados por Aguiar y Guerrer o (2006) 
a a b 

0.00 2.070 0.381 

0.05 1.247 0.248 
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Tabla 3 Valores de “a” y “£>” encontrados por Nassar y Krawinkler 



a 


a 


b 


0.00 


1.00 


0.42 


0.02 


1.00 


0.37 


0.10 


0.80 


0.29 



En la figura 3 se presentan tres graficos para ductilidades de 2, 3 y 4. En cada uno de 
ellos se indican los valores medios obtenidos con los sismos de la tabla 1, las curvas que se 
obtienen con las ecuaciones de ajuste de Aguiar y Guerrero (2006), Miranda (2000) y Nassar y 
Krawinkler (1991). Se aprecia que existe una muy buena aproximacion a pesar de que Aguiar y 
Guerrero (2006) trabajaron con sismos registrados en Sur America y los otros investigadores 
con sismos registrados en los Estados Unidos de Norte America; de tal manera que se podria 
indicar que la ubicacion geografica de los sismos no influye en la relacion A INE / S d . 

Los registros de la tabla 1 en un buen porcentaje fueron registrados en suelo firme por 
este motivo es que existe una buena correlacion con los otros trabajos ya que fueron en suelo 
similar. 



3 INFLUENCIA DEL TIPO DE SUELO 

El tipo de suelo tiene influencia en la determinacion del parametro /? 3 como lo han 

demostrado Ordaz y Perez (1999) , Ruiz-Garcia y Miranda (2003), Chopra (2005), entre otros. 
Lamentablemente en America del Sur no se dispone de suficientes registros sismicos, en 
funcion del tipo de suelo. 

Al no disponer de registros sismicos se decidio como, segundo trabajo de 
investigacion, generar registros ratificales en funcion del tipo de suelo y determinar la relacion 
entre el desplazamiento maximo inelastico con respecto al desplazamiento maximo elastico. 

Aguiar y Gonzalez (2006) generaron 28 acelerogramas sinteticos los mismos que son 
compatibles con los espectros del Codigo Ecuatoriano de la Construccion CEC-2000 para los 
cuatro perfiles de suelo. Se generaron registros con una aceleracion maxima del suelo igual al 
40% de la aceleracion de la gravedad, la duracion de los eventos fue variando desde 20 s., 
hasta 50 s. 



En la figura 4 se presentan los valores medios de la relacion A INE / S d encontrados por 
Aguiar y Gonzalez para los cuatro perfiles de suelo del CEC-2000. En forma general se puede 
indicar que el perfil SI corresponde a un suelo firme de gran resistencia a la compresion; por el 
contrario el perfil S4 corresponde a un suelo muy blando. El perfil S2 a un suelo de dureza 
intermedia y el S3 a un suelo blando. 



Para cada perfil de suelo se tienen cuatro curvas para ductilidades de 2, 3 y 4, y la 
curva de valores medios. En base a estas curvas se hallo la ecuacion de ajuste la misma que 
esta definida por: 



A - 1 + 



- + c 



A 4 j V 



n-i 






\ T sJ 



Donde los coeficientes “a”, “b”, “c” y “d” fueron obtenidos del ajuste para cada tipo de 
suelo y estan indicados en la tabla 4; T s es el perfodo caracterfstico del suelo que se 

considero igual al valor donde empieza la rama descendente del espectro del CEC-2000, los 
valores estan indicados en la ultima columna de la tabla 4. 
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0.00 0.50 1.00 1.50 2.00 2.50 

PERIODO (s) 

| NASSAR Y KRAWINKLER (1991) VALOR MEDIO -m- MIRANDA (2000) — Propuesta| 

Figura 3 Valores medios de j3 3 y ecuaciones de Miranda, Nassar y Krawinkler, 
y propuesta por Aguiar y Guerrero (2006) para a - 0 
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Figura 4 Valores medios de /? 3 encontrados por Aguiar y Gonzalez (2006) 



Tabla 4 Valores encontrados para hallar /? 3 . Aguiar y Gonzalez (2006) 



Perfil de Suelo 


a 


b 


c 


d 


T 

1 s 


SI 


30.00 


1.34 


-1.49 


0.60 


0.50 


S2 


71.80 


2.00 


-1.50 


0.50 


0.52 


S3 


81.04 


2.00 


-2.55 


0.50 


0.82 


S4 


86.00 


2.10 


-2.60 


0.48 


2.00 



4 REGISTROS SISMICOS 

Para el tercer trabajo desarrollado en el CEINCI se generaron 112 acelerogramas 
artificiales; 28 con aceleraciones maximas del suelo igual a 0.15 g; igual cantidad de 
acelerogramas se obtuvieron para 0.25 g; 0.30 g., y 0.40 g. Las aceleraciones maximas del 
suelo corresponden a las cuatro zonas sismicas del Ecuador. 
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De los 28 acelerogramas artificiales generados para cada una de las zonas sismicas, 
se obtuvieron 7 acelerogramas para cada uno de los cuatro tipos de suelo estipulados en el 
CEC-2000. 

Los datos mas relevantes de los sismos reales con los cuales se trabajo se indican en 
la tabla 5, se aprecia que son sismos de baja magnitud y con aceleraciones maximas del suelo 
menores a 0.1 g., por lo que se procedio a normalizar estos registros a las aceleraciones 
maximas de las cuatro zonas sismicas del Ecuador. Se destaca que todos registros indicados 
en la tabla 5 fueron obtenidos en la ciudad de Guayaquil. 



Tabla 5 Datos de los acelerogramas utilizados en el estudio. 



ESTACIO 

N 


PERFIL 

DE 

SUELO 


FECHA 


COMP. 


MAGNIT 

UD 


LATIT. 


LONG. 


PROFUNDIDAD 

(Km.) 


ACELERACION 
PICO (cm./s 2 ) 


DISTANCIA 

FOCAL 

(a.m.) 


U. 

CATOLICA 


SI 


11/06/91 


N-S 
E- W 


5,4 mb 


-3,1636 


-79,3767 


113,80 


4,274 


166,96 


U. 

CATOLICA 


SI 


30/05/90 


N-S 
E- W 


4,5 mb 


-2,2288 


-79,9963 


4,87 


3,401 


13,47 


U. 

CATOLICA 


SI 


21/11/91 


N-S 
E- W 


4,2 mb 


-2,3130 


-80,1520 


62,00 


3,693 


69,83 


U. 

CATOLICA 


SI 


21/11/91 


N-S 
E- W 


4,1 mb 


-1,9000 


-79,7900 


10,00 


1,816 


35,49 


U. 

CATOLICA 


SI 


25/04/96 


N-S 
E- W 


5,2 mb 


-9,7060 


-80,1890 


36,00 


8,218 


835,71 


BARRIO 

DEL 

CENTENA 

RIO 


S3 


26/04/93 


N-S 
E- W 


Registra 

do 


? 


? 


? 


10,994 




BARRIO 

DEL 

CENTENA 

RIO 


S3 


06/01/95 


N-S 
E- W 


5,4 mb 


-1,9005 


-78,0854 


180,58 


4,520 


271,60 


BARRIO 

DEL 

CENTENA 

RIO 


S3 


25/01/94 


N-S 
E- W 


4,8 mb 


-1 ,7022 


-78,2391 


198,64 


2,723 


275,63 


BARRIO 

DEL 

CENTENA 

RIO 


S3 


08/08/93 


N-S 
E- W 


? 


? 


? 


? 


2,939 




COL ING. 
CIVILES 


SI 


11/06/91 


N-S 
E- W 


5,4 mb 


-3,1636 


-79,3767 


113,80 


27,699 


166,96 



5 ESPECTROS ELASTICOS DE DESPLAZAMIENTOS 

Se han agrupado los registros sismicos por aceleracion maxima del suelo y por el perfil 
del suelo. Asi por ejemplo para los registros con aceleracion maxima del suelo igual a 0.15 g., y 
para un perfil de suelo SI, se tienen 12 acelerogramas reales y 7 acelerogramas sinteticos. 
Existe una diferencia marcada en la forma de los espectros de los sismos reales y de los 
sismos sinteticos como se aprecia en la figura 5. Los espectros de los sismos sinteticos se 
hallan agrupados en la parte superior y los espectros de los sismos reales en la parte inferior. 
Los espectros sinteticos estan identificados de Sdl a Sd7. 
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Sd (cm) 



ESPECTROS ELASTICOG DE 33UlOSREflLESYS»fTETICOS 



1.8CE+OI 

1.6QE+0I 

1.4CE+01 

1.ZE+01 

I.OCE+OI 

SOCE+GO 

&0CE+G0 

4.0CE+C0 

20CE+C0 

QOQE+CO 




Figura 5 Espectros elasticos de desplazamiento con los registros de 015 g., en suelo SI 



6 ESPECTROS INELASTICOS DE DESPLAZAMIENTOS 

En la figura 6 se presentan los espectros inelasticos de desplazamientos, para los 
registros con aceleraciones maximas del suelo igual a 0.15 g., sobre un perfil de suelo SI, para 
una ductilidad de 4. El comportamiento es similar al de la figura 5, en el sentido de que los 
espectros correspondientes a los sismos reales normalizados tienen valores bajos y estan 
agrupados. En cambio los espectros de los registros artificiales tienen valores altos y estan 
agrupados. Con Itnea gruesa se indica en la figura 6 la curva de valores promedios. 

Este comportamiento de que se agrupan los espectros de los registros artificiales y de 
los sismos normalizados, se observo en todos los grupos de datos. No es unicamente con los 
registros cuya aceleracion maxima del suelo es 0.15 g., sino con los de 0.25 g., 0.30 g., y 0.4 g. 

Los espectros inelasticos se obtuvieron para un modelo elasto perfectamente plastico, 
es decir para a- 0 . 



114 



Roberto Aguiar, Eduardo Aragon, Marcelo Romo 




( cm ) 



ESRECTROS INELASTICOS DE SISMOS REALES Y SINTETICOS 




— sdl 

— sd2 

— sd3 
sd4 

— sd5 

— sd6 

— sd7 
sdIO 
sdl 1 
sd12 
sd13 

— Sd14 

— sd15 
sd16 

— sd17 

— sd18 

— sd19 

— sd8 

— sd9 

PROMEDIO 



T(s) 



Figura 6 Espectros inelasticos de desplazamiento con los registros de 015 g., en suelo SI , para 

ductilidad de 4. 



7 RELACION ENTRE DESPLAZAMIENTOS NO LINEAL Y LINEAL 

Al tener una marcada diferencia en la forma de los espectros elasticos e inelasticos de 
los sismos reales y de los sismos sinteticos. Es conveniente no agrupar los datos y trabajar por 
separado la relacion A INE / S d . 

Con los espectros de respuesta elastica y con los espectros de respuesta inelasticos 
para ductilidades de 2, 3 y 4. Se procedio a encontrar la relacion entre el desplazamiento 
maximo inelastico con respecto al desplazamiento maximo elastico A INE /S d , para cada una 
de las aceleraciones maximas del suelo y para cada una de las ductilidades. Asf por ejemplo 
en la figura 7 se presenta esta relacion para los registros cuya aceleracion maxima del suelo es 
0.15 g., en un perfil de suelo SI y para ductilidad de 4 pero con los registros de los sismos 
artificiales. Con Ifnea mas gruesa se ha indicado la curva de valores medios, en la cual se 
aprecia que aproximadamente para valores menores a 0.6 s., la relacion A INE /S d tiene 
valores mayores a la unidad y que para cuando el perfodo tiende a cero esta relacion tiende a 
la ductilidad. La curva de valores medios con la que se trabaja es aquella en la cual la 
probabilidad de excedencia de los datos es del 50%. 

En la figura 8 se presenta la relacion A INE / S d para ductilidades de 2, 3 y 4. Para cada 
uno de los cuatro perfiles de suelo del CEC-2000, obtenidas con registros artificiales cuya 
aceleracion maxima del suelo en roca es 0.1 5 g. 
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A mE / S( j Ao=0.15 <j: Peifil SI; p.= 4 

4.0CE-tO0 
3.9CE-IOO 
3.0CE-tO0 
2.5CE-tOO 
2.0CE-tO0 
1 9CG-tOO 

i.oce-too 

5 OOE-OI 

o.oce-too 

OOOOOOOOOi-i-i-t-i-i-i-i-(MCM<MN<M<'l<MN 




— sdl 

— sd2 

— sd3 
sd4 

— sd5 

— sd6 

— sd7 

PRCMEDIO 



Us) 



Figura 7 Relacion A INE / S d para registros con aceleracion maxima del suelo 0.15 g., en perfil de 
suelo SI y para una ductilidad de 4. Obtenido con sismos ratificales. 



8 AJUSTE DE CURVAS 

Se detalla la forma como se obtuvo la curva de ajuste en base a los resultados 
encontrados para un perfil de suelo SI, para una ductilidad de 4 y para una aceleracion 
maxima del suelo de 0.15 g. Se destaca que se siguen los lineamientos generales propuestos 
por Chopra (2005). 

La figura 7 se presenta nuevamente en la figura 9 pero en el eje de las X, se coloca el 
logaritmo en base 10 de la relacion del perfodo T con respecto al perfodo caracterfstico del 
suelo T s y en el eje de las Y, se coloca el logaritmo en base 2 de la relacion A /NE / S d . Para 
el perfil de suelo SI se considero T s = 0.5 s. 



En la figura 9 con Ifnea un poco mas gruesa se han dibujado los ejes de coordenadas. 
Se desea ahora encontrar una curva que mejor se ajuste a los puntos indicados para lo cual 
como punto de partida se considera que esta tiene la siguiente forma: 



c *jm = JL = i + 

u 5, R„ 



f \ 
a 






f j \ c 



-i-i 



T 



En lugar de denominarle j3 3 se lo ha llamado C . Con esta acotacion se describe la 

ecuacion (10). Las dos primeras igualdades ya fueron indicadas restando manifestar que los 
coeficientes “a”, “b” y “c” se deben determinar del ajuste. Para los datos del suelo SI y /u-4 
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que se esta analizando los coeficientes son: a = 1.30; Z? = —1 .50; c = 1.76. Con estos 
valores en la figura 10 se indica la curva de valores medios y la curva que se obtiene del ajuste 
con la forma de la ecuacion (10 ). 



A IIIE / S(l 



35 

3 
25 

2 

15 

1 

05 

0 

^ llll 
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3 

25 

2 

1.5 
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Q5 1 15 



25 



A ,„ E / Sd S2 





T ( s ) 



T(s) 



CUCTIUDAD2 DUCTILIDAD3 



DUCT1UDAD4 



Figura 8 Relacion A INE / S d encontrada en perfiles de suelo SI a S4. 




Figura 9 Relacion logarftmica de la relacion A INE / S d en funcion de T !T S . Acelerogramas sinteticos. 
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Se observa en la figura 10 que la curva obtenida con la ecuacion ( 10 ) tiene un 
desfase hacia la izquierda con relacion a la curva de valores medios. Para corregir este 
problema se complementa la ecuacion de ajuste, quedando de la siguiente forma: 

r -i-l 



c„= 1 + 



X 



f \ 

a 






¥ 



V 



s y 



Donde X es un factor de traslacion de coordenadas y y/ es un factor de escala 
horizontal. En la figura 10 se aprecia que la rama horizontal inicial es mas grande que la rama 
horizontal de los valores medios para que estas dos rectas sean aproximadamente iguales se 
ha introducido el factor y/ . 




Figura 10 Valores medios y curva preliminar ajustada. Acelerogramas sinteticos. 



Para un perfil de suelo SI y para una ductilidad de 4. Los coeficientes y factores que se 
obtienen del ajuste son los siguientes: 



a = 1.30; & = -1.50; c = 1.76; /l = 0.03; ^ = 0.25 



En la figura 1 1 se muestra la curva de valores medios y la curva de ajuste obtenida con 
la ecuacion ( 11 ) para valores en los cuales la relacion A INE IS d es mayor a la unidad, se 
aprecia que con la inclusion de los factores X y y/ se tiene un muy buen ajuste. Para el lado 
derecho la tendencia del ajuste es hacia la unidad. 

Al trabajar con los acelerogramas reales normalizados a 0.15 g., que son las curvas 
que se vienen presentando se obtiene un mejor ajuste como se observa en la figura 12. Los 
coeficientes para este caso, son: 

a = 0.05; Z? = -12.50; c = 6.49; /l = 0.02; ^ = 0.70 

Se observo en el estudio que la relacion A INE / S d es la misma para el grupo de datos 
correspondientes a los sismos normalizados a 0.15 g.; 0.25 g.; 0.30 g.; y 0.40 g. Esto talvez se 
deba a que se trabajan con registros artificiales y con registros de sismos reales muy bajos que 
han sido normalizados. Como las relaciones son las mismas para las diferentes aceleraciones 
maximas del suelo, en lo posterior no se hace referenda a la aceleracion maxima del suelo en 
roca. 
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Figura 1 1 Valores medios y curva ajustada. Acelerogramas sinteticos. 




Figura 12 Valores medios y curva ajustada. Acelerogramas reales. 



9 RESULTADOS 

Se ha explicado con detalle la forma como se obtuvo la ecuacion de ajuste para el 
grupo de datos correspondiente a un perfil de suelo SI y para una ductilidad de 4. En forma 
similar se procedio con los restantes perfiles de suelo y ductilidades. 

El ajuste con la ecuacion (11 ) es valida siempre y cuando se tenga en el eje de las X, 
logaritmo de base 10 de la relacion T/T s y en el eje de las Y logaritmo en base 2 de la 
relacion A INE / S d . Cuando se trabaja con ejes normales, que no sean logarftmicos la ecuacion 
( 1 1 ) se convierte en: 
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Con esta aclaracion los resultados obtenidos se resumen en la tabla 6. Unicamente 
para los sismos sinteticos. 



Tabla 6 Valores obtenidos en el estudio para diferentes tipos de suelo y ductilidades. 



Ductilidad 2 


Perfil de Suelo 


a 


b 


C 


2 


V 


SI 


0.35 


-3.50 


1.40 


0.17 


0.87 


S2 


0.60 


-2.90 


1.31 


0.17 


0.82 


S3 


3.40 


-1.00 


1.50 


0.21 


1.00 


S4 


2.10 


-1.40 


1.00 


0.12 


1.70 


Ductilidad 3 


SI 


1.00 


-2.70 


1.40 


0.04 


0.15 


S2 


1.00 


-1.20 


1.40 


0.05 


0.49 


S3 


3.00 


-1.00 


1.80 


0.07 


0.73 


S4 


15.00 


-0.08 


1.40 


0.07 


0.30 


Ductilidad 4 


SI 


1.30 


-1.50 


1.76 


0.03 


0.25 


S2 


7.80 


1.00 


1.40 


0.02 


0.50 


S3 


1.30 


-0.20 


1.41 


0.01 


0.93 


S4 


0.23 


-0.60 


1.80 


0.04 


2.91 



En la figura 13 se presentan los valores medios encontrados en el estudio, que 
corresponde a la figura 8 y las curvas de ajuste que se hallan con la ecuacion ( 1 2 ) y los datos 
de la tabla 6, para un perfil de suelo SI. Se aprecia que existe una buena correlacion 
especialmente para periodos menores a 1.5 s. Entre 1.5 s., y 2.0 s., el ajuste no es tan bueno 
pero el criterio del ajuste es que las curvas tiendan a la unidad lo que implica el cumplimiento 
de la regia de igual desplazamiento. El comportamiento de los suelos S2 y S3, es similar a la 
del suelo SI . 

La figura 14 es similar a la figura 13 pero con los datos correspondientes al perfil de 
suelo S4. En este caso, para valores de perfodo comprendidos entre 2.0 y 2.7 s. se tiene que la 
relacion A INE /S d es menor que la unidad, para los restantes valores el ajuste es bastante 
satisfactorio. 



10 COMPARACION DE CURVAS DE AJUSTE 

Interesa comparar las curvas que se hallan con las ecuaciones propuestas por Aguiar y 
Guerrero; Aguiar y Gonzalez y las encontradas en este estudio para los diferentes perfiles de 
suelo. En la figura 15 se presenta esta comparacion para un suelo S2 y ductilidad 4. Se aprecia 
que los valores obtenidos con la ecuacion de Aguiar y Guerrero reporta valores muy bajos. Los 
valores que se hallan con las ecuaciones de Aguiar y Gonzalez son bastante parecidos a las 
encontradas en el presente estudio pero ahora la curva decae en forma mas uniforme en la 
zona de periodos cortos. 
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Figura 13 Comparacion de valores medios y curva de ajuste en suelo SI . 



a IHE / Sd Peifil S4 ; M-= 4 




CUCT1UD6D2 DJCTIUDAD3 DLCT1UDAD4 

AJUSTE 2 AJUSTE 3 AJUSTE 4 



Figura 14 Comparacion de valores medios y curva de ajuste en suelo S4. 

Para un perfil de suelo S4 y ductilidad 4 los valores que se obtienen con la propuesta 
de Aguiar y Gonzalez son bastante altos para periodos cortos, como se aprecia en la figura 16 
y no cumple con la Dinamica de Estructuras que establece que para cuando el periodo tiende a 
cero la relacion A , NE / S d debe ser igual a la ductilidad. En este sentido se aproxima mejor la 
ecuacion encontrada en el estudio. 
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PERFIL DE SUELO S2. DUCTILIDAD = 4 

A|„ E /Sd 




Figura 15 Comparacion de curvas de ajuste para perfil S2 y Ductilidad 4. 



PERFIL DE SUELO S4 DUCTILIDAD 4 




11 CURVAS DE AJUSTE EN FUNCION DEL TIPO DE SUELO. 

En la figura 17 se comparan las curvas que se obtienen para los cuatro tipos de suelo, 
para una ductilidad de 4, con la ecuacion (12) encontrada en el estudio. Se aprecia que las 
curvas para los suelos SI, S2 y S3 son bastante parecidas y que la curva para un perfil de 
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suelo S4 presents valores mas altos. Un comportamiento similar se observo para las otras 
ductilidades. 

El tener valores mas altos en un perfil de suelo S4, implica que los valores 
correspondientes al factor de reduccion por ductilidad R u son mas bajos. Por otra parte en 

todos los casos los valores hallados en el presente estudio son mas altos a los que se 
obtuvieron al no considerar el tipo de suelo, que corresponde al trabajo de Aguiar y Guerrero y 
menores a los que se hallan con las ecuaciones obtenidas por Aguiar y Gonzalez. 



DUCTILIDAD 4 




Figura 17 Variacion de la relacion A [NE / S d en funcion del tipo de suelo. 

Con el proposito de visualizar mejor los resultados, para la zona de perfodos bajos, en 
la figura 17 se presenta la relacion A INE IS d para perfodos menores a 0.5 s. Para perfodos 
mayores a esta cantidad se mantiene la tendencia de las curvas y para perfodos largos todas 
ellas tienden a la unidad. 



12 CONCLUSIONES 

Se han encontrado varias ecuaciones para determinar la relacion entre el 
desplazamiento maximo inelastico con respecto al desplazamiento maximo elastico, en un 
sistema de un grado de libertad, trabajando para ello con acelerogramas sinteticos, que son 
compatibles con los espectros de diseno del Codigo Ecuatoriano de la Construccion, para los 
cuatro perfiles de suelo estipulados en el Codigo y con acelerogramas reales normalizados. 
Del estudio realizado se desprenden las siguientes conclusiones: 

• Existio una marcada diferencias en la forma de los espectros elasticos e inelasticos 
obtenidos con los acelerogramas sinteticos y con los acelerogramas de sismos reales, 
por lo que se decidio trabajar unicamente con los registros sinteticos. Se descarto los 
sismos artificiales debido a que todos ellos fueron registrados en una ciudad. 

• No existe una diferencia notable entre los valores obtenidos para los diferentes tipos de 
suelo con las ecuaciones halladas por Aguiar y Gonzalez, con respecto a las 
encontradas en el presente estudio. Lo que si se debe reconocer es que para perfodos 
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que tienden a cero las ecuaciones halladas reportan valores que se aproximan a la 
ductilidad, de mejor manera a las que reportan las formulas de Aguiar y Gonzalez. 
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